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1}  ASSfRACt  rMAi4MT«wm  TOO  oKvee 

This  final  ropryi  covers  the  analysis  pertormed  under  Task  3  of  this  contract,  usaig  the  raw  station  mag^ 
nrtudes  collected  under  Task  t  35.21 1  staton  values  feom  e  maior  nuctear  test  sues  have  been  ted  mto 
»te  GlM  (General  Linear  Model]  inversion  »  dctcrnMne  joindy  the  cvera  sizes,  the  slabon  corrections,  and  the 
speofic  path  corrocfion  fgr  each  souco  station  paa  the  smHiHwieousiy-inlefTed  path  and  statKMi  corrections 
«re  related  to  known  geological  /  geophysical  features  Applying  tecsc  path  ^  station  corrections  to  the  raw 
station  magnitudes  of  any  individual  expioston  yields  a  systematic  reduction  in  the  ttuctuabonai  vanation  of  sta- 
«on  magnitudes  across  the  whole  network  with  a  reduction  factor  ran^ng  from  1  2  to  3  tor  an  Soviet  events  in 
our  set.  These  station  corrections  have  been  used  n  estimatsig  the  event  magnitudes  of  those  isolated 
explosions  which  were  not  included  n  the  GlM  inverston  With  Vi^e  pate-oomected  /  station-corrected 
^  ihnlP***)  rnj(Ljj)(NORSAR)  teas  between  the  southvrest  and  northeast  subregions  of  the 
Soviets  Balapan  test  site  is  assessed  as  0  06  magnitude  umt  |m  g  J  Magnitude-magnitude  regression  were 
contacted  wUh  the  standard  errors  in  the  event  magnitudes  toty  courted  tor  We  observed  some  incon¬ 
sistency  in  toe  scaling  relationships  The  systematic  teas  «  toe  L,  scale  can  be  removed  by  revising  toe 
absoMe  4  excrtatwi  level  tor  Novaya  Zemiya  eiptoskrts  That  is.  toe  ^  scale  suitable  for  Iranian  Rateau 
might  be  more  appropriate  tor  computing  of  Novaya  Zemlya  explosions 


IS  NUMRER  or  RACES 


»*  SUBitCr  TERMS 


Linear  Model. 

Maximum  itkefihood.  Rec  ession  Model.  Uncertamfy 


•  7  SICURlTt  ClASSmCAflON 
OR  RERORT 

Unclassified 


It  SECURITV  CIASSIRICATION 
OR  THIS  RAGE 

Unclassified 


It  SECURJTT  ClASS«<ATlON 
OR  ARSTRACT 

Unciassifiert 


1«  RRKf  CODE 


70  UMlTATlON  or  ABSTRACT 


VSN  iSdO  O'  iSO  SXOO 


ro''"  i’98  fRfv  ?  89) 

•  •rv-'RMMt  tH  4*»V  *.»d  /  fg  If 

;4«  -1} 


(THIS  PAGE  INTENTKXAiiY  LEFT  BLANK) 


ii 


SUMMARY 


Se^mic  yield  estimation,  which  n  ttased  on  the  seismic  magnitudes,  is  a  function  of  source  region 
effictency,  near-source  effects,  path  atlenuatiofv  wf  near-receiver  effects.  The  interplay  among 
these  factors  introduces  uncertairlies  in  the  esSmates  of  yields  of  underground  nuclear  explosions  using 
the  seismic  techrrique.  For  explosions  below  100  kt.  the  geophysicai  and  geoiogicai  corxfitions  of  ttte 
source  re^on  piay  a  critical  role. 

Teledyne  was  awarded  this  contract  in  Jitfiuary  1991  to  transcribe  arto  measure  the  (fgNal 
waveforms  from  six  mafor  nuctear  test  si**s.  The  objectives  of  the  project  are  to  investigaie  seismic 
magratude-yieid  regression  modeto  for  u  '^ophysical  and  geologtoal  regiorv  and  their  associated 
errors,  to  study  the  eartfportalion  of  thuw  re^  .jtion  models  between  nuctear  explosion  test  sttes  and 
estimation  of  assoo^ed  errors,  and  to  devcic  icaTtg  laws  that  aNow  yield  estknation  of  nuclear  explo¬ 
sions  emplaced  m  media  wdh  (ffferent  geophysicat  arto  geoiogicai  paramelars.  The  project  invotvas 
three  ma)or  tasks: 

(1|  Collection  «to  measurements  of  dgilal  seismograms  tor  waveforms  from  al  six  data  sets.  Due  to 
the  large  amount  ot  waveforms  to  be  oollecled  and  meesurtto.  Task  i  *'%  the  mm  labor-ittensive 
task  throughout  the  whole  contract  period  There  were  564  events  in  total  r'^easurad  /  processed 
The  measurements  arto  toe  station  magnitudes  of  al  1 1  phases  tor  al  evatif.  are  summeiized  in 
an  accompanytng  tinel  report  78E'4S17-1 1  TGAL‘94-C1  (Baumstark  arto  WugnBt.  1994). 

(21  Construction  ot  Geophysicai  arto  Geoiogicai  {G&G)  dtoabase.  The  results  derived  under  this  task 
are  described  m  a  seperaie  volume  of  the  Itoei  report 

|3|  Data  reduction  arto  interpretation  Thta  task  would  mctode 

(a)  the  construction  of  magniluda-magnaude  regression  models  (tor  al  test  sitos)  as  wal  as 
magrvtude-yieid  regression  models  (tor  NTS). 

(b)  investigaie  the  seeing  and  transportaMly  of  the  regressmn  models 

This  final  report  covers  the  unclassified  work  performed  ixtoer  Task  3  usmg  the  raw  stabon  magni- 
hides  collected  under  Task  t  No  aliampi  was  made  to  utSize  any  G&O  mtormation  obtained  urtoer  Task 
2.  due  to  the  classSied  nature  of  toat  task  Tha  report  Uso  desertoes  the  inversion  arto  regression  elgo- 
rtthms  used  m  tfvs  protect 

Two  magnsude  determinoaon  schemes.  OLM  (Generai  linear  Model)  and  EMILS  {Expacdorv 
Maximuation  *  Hertoive  Least  SqMres)  were  used  tor  clustered  arto  isoiated  nuctear  mds,  respectively. 
The  QlM  aigoriihm  not  only  produces  the  network-averaged  megnaudes.  which  are  optimal  in  the 
mextmum  likelihood  sense,  a  also  genertoes  two  types  of  by-products  the  stabon  oorrecions  and  lha 
path  corrections  This  QLM  roubna  has  been  eppsed  to  al  bme-domain  magnitudes  used  in  this  study, 
although  much  attenlon  was  paid  to  .  Another  norvstandard  elgorahm  we  used  in  this  study  is  a 
magnitude- magnaude  regression  twAm  which  accounts  lor  the  unoeruwtoes  in  the  evera  magnttudes 
fvrnuc^  ktonte-Carto  resampitog.  A  feature  of  automatic  oufier  refection  has  be»i  added  to  this  routtoe. 

The  huge  detabese  of  statton  values  based  on  short-period  veriici^<omponeni  (SPZ)  record¬ 
ings  has  been  fed  into  the  OLM  inversion  to  determine  simua^toousiy  the  event  sizes,  the  station 

W 


ooirtctiont,  and  spacMc  path  oorrection  lor  §dct\  sourca-station  pair.  Tha  linaar  tystam  Inwolvaa 
35,211  aqualiona  artd  3,306  unknowns.  Tho  simuRarteousiy-Marrad  path  and  station  oorractiona  aia 
raiatod  to  known  gaoiogicaVgaophysicai  faaturas.  Applying  thasa  path  and  station  oorractions  to  tha  raw 
station  magniludas  of  any  indMdual  axpiosion  ylalds  a  systamatic  raduction  in  tha  ftuduational  vartation 
of  station  magnihidas  acroaa  tho  whoia  nafwork  with  a  laducbon  fackx  ranging  from  1.2  to  3  lor  al 
Soviat  avanta  In  our  sat  Thasa  station  oorractions  hava  baan  usad  In  computing  tha  avart  magni- 
tudas  of  thoaa  isolatad  axploaiona  that  wara  not  indudsd  in  tha  QLM  invarsioa  With  thasa  pativ 
oorractod/statiorvoorractad  n\(P^ ,  tha  bias  balwaan  tha  southwast  and 

northaaat  subragiono  of  tha  Sowiat  s  Baiapan  tast  sita  is  assassad  as  0.07  magniiudo  uni  {m  u.),  which 
is  significartly  amaUar  ttian  mat  of  pravtous  studlas. 

Tha  raauts  praaaraad  in  this  raport  ara  vary  oonsistara  with  mosa  of  our  pravtous  study  (cf. 
Jih  af  a/. ,  1993)  using  assantialy  tha  sama  tachnologias  and  <Sftarani  explosion  data  sat.  In  oonrparlng 
n\.  Ms  .  and  n%lLf) ,  wa  obsarva  soma  inoonsistancy  in  dia  scaing  ralationsl^.  Tha  systamatic  bias 
in  tha  scale  can  be  ramovad  aUw  by  raisirrg  Nuttl's  path  attanuation  ooalftciarts  (which  wara 
dartvad  with  tha  coda  O  matfiod).  or  by  revising  tha  absokaa  Lg  excitation  Itval  lor  Novaya  Zarnlya 
sxplosiona.  Whichavar  approach  adoptsd.  there  would  be  soma  question  about  whether  me  Lg  scala  Is 
realty  aa  irvisportabia  as  Nusi  suggastsd  Tha  good  raauts  of  applying  his  NTS  Lg  -yiald  formula  to 
Samipaialinsk  axptosiona  may  turn  out  to  be  a  special  case  radtar  than  a  gsnaral  njia.  M  tha  Lg  acaia 
sutabla  tor  Iranian  Plateau  wara  usad  in  computing  tha  t%iLg)  of  Novaya  Zarnlya  axptosiona,  than  al 
tha  magnihida-yiald  ratattonshpa  lor  ms  three  phases.  Mg.  and  Lg  would  be  conslstart. 

To  tuiy  digasi  the  huge  aa^ioMan  daiasaf  would  raquira  a  signiKcais  amount  of  Urns  and  labor 
sdNch  are  certainly  beyond  ma  isval  of  adorn  saocaiad  lo  thts  Task  3  Tha  raauts  prasaniad  In  mia 
raport  should  be  ragardad  aa  praarrsrwry  In  nature  wth  many  topics  remain  to  be  axplorad.  This  rapoil  Is 
a  pitot  skxty  which  sarvas  to  itoistrata  how  ma  caraMy  moasurad  paramatars  (ooiactad  under  Task  i) 
can  be  usad  in  soma  magniuda-yiaid  analysia  uUkting  ma  stsiiaiical  software  package  davalopad  at 
Taisdyna 
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Caption  Pago 

ScaPof  piol  of  3  trnmmt  lypaa  of  Nation  n%ttot  Balapan  ovant  8M014  (JVE).  29 
Tha  26  good  tacofdlngi  and  l  dip  ara  ahoam  wth  Mad  dtciat  and  upwaid  ar- 
nm.  raapadNatjr.  Tha  raw  atation  i%  a  (lop)  hava  a  alandaid  daviadon  of  0.30 
m4i.  Appiying  ttta  ‘talmary*  atation  oowactiona  raduoaa  9ia  aeattar  to  0.21  m.u. 

Applying  tw  pmpoaad  “aaoondanr  oorractiona  to  oouri  lor  tha  path  aflacta 


raducaa  tha  aeattar  hjrthar  down  to  0.14  nui. 

Scattarplolof  3  dMaianttypaa  of  atation  <1%  a  tor  Balapan  avara  Ml 019.  90 

Scattar  plot  of  3  dHarani  typaa  of  atation  nv  a  lor  Oagalan  avari  710425.  31 

Scaaar  plot  of  3  dWarant  typaa  of  atation  nv  a  lor  Dagaian  avant  761230  32 

Scatter  plot  of  3  dMtarara  typaa  of  atation  ri\a  lor  Novaya  Zamlya  awara  33 

601014. 

Scatter  ptot  of  3  dtttorani  typaa  of  atation  /i%a  tor  Novaya  Zamlya  avert  94 
761020 

Map  of  ri%  atation  larma  Marred  Irom  a  QLM  Mvaraion  aolvaa  tor  3.306  unit-  42 


nowna  witti  96JI11  Inaar  oqtiationa  The  Ngh  oonatation  baiwaan  the  laclonic 
type  and  9m  atation  tarma  auggaatt  twi  ttwaa  an^Arteaf  oorractiona  do  rattaot 
the  (gpar  mania  conjuona  undamaath  tw  raoaivara.  Tha  4  dartianad  atara 
rapraaart  aoma  of  tw  laidaar  laai  attaa  uaad  In  taa  atudy. 

Tha  map  ahowtng  tw  Ipura  propagation  attacr  pop>  and  tw  complnad  attHon  43 
ampMcatlon  (botton^  dattrwd  aa  tw  aum  of  tw  tacalvar  larm  (Fgura  7)  and  tw 
path  attact  lor  NE  Baitawn  mtptomam. 

in«  iiVIfl  WOTig  WW  pfQptQVOn  iRMi  |10P|  §nO  WW  OOmDinOT  nBOn  44 

ampMcatlon  (bottong  dattnad  aa  tw  aum  of  tw  racaivar  larm  (Figura  7)  and  tw 
pat»  attact  lor  SW  flifapai  taai  atta 

The  map  ahowtig  tw  ipura  propagation  attacr  (tort  and  tw  oorrblnad  atatton  46 
ampMcatlon  (bottom)  dattrwd  aa  tw  aum  of  tw  racaivar  wrm  (Figura  7)  and  tw 

pati  attact  lor  TZ  Balapan  leal  atta 

in4  •noMTig  vit  picpigMon  mmcr  (lopi  ino  fi4  oonoviio  sanon  4o 

ampMcatlon  (bottom)  dattrwd  aa  tw  aum  of  tw  racaivar  larm  (Figura  7)  and  tw 
patt  alteci  lor  Oagalm  taat  atta 
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LM  Of  FlowM  (Continutd) 


Flgurt  No.  Caption  ^igo 

12  Tho  mop  showino  tio  *^0  pwpooation  olfoer  (lop)  and  fit  combinad  Nation  47 

ampMcadon  (bo&on^  doinad  as  tw  «an  o<  tho  meoNor  lorm  (figun  7)  and  tho 

potti  oCod  tor  Noflhom  Novaya  Zomlya. 

13  Tho  map  fhouNtoQ  tio  ‘turo  prepagahon  oHocT  (tort  and  too  oombinad  Nation  48 

ampNtoalton  (boaom)  dainad  at  too  Non  ol  too  mootoar  lomi  (ftpuro  7)  and  lha 

pato  oliaN  tor  Pa^Mo  Uoaa. 

14  Tho  map  ihoiMing  too  “}pun  propigNion  ottoo''  (top)  and  too  oombinod  NNton  48 

ampMcadon  (bottom}  dainod  aa  too  turn  o<  tho  loooMor  lorm  (Rguro  7)  and  too 

path  ottod  tor  NaMor  lloaa. 

15  Tho  map  ahowing  too  ‘Ipuro  propagation  oNocr  (lop)  and  too  oombinod  Nation  SO 

ampMcalton  (bottoro)  dattnod  aa  too  torn  oi  itm  looaixor  larm  (figuto  7)  and  too 

poto  ottod  tor  Yucca  FlaL 

*9  ¥9m  mWMBWlQ  WW  ^Uiw  prapSQBUOfi  •■Kl  (wpl  n  mm  OQfnPVlM  MflDOn  m 

ampMcalton  (bottom)  doSrtod  aa  too  aum  ol  too  mcoNor  lorm  (Rguro  7)  and  too 
path  ottod  tor  Franch  SaharN 

17  Tho  apNiai  poBom  rt  i%  Lf  roaiduatti  oi  CamtpNMnaii  oaptoatona  uatog  too  S7 

iftQ  (F*mm)  <tt  toia  atody  (Tabio  S)  and  miS  vatooa  raponod  by  NingdN  m  at 
(1982)  Tho  tocattona  art  btaod  on  Bocharoo  N  N  (1988)  and  Thutbor,  Quto. 

and  Richarda  (1999).  and  too  tododca  art  baaod  on  Bortoom  N  at  (1980)  and 
Loan  (1987)  Tho  potttm  auggooN  aomo  dWaranco  to  too  aouroo  mtdtom 
acroaa  tot  ChirvNi  lautt  aoparNwg  too  nortooaalom  and  aoutooaNam  pottton 
cttooiodaat  Tho  moan  o% btaa  btioton  sw  and  ne  Batapon  la  about 

0  08  m.u. 

18  Companion  ot  ttto(^B»)  agatod  miS  (NORSAR)  raponad  by  RtogdN  N  at  68 

(1882).  Tho  atopo  la  Itoad  N  1 
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I.  INTRODUCTION 


Th«  8cop«  of  this  project  Is  to  provide  (or  seismological  studies  to  investigate  seismic  yieid  scaling 
artd  transportation  oi  magnitude-yield  models  for  underground  nuclear  explosions  below  100  kilotons  (Id) 
between  different  test  sites.  The  objectives  of  the  project  are  to  investigate  seismic  magnitude-yield 
regression  models  for  distirrct  geophysical  and  geological  regions  and  their  associated  errors,  to  study 
the  transportation  of  these  regression  models  between  nuclear  explosion  test  sites  and  estimation  of 
associated  errors,  and  to  develop  scaling  laws  that  allow  yield  estimation  of  nuclear  explosions 
emplaced  in  media  with  different  geophysical  and  geological  parameters. 

Seismic  yield  estimation,  which  is  based  on  the  seismic  magnitudes,  is  a  function  of  source  region 
elastic  efficiency,  near-source  effects,  path  attenuation,  aivj  near-receiver  effects.  The  interplay  arrxing 
these  factors  IrSroduces  uncertairties  In  the  estimates  of  yields  of  underground  nuclear  explosions  using 
the  seemic  technique.  This  project  investigates  yield  estimation  for  explosions  below  100  kt,  for  which 
the  geophysical  and  geological  conditions  of  the  source  region  play  a  critical  role.  Path  attenuation  and 
near-receiver  effects  m  previous  studies  or  derived  from  this  project  are  used.  The  emphasis  of  this 
project  is  the  development  of  scaling  laws  that  reflect  source  media  with  different  geophysical  and  geo¬ 
logical  parameters. 

To  achieve  these  goals,  AFTAC/TTR  issued  a  contract  to  Teledyne  Geotech  Alexandria  Laboratory 
(TGAL]  In  January  '^Sl  to  transcribe  and  measure  the  digital  waveforms  from  six  major  data  sets: 
Novaya  Zemtya  (NZ),  Nevada  Test  Site  (NTS).  French  Sahara  (FS),  Balapan  (Shagan  River,  SR), 
Oegeien  Mourttakt  (DM),  arxt  Peaceful  Nuclear  Experimerts  (PNEs).  The  project  involves  three  major 
tasks 

(1)  Collection  and  measurements  of  cSgital  seismograms  (or  waveforms  from  all  six  data  sets. 

(2)  Construction  of  Geophysical  arxl  Geological  (G&G]  database  (for  all  test  sites  except  FS). 

(3)  Data  reduction  and  interpretation.  This  task  would  include 

(a)  the  construction  of  magnitude-magnlude  regression  models  (for  ait  test  sites)  as  well  as 

magniiude-yieid  regression  models  (for  NTS). 

(b)  Investigale  the  scafing  and  transportabUty  of  the  regression  models. 

Due  to  the  bvge  amount  of  waveforms  to  be  coSected  arxl  measured.  Task  1  was  the  rrxtst  labor- 
intensive  task  throuf^KJut  the  whole  contract  period  There  were  584  events  in  total  measured  /  pro¬ 
cessed  (Baumstark  and  Wagner,  1994): 


US  data  set 

246  events 

Shagan  River  dsAa  set 

97  events 

Oegeien  Mountain  set 

99  events 

.ovaya  Zemelya  data  set 

35  events 

PNE  (Me  set 

94  events 

French  Sahara  data  set 

13  events 

The  processed  dgitaf  waveforms  (converted  to  CSS  3.0  format),  instrument  responses,  as  well  as 
all  the  measurements  made  by  the  analysts  uskig  an  interactive  software  package  "geotool"  developed 


1 


by  Henson  and  Coyne  (1993)  were  delivered  to  AFTAC/TTR  in  six  8-mm  tapes  in  late  1993.  The  meas¬ 
urements  and  the  station  magnitudes  of  al  1 1  phases  tor  all  events  are  surrvnarized  in  an  accompany¬ 
ing  final  report  TBE-4617-VTGAL-94^1  (B^rrtsiartt  and  \NaQnK.  1994).  The  results  derived  under  Task 
2  are  descrtoed  In  a  separate  volume.  This  volume  of  the  final  report  summartzes  the  results  obtained 
under  Task  3  using  the  data  oolected  under  Task  1.  This  report  also  descrtoes  the  inversion  and  regres¬ 
sion  algorithms  used  In  this  prolecL 

Two  magnitude  determiruilon  echemes.  GLM  (Gerteral  Linear  Model]  and  EMILS  [Expectiorv 
Maximization  *  Iterative  Least  SQuares)  were  used  for  ctostered  and  isolated  nuclear  tests,  respectively 
A  criterion  was  tested  to  reject  those  stations  that  recorded  less  than  five  events.  As  a  result,  the  final 
number  of  evenu  led  to  the  GLM  inversion  was  slighlly  difterenl  than  the  number  of  events 
processed/measured  under  Task  1 .  The  station  oorrecttorM  derived  by  the  GLM  inverstons  were  appfied 
to  those  recordings  of  isolated  nudear  tests  before  the  station  magnitudes  are  ted  into  the  sirx)ie-event 
averaging  scheme,  EMILS.  Those  station  reoorcSrxjs  without  GLM  station  corrections  were  not  included 
In  the  EMILS  oomputalion.  Therefore,  the  number  of  everts  used  in  the  EMILS  computation  Is  also 
somewhat  less  than  the  total  number  of  events  actualy  measured 


Exploslone  Used  in  This  Study 


Pahute  Mesa,  NTS 

58 

RaMar  Mesa.  NTS 

22 

Yucca  Flat  NTS 

148 

Northeast  Balapan  (Shagan  River) 

30 

Southwest  Balapan  (Shagan  River) 

48 

Transition  Zone,  Balapan  (Shagan  River) 

19 

Degeien  Mourtain.  E.  Kazakh 

98 

Northern  Novaya  Zemiya 

30 

Ahaggar.  French  Sahara 

9 

(462  GLM  everts) 

Explosions  Not  Used  in  This  Study 

Amchika  Isiwid.  Aleutiane,  U  SA 

3 

CXitside  Nevada  Test  Sle.  U  SA 

6 

Southern  Novaya  Zemiya 

4 

PNEs,  U  S  S  R. 

94 

(107  EMILS  everts) 

This  report  oontaine  three  parts.  Sections  R  wid  Rl  descrtoe  the  aigorfihms  we  used  In  oomputirx) 
the  maxImum-BreVtood  evert  magnliudes.  The  krverston  resuRs  are  strtmrartzed  In  Section  IV.  In  Sec¬ 
tion  V,  we  present  some  preUminary  analysis  and  interpretation  of  the  resuts. 


i.  STATION  MAGNITUOE  COIIPUTATION 

■.1  FofimilM  UMd  In  Station  Magnltuda  ComfxAailon 

Tha  maaniramanla  kicbda  tha  dbpUcamam  ampttuda  In  nm.  tha  partod  in  aaoonda.  tha  phna 
nwna  (#.9. .  f%.  kl$,  L^.  PS),  lha  MonnNion  which  indudM  ttw  apicartar  and  tha  avanl  nama. 
Each  phata  haa  a  tpadHc  tomula  lor  datarmining  tha  magnibjda.  and  hanca  dWarant  paramatara  miohl 
baracMrad.  Tha  lonnuiaa  usad  In  Ms  study  art  daachbad  brially  In  tha  following; 

(1)  n\  m  log(A^  ♦  B(A)  lor  20*  <  A  <  05*.  wh^  B(A)  is  lha  dManca  normaizar  darivad  by  Vafth  and 
Clawaon  (1072). 

|2)  /hk  (P« )  ■  togCA)  ♦  2.42  100(A)  •  3.96  lor  A  <  10*  {et.  Vargino  and  Mansino.  1990). 

(3)  For  ,  two  dMarant  lornsilaa  ara  utad; 

M  A  >  25*  .  ■  loo(A^  ♦  1.M  100(A)  ♦  3.30  {fit.  Quiarbaio.  1945). 

N  10*  <  A  <  25*.  Ms  •  IO0(A/T)  ♦  1.07  loo(A}  ♦  4.16  {ct.  NuOl  «td  Kbn.  1975). 

Al  Ihraa  tono-padod  oomponanls  wart  lolalad  usino  lha  undassiitd  location  batora  an  LO  rras 
pidiad.  M  tha  roialad  wavatomv  axh/bt  a  station  10  avam  azinsith  which  is  sionificanily  dMarem 
from  tha  pradktad  atinsjih.  Vtan  tha  pickad  LR  anchor  LO  ara  traalad  as  noisy  raoordhtos. 

(4)  For  n^{Lf) .  Jh  and  Lyrmas  (1993)  suppast  tha  lolowino  formuta: 

Athouoh  tNa  r\{Lf)  formula  migN  appaar  to  ba  dMarant  fcom  moat  othar  lonnulaa  ourrandy  In 
uaa,  this  aquation  la  actualy  aquivalant  10  NuttTs  (I986ab.  1987).  For  inaumoa.  a  aaismic  touroa 
inaastamU.S.  wOhl-tacLs  amplluda  of  110  lunai  10  km  apicartraldbtanca  would  corraspond 
10  a  m^(Lf)  of  4  0272  ♦  2  0414  ♦  0  3333  •  1  4019  ♦  0.0000  •  5.000,  tfta  sama  valua  that  NuM  t 
original  2  stap  formiaa  would  giva.  Tha  ^  and  q  vakias  utad  in  Ms  proiacl  ara  Wad  In  Tabias 
1  through  3.  ankrala  ara  pickad  batwaan  0*  and  40*.  Two  dMarant  valodly  windows  ara  usad 
MpicMnotha  anrkrals: 

(a)  4.0  knVsac  10  2.5  knvsac  (lor  a  <  10*), 

(b)  4.0  knwsac  to  3.0  knViac  (tor  10*  <  A  <  40*}. 


iL2  Tlm^Hkwnaln  0«t«nnlfUMion  of  Ptf  h  Conoctlon 

Thor*  art  ttvtrai  dWtrtnt  approachtt  that  havt  bttn  proposed  lo  determine  the  p^h  Ooh|; 

(A|  Apply  the  ooda-0  method  of  Herrmann  (1960).  as  dd  NuttI  (1986ab.  1987, 1988). 

(q  Synlheabe  the  atong  itm  great-cirde  pdh  between  the  source  and  the  receiver  using 

the  Z-dbnenaionai  (Vh  map  of  that  region  (Jit  and  Lynnes,  1993). 

(C)  Apply  OLAI  (General  Linear  Modal]  or  LSMF  (Least  square  Matrix  Fadorizabon]  inversion  to  Mer 
ttte  path  oorrectiona  along  w«h  the  source  terms  (Jh.  1992). 

Approach  (C|  would  perform  very  wel  when  some  extra  relabia  Information  about  the  evwts  (a.g. . 
die  average  of  or  values)  is  available  to  constrain  the  join  inversion  (Jti.  1992).  Here  we 

provide  another  approach  which  is  very  simlar  to  (C]  except  that  die  stations  are  caKirated  irnSvidUaly 
wlh  those  evenis  tor  which  Niddi  (1988)  already  determined  the  values. 


In  processing  the  data  set  asserrbied  under  this  contract,  the  “sustained  maximum  motion“tit 
phase  is  measured  h  a  manner  identicai  to  that  wNch  Nutdi  (I986ab.  1987, 1988)  proposed.  That 
is.  die  ampMuda  equaiad  or  exceeded  by  die  three  targest  ampSude  waves,  of  the  verticaFoomponeni 
Lf  waves  with  period  around  1  second  was  picked.  The  station  ampiaude  reading  is  first  corrected  for 
the  effects  of  geometrical  sprearfing  and  dfopersfon  wlh  the  formula  appropriate  lor  the  Airy  phase:  the 
residual  (relative  to  Nultrs  r%{Lf))  la  then  regarded  as  completely  due  to  the  anelastic  attenuation 
along  the  path: 


lll.i(knvdeg) 


where  A  is  the  epicenirai  dislanos  in  km.  A(a)  is  the  observed  L,  ampfituds  measured  In  the  time 
domain  in  pm  (miciona)  at  the  epicenirai  dstance  of  A  km.  The  oorrespondfog  0(f)  can  then  be  deter* 
mined  in  a  straightforward  manner 


o(n- 


a  f 

T  U  ' 


PI 


where  U  is  the  group  vefodly.  Once  a  suie  of  0(f)  vakies  is  avalabie  for  a  station  of  MeresL  a  Inear 
regression  is  dien  conducted  to  find  the  maximum-Bielhood  estimate  of  the  quafity  factor.  Oo.  as  wel  as 
the  trequency-dapendsncy.  i|.  via  the  model: 


Q(f)  -  Oo  n .  or 


togCQWl-»og(Oo|  +  ii  t. 


W 


The  maximum>BiaNhood  regression  scheme  tor  the  censored  case  is  described  in  Appendx  B. 
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N4  Path  Corrtcllont  for  Novayt  Zomlya,  S«m4)alatinsk,  and  NTS 

Tha  resulino  path  corrections  tor  Novaya  Zemiya  test  site  are  listed  in  Table  1,  along  with  those 
corrections  of  NuttU's  (1988).  The  match  is  (airly  good.  This  simply  suggests  that  Teiedyne’s  arrpl- 
lude  measurements  lumished  by  Baumstark  and  Wagner  (1994)  are  very  oorrsistent  with  NuttM's.  k  is 
Interesting  to  note  that  1ST  (Istanbul,  Turkey)  and  TRI  (Trieste,  Italy)  dkt  record  Lg  phases  from  large 
hialoricai  Novaya  Zemiya  events.  Along  with  the  7  WWSSN  stations  lor  which  NuttH  (1988)  already  pub> 
■shed  the  Oo  values,  now  we  have  a  total  ot  12  paths  cailirated  lor  Lg  waves  from  Novaya  Zemiya. 
Statione  KON  (Konsberg,  Norway)  and  KBS  (Kingsbay,  Svatt>ard)  are  not  wel  corrstrained  due  to  the 
Imited  data  sbe,  and  hence  Nuttrs  (1988)  Ooh|  would  have  to  be  retained 


ML  IIAXmUII-UKEUHOOO  NETWOlUC  MAGNrTUDE  ESTIMATORS 

NL1  Slngl«^«nt  MaxiimnHliMlhoocl  Eatbnator 

Ttw  problwn  at  Mtimating  body-wa/«  magnMude*  (/n^ )  using  amplitudes  read  at  a  number  of 
leooRing  staMone  is  trequently  oonplicated  by  the  fact  that  the  data  may  be  heaviy  censored.  This 
arises  either  because  of  dpping.  where  al  ampitudes  can  be  determined  only  to  exceed  a  given  lower 
bound  (/.a.  ttw  right-censored  case  in  statistical  terms),  or  because  the  signals  are  weaker  than  the 
ambient  noise  level  and  hence  are  not  detected  (L0.  the  left-censored  case).  If  one  simply  averages  tfie 
magnftudes  tor  those  staftons  which  detected  an  evert,  without  regard  tor  those  that  dpped  or  dU  not 
record,  serious  biases  may  resul  in  the  evert  magnftude  estimated. 

For  single-evert  network  n\  determination,  at  least  three  types  of  station  magrtiude  ought  to  be 
considered: 

[0]  the  station  magnitude,  X.  is  Imown  as  Xq. 

[1]  X  is  only  known  to  be  less  than  certain  level,  say.  ti, 

(2)  X  is  only  known  to  be  larger  than  certain  level,  say.  Ij. 

We  assume  that  the  observed  station  magnitude.  K  can  be  represerted  as  the  sum  of  the  unk¬ 
nown  evert  magnitudo,  p.  and  a  perturbing  random  noise,  v. 

X-p-i^v  IS) 

where  v  Is  assttfned  to  be  a  Gaussian  random  variable  wlh  mean  zero  and  starKlard  deviation  o. 
Elegart  maxlmunv-llieiftwod  theory  can  be  derived  tor  this  Inear  model.  S<a)pose  there  are  no.  n,.  and 
n,  stalton  reoordtogs  for  each  type,  respectively.  The  condRionai  Btelhood  function  of  the  censored 
observations  (  Xo,  ti.  ti)  given  the  network  magnitude  p  and  o  is 

L(Xo.t,.ijip.o)-np(x,«xqip,o)-np(Xj<t,jip.o)-np(><j>*a'tLo) .  m 

1-1  j»i  1-’ 

and  toe  log-Btettwod  fiawtion  is 

lnL(Xo.l,.talp.o)--^ln(2i»*)---~  J(xc-p)»t^  £ln®(z,^-*^  £ln4K-Za)  .  U) 

where  Zt  ■  (ti  •  p  )/o  :  Xq,  t,.  and  to  are  ootocttons  of  the  observed  station  magnftudes  of  each  type, 
respect^eiy.  and 

4(u)  -  ^exp(^) .  dKu)  ■  £a(x)dx  .  P) 

are  the  probafaWy  densfty  furwtion  and  probabSty  dtotrtwtion  function,  respectively,  of  the  standard  nor¬ 
mal  random  variable. 

Solving  ■  0  imples  that  d,  the  optimal  estimate  of  o,  must  satisfy  the  foBowIng  necessary 

OO 

condtoon: 


0 


SoMng  ■  0  hnpfes  that  A.  tha  optimal  tsUmato  of  p.  must  satisfy  tha  folowHng  naoassary 
oondWon: 


H  H 


AtkMng  (n,  -f  nj)  A  to  botti  sldaa  of  (101,  ^  than  dMdng  both  skies  by  (r^  -f  n,  -f  yiakls 


[101 


h#  ti^  + 1^  ^ 


fill 


Tha  righl-^iand  skJa  of  Equation  (ill  happans  to  ba  tha  sampla  maan  of  “ar  data  with  tha  cansorad 
maasuramants  rapiacad  by  thair  oofraspondng  bast  fiO-ln  (saa  Appandbi  A): 


ir— n;-(lEIXlX.x^l  +  £E(XlX<t,,l  +  2E(XlX>tal).  (121 

"0  ♦  "1  ♦  •a  pi  pi  pi 

Conaa<|jantly,  wRhin  tha  oontaxt  of  Gaussian  asstanption.  ona  can  translata  thosa  seamingiy  not-that* 
pradsa  statamants  of  X  >  t  or  X  <  t  Mo  quantitativa  oonstraMs  which  can  oouple  wlh  other  measura- 
manta  of  type  0  easily.  Thus  Equations  (ill  vxt  (121  ptovlda  tha  theoraticai  (ustification  of  an  fearidon 
prooadura  to  ba  tflscussad  below. 


An  itarniva  procadura  calad  "EM  atgodthm"  (EspactatiorvMaximttation  algorthm]  (Dempster 
ar  at .  1977)  can  ba  applad  to  solva  for  p  and  o  in  a  vary  straightfotwatd  manner.  To  start  tha  laration. 
ona  needs  an  MUai  guess  of  o  and  p.  A  good  Miiai  vaiua  of  p  is  tha  sampla  maan  of  Ml  typa^>  station 
magniiudss.  Since  buMin  typically  axhMs  a  o  (of  single  obsatvMion)  around  0.3  magnitude  uni. 
this  vaMa  can  serve  as  tha  inWai  value  of  o.  Tha  laration  procadura  folows: 

(1}  Based  on  tha  current  astimatas  of  p  and  o.  tsfiece  al  tha  cansorad  data  wlh  their  oorra^rondbig 
ooncPbofUti  expectations  (cf.  tha  ligN-hand  side  of  Equation  (11]}.  This  Is  tha  so-calad  "E  step"  of 
tha  EM  aigonhm. 

(2]  Compute  A  as  tha  sampla  mean  of  these  "laflnad  observations.'* 

(31  Update  tha  astimala  of  o  using  EquaHons  (91. 

(4]  Repeal  (1H31  untl  soma  oonvarganoa  criarion  ia  maL 

Steps  PI  and  (31  oorratiuta  tha  "U  step"  of  tha  EM  algoitthm.  Note  that  in  tha  norvcanaoring 
case,  f.a. .  ni  >  0}  •  0,  A  MXf  ^  would  reduce  to  the  regular  sample  maan  and  the  RMS  residual,  respec¬ 
tively: 


C101 


2>«, 

H _ 


rV) 


.o* 


tl _ 


no 


Tlw  algoithm  dMCftMd  abovt  has  baan  impiamartad  as  a  ulitty  program  “amfe*  vvhtch  axpacts 
to  raad  )ust  tm  ootumns  of  data  raprasaniing  tha  data  typa  or  '*>")  and  tha  actutf  data.  Taka 

Novaya  Zamiya  avani  NZ661027  (Oclobar  27,  1966)  as  an  axampla.  Tabia  4  Ists  tha  stMion  r%(Lf) 
vahMS  of  this  Novaya  Zamiya  avant  basad  on  our  formula  as  wal  as  tha  patfi  oonractions  wa 

Instalad 


Tabia  4.  Station  Racordings  of  Novaya  Zamfya  Explosion  661027 


Station 

A* 

Ampituda  (rvT^ 

Pariod|sac] 

mm 

q 

Vaiocify 

ntk(4) 

COP 

24.57 

870S 

131 

668 

0.41 

33 

6.341 

KEV 

9.48 

<1833.6 

038 

249 

0.74 

3.7 

<6306 

NUR 

867.5 

1.08 

433 

0.42 

3.6 

6.389 

8TU 

31.67 

234.3 

130 

550 

0.55 

33 

6.514 

UME 

15.58 

1168.2 

120 

397 

0.82 

33 

6.525 

ESK 

29.23 

155.7 

1.68 

3.6 

6.423 

1ST 

34.70 

49.5 

0.93 

561 

0.64 

3.6 

6.464 

KON 

21.91 

7893 

122 

496 

EES 

3.6 

6.518 

TRI 

33.38 

163.1 

2.09 

417 

0.24 

3.6 

6.221 

Thara  ara  8  good  slgn^  and  1  noisy  maasurament: 

VSSOS 

*w’aias 

•W-SL5M 

‘^’sjas 

•w-asaa 

V$.404 

vssfs 

“m-9221 

N  tha  cansorad  laooiding  of  6.506  at  tha  station  KEV  is  ignorad.  tha  avam  magnkuda  would  ba 
6.424(10.037.  Tha  program  'lamas"  givas  tha  maximunv4kaNiood  astimata  as  6.42010.034.  using  al  9 
obaarvattons.  Basicaly,  what  tha  maximunHkelhood  mathod  doas  is  to  utBza  ttM  cansorad  Mormation 
0*  nV(^)(KEV)  <  6.506  as  axtra  constraint  to  rafbta  tha  Manad  porametar  obtained  wlh  tha  stan¬ 
dard  least  squares.  For  this  avert,  NuttH  (1988)  gave  a  n^{L^)  of  6.45.  Oitf  QLM  nib(^)  of  overt 
Is  6.4210.03  (cf.  Table  5). 
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lil^  SUnuIttntous  Inversion  of  Evertt  htegnltudes.  Station  Terms,  and  Patli  Tenna 

M^.1  General  Concepts  of  Joint  biverslon  Model 

As  descrt)ed  in  Section  II,  the  conventional  defirvtion  of  the  station  magntuds  is  computed  as 

mb  ■  IO0io(Ani  +  B{A)  .  (141 

where  A  is  the  dlqjiacemem  ampMude  (in  nm)  and  T  Is  the  predominart  period  (in  sec)  of  the  F  wave. 
The  B(4)  is  the  dtatance'correction  temn  that  oompensales  tor  the  change  of  P-wave  ampOtudes  wfth 
dbtance  (e.g. ,  Gutenberg  and  Richter,  1956:  Veith  and  Clawson.  1972).  n\  in  (14]  is  also  denoted  as 
m,  in  Marshal  ef  ad  (1979).  The  ISC  (IrSemational  Seismograph  Centre]  buletin  riv  is  Just  the  network 
average  of  these  raw  station  rr%  values  wlhoul  any  further  adjustment.  Implicitly,  the  ISC  buletin 
assumes  a  Inear  model  as  the  tolowing: 


(D  «  E  ♦  v(j}  ,  (15] 

where  (D  the  stabon  magnitude  recorded  at  the  station  ]  lor  the  event  of  size  E.  and  v  is  the  rar>> 
dom  perturbing  term. 

Now  cortsider  Nc  explosions  detonated  at  f%  source  regions  that  are  recorded  at  some  or  al  of  14$ 
stations.  In  LSMF  (Least  squares  Matrix  Factorization)  and  the  standard  GLM  (General  Linear  Model] 
schemes  (e.g. ,  Douglas.  1966:  Blandlord  and  ShunrMray,  1962:  Marshal  ef  ad .  1964;  Uwal  ef  ad , 
1966;  Jih  and  Shumway,  1969;  Murphy  ef  ad ,  1989).  K  is  assumed  that  the  observed  station  (ij)  is 
the  Stan  of  the  true  source  size  of  the  i-th  event.  E(i).  the  receiver  term  of  the  Hh  station.  $(]),  and  the 
random  noise,  v(ij}; 


r?%  04)  -  E(l)  +  S(D  ♦  v(IJ)  .  (16] 

The  receiver  tenn,  S(D,  Is  constant  with  respect  to  al  explosions  from  tflierent  test  sles,  and  hence  I 
would  inherently  rebect  the  "averagerr  receiver  effect  —  provided  the  paths  reaching  the  station  have 
broad  azimuthai  coverage.  When  worldwide  explosiorts  are  used,  ttie  standard  deviation  (o)  of  the 
noise  V  in  (16]  is  typictfy  around  04  m.u. 

If  LSMF  or  GLM  is  appled  to  everds  wlhin  a  smaller  area  of  source  region,  ttten  the  o  of  v  hi  (13] 
could  reduce  to  O.lS-04  m.u.  However,  the  resul  of  such  "singie-test-sle  (3Uyr  approach  should  be 
Interpreted  or  uMzed  cautiously.  The  event  n\  values  (de..  the  "E'*  tonn  hi  (lOD  so  determined  are 
exceleni  esbirades  of  the  ”relative  eouroe  size"  for  that  lest  sle  only.  I  this  '‘single<test-sle  GLM" 
Inversion  is  appled  to  several  test  sites  separately.  I  may  not  be  easy  or  obvious  to  find  a  oonsistem 
baselne  for  estimating  the  "absolito  yield.”  since  the  recortlng  network  is  typicaly  dMerenl  from  one 
test  sle  to  another,  and  hence  the  station  terms  are  inevlably  hioonsistent.  Fiathermore,  the  station 
terms  derived  by  the  "shigle-lest-sile  GLM"  may  not  necessarfy  represent  the  aOenuation  undemeatfi 
the  receiver  side  alone.  They  could  be  "oontamlnrted"  or  eomethnes  even  overwhebned  by  the 
pMhfnear-souroe  effects  shared  by  the  explosions  confined  hi  a  narrow  azimuthai  range.  This  could 
explain  the  once  puzzing  and  controversial  phenomenon  Butler  and  Ruff  (1960)  (also  Butler.  1981;  Bur- 
dkk,  1981)  reported,  namely  that  using  Soviet  explosions  from  one  test  sto  Mone  may  fal  to  (Sscem  the 
attenuation  dHferential  between  the  eastern  and  western  U.S.  There  is  no  doubL  however,  that  the  GLM 
or  LSMF  type  of  methodology  can  infer  the  station  terms  which  are  strongly  oorrelaled  with  the  upper 
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mantle  attenuation  underneath  the  stations,  provided  the  seismic  sources  have  a  broad  spatial  coverage 
as  did  those  in  North  (t977),  Douglas  and  Marshall  (1983),  Litwal  and  Neary  (1985).  Ringdal  (1986),  Jih 
and  Wagner  (1991),  and  many  others.  The  even  magnitude  derived  with  Equation  (3)  is  hereby  denoted 
as  In  Marshal  H  at.  (1979).  a  priori  inlormation  about  the  velocity  underneath  each  station  is 
used  to  determine  Is  associated  ‘‘determmistic''  receiver  correction,  S(D.  and  the  network-averaged 
ma(^)itude  based  on  the  st^ion-corrected  magnitudes  is  caled  m, .  The  receiver  oorrectiorrs  as  derived 
in  Equation  [16],  however,  are  Merred  Jointly  from  a  sule  of  even-station  pain,  and  no  a  priori  geophy¬ 
sical  or  geologicai  oorkilion  Is  assunred  (and  herrce  the  dlteren  notation  m;^).  The  high  correlation 
between  the  tenonic  type  and  the  GLM  station  terms  suggests  that  the  empirical  station  corrections  do 
reflect  the  averaged  upper  manie  oorxtlions  urrdemeath  the  receiven,  I  the  azimuthal  coverage  at  each 
station  is  broad  enough. 

Jih  arkf  Wagner  (1992ab)  and  Jih  ef  at  (1993)  propose  to  reformulate  the  whole  model  [16]  as 

X(ii)  »  E(l)  ♦  S(0  ♦  F(kj)  ♦  v(i  j)  .  [17] 

where  F(k.j)  is  the  correction  term  at  the  i-th  station  lor  the  propagation  effect  or  the  near-source 
focusin^defocusing  effect,  which  is  constan  for  al  evens  (induing  this  Ur  even)  in  the  k-th  “geologi¬ 
cally  arxf  geophysicatly  uniform  regioa“  For  each  seismic  station,  this  F  can  be  regarded  as  Is  azimu¬ 
thal  variation  around  the  mean  station  term  S.  However,  as  explained  previously,  I  would  be  more 
appropriate  to  consider  F  the  path  or  near-souroe  term  because  the  back  azlmiths  at  the  station  could 
be  nearly  Idenical  lor  ac^acen  test  sites  (such  as  Oegelen  and  Murzh8(),  and  yet  the  “F“  temns  could  be 
very  differen.  By  incorporating  the  F  term  Mo  the  model,  the  o  for  world-wide  explosiorts  Is  reduced  to 
about  0.2.  roughly  the  same  level  that  which  a  “sinf^test-sle  0114“  could  achieve,  intuitively,  the 
presen  scheme  (Equation  [17]}  provides  a  more  detafled  (and  hence  better)  model  than  that  of  Equation 
(16)  In  desertring  the  whole  propagation  path  from  the  source  towards  the  receiver.  Simply  put.  Equa¬ 
tion  [16]  yields  a  strorxjer  fluctuation  in  the  source  terms.  E.  as  wel  as  a  larger  standard  deviation  of  v 
because  each  term  in  the  rign-hand  side  of  Equation  [16]  would  have  to  “absorb"  part  of  the  missing  F 
term.  The  resuUng  new  even  magnitude  (viz. .  E(i)  in  [17])  is  hereby  called  nio  to  avoid  confusion  with 
the  mj  defined  In  Marshall  ef  aL  (1979)  that  corrects  for  the  source-region  attenuation  and  station 
terms  solely  based  on  pubSshed  velocihr 

Roughly  speaking,  the  model  deserttred  in  [17]  has  the  loBowing  advantages; 

•  It  provides  more  stable  measuremens  across  the  whole  recording  network,  as  compared  to 
the  conveniiortal  GLM  or  LSMF  procedure  which  only  corrects  lor  the  station  terms.  The  reduction 
in  the  standard  deviation  of  network  from  m,  to  fflo  could  reach  a  factor  of  nearly  3.  As  a 
result,  the  scatter  in  mg  versus  log(yield)  is  smafler  than  that  for  (Xher  n\ . 

•  The  separation  of  the  path  effect  from  the  sttf  ion  effect  is  a  crudai  step  to  investigate  the  various 
propagdion  phenomena,  which  in  turn  would  Improve  our  understanding  of  the  seismic  soiace  as 

wel. 

We  have  applied  this  model  to  462  worldwide  explosions,  and  the  resuHng  mg  values  of  these 
explosions  are  Isted  in  Table  5.  The  194  stations  are  selected  such  that  each  dation  records  5  or  more 
good  explosion  signals.  Most  of  the  earSer  magnitucle-ylk.^  study  at  Tetedyne  has  been  relying  primarily 
on  WWSSN  as  the  core  recordkig  network  (Biandford  ef  aL.  1963:  a.g.,  Jii  ef  1993:  and  rmny 
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others),  IMS  s«(  contains  signiticari  traction  ol  data  conirtiutsd  by  USAEDS  limited  States  Atomic 
Energy  Detection  System)  seismic  network  (Saumstark  and  Wagner,  1994).  In  this  data  set.  there  are 
^,890  signate,  1.169  noise  measurements,  and  7.152  dipped  measurements  from  10  test  sites  that  are 
used  to  invert  for  the  3306  unknovm  paramMers  wSh  the  maiimurrvtiielihood  approach.  The  standard 
deviation  ol  v()J)  in  (14)  is  0  168.  as  compared  to  that  of  0.240  i  the  oonventiond  GLA4  (Equation  (13))  is 
appMed  to  the  same  data  set.  The  algorithm  and  sample  input  fles  are  descrt>ed  in  the  rtext  sectbtt 

18.2  J  QLM  Iteration  Procedure 

Equation  (17)  is  a  special  case  of  the  general  linear  models  (GLil).  An  keralive  procedure  based 
on  the  EM  algorithm  is  presented  below  The  basic  ideas  are  very  similar  to  those  underlying  the  sin(^ 
event  network  averaging  presented  in  Section  lil  t. 

StepO 

Set  up  initial  conditiorts  as  follows; 

(1)  o  «  0.3  magnituds  uni, 

(2) S(j)  •  0  lor  j-  1.2 . N,. 

(3)  F(kJ}  -  0  lor  )  .  1.  2 . and  k  .  1.  2 . Nf. 

Step  1 

Compute  everk  magrvtudes.  E(i),  lor  I  ■  i . Nc  as 

E(0  •  -  S(i)  -  F(kj)) . 

where  •())  Is  the  number  of  stations  that  ‘Teoordetr  the  everk  1. 

Slep2 

Compute  dation  correctiorw.  SO),  lor )  •  t . Ng  as 

SU)  -  • 

where  f(i)  is  the  number  of  events  "recorded"  at  station 
Step3 

Compute  path  corrections.  F(k.i).  for )  •  t . Ms;  k  •  t . Mr  as 

■ 

where  f(  (kj) )  is  the  number  cl  paths  from  the  test  ska  k  (where  the  event  i  is  locaied)  to  the  sta¬ 
tion  ).  Thie  step  is  skipped  if  options  0  through  2  are  chosen.  Conseqiently,  F(k.  0  wfl  remain  0 
lor  al  k  and )  when  is  the  desired  event  magniude. 

Steps 

Remove  the  mean  of  SO)  from  each  station  term  so  that  £  SO)  -  0  . 

I 

Steps 

For  each  source-station  pair,  (l,  D.  compute  M(i.£  ■  E(i)  ♦  SO)  ♦  F(ki). 

Steps 

For  maaknum-Htelihood  estknetors.  compute  otMLE)  via 
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H*oi  -  m  )* 

- ti - 

For  toralivt  toatt-iquarM  Mtimalors.  oomputo  o<ILS)  via 

o*. - C! - 


whara  andtaaia  oolaciioflB  of  tha  obaarvad  station  magniudM  of  aach 

typa.  raapactMfy.  and 

1  -U*  r“ 


♦(U)  • 


ara  tha  probabiity  dansily  function  and  probabity  dWrtwUon  function,  raapactivaly.  of  tha  stan¬ 
dard  normat  random  variabta. 

If  no  oanaorad  lacordbiQ  la  invofvad.  |181  and  |i9i  would  ba  aqualy  appicabia  sinca  n,  •  n^  >  0. 
Slap7 

Rapfaoa  oanaorad  and  mfaslno  obaarvatlona  mm  «dm  via  oorraspondino  oontftionai  aipactationa: 
For  typa-1  pa0«:  E|XlX<l^i«  p  - 

For  lypa-2  pad«:  E|XlX>m  -  u* 

•i”«af 


For  lypa-^  paitia;  E  {  X  I  X  Is  misafno  i  •  P(U)  -  E(f)  ♦  S(0  ♦  F(liJ  . 

Thasa  oondUonaf  axpactallorw  ara  twn  uaad  as  mm  siapa  1  through  3. 

Slaps 

Ftapaai  alapa  to  updaia  E.  8.  F.  and  o  uml  oonvarganoa. 

In  Sia  IrM  sariSon.  only  typa  0  data  ara  uaad  in  siapa  1  through  3.  StarPng  Irom  tia  taoond  loop^ 


paaudo- obaarvatlona"  aa  daacrtiad  In  Map  t.  tn  oStar  words,  tia  ayrrbof  X(VA  In  oMpa  1-3  acfcisly 
rapraaanla  tha  oondSlonal  aspactiaion of  X  ghran Sia  oanaortrig or  norvcarabrfnQ  asaurnptlori  For^rpa- 
0  data.  E|XlX«S(||*i||,and  hanoa  tha  aduafy  obaarvad  mognauds  Is  uUKtad  in  aaoh  SaraSon 
wShout  changa.  For  oSiar  typM  of  data,  howavar.  vw  "aspaolad*  obaanraUon  ««  ba  varying  as  tha 
tar  at  ions  procaada.  sinoa  tha  opdmal  aattnala  of  o  and  al  ofhar  paramaisrs  aHB  clwnga  at  aaoh  alsp. 


Onoa  tha  "E  slop"  (vtr. .  sMpa  5  and  7)  Is  asacuiad.  tia  H  slap"  (vCr. .  aiapo  i  tfvough  4)  in  aach 
laraton  loop  can  ba  rapfaoad  wfti  standard  naaib  in»  anion  tachnfqMa  such  aa  Stputor  Valia  Otoo/^^^ 
poattort  isvcq  or  OmmJm  9Mmon  madiod.  To  do  so.  typa-3  patis  shoiAf  ba  asdudsd  bom  slap  7. 
Numaricat  ataorthma  Bta  SVD  and  Qausaian  afknfnsdon  ara  caSad  dbacr  iwtriatw  iiowavar.  draol 
madiods  can  ba  irnpracitcai  f  tfia  dsaign  mabti  is  larga  and  sparaa.  in  our  casa.  tia  Inaar  syalam 
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invoivM  3,306  unknown  parameters  and  35,21 1  station  nvignitudes.  For  these  types  ot  problems,  Kara- 
live  methods  are  superior  to  Gaussian  elimi/usion  ar>d  matrix  factorization.  The  largest  area  for  the 
application  of  iterative  niethods  is  that  of  the  lirtear  systems  arising  in  the  numerical  solution  of  partial 
differentiai  equations.  Systeritt  of  orders  10,000  to  100,000  are  not  urKiSual  in  aerospace  sciences, 
although  the  majority  of  the  coefficients  of  the  systems  are  typically  zeros. 


IS 


IV.  GLM  INVERSION  RESULTS 
IV.1  GLII  Evmm  MagnitudM 

Table  5  gives  the  Istkig  of  event  tnagnibKles  based  on  six  time-don^  measurements:  . 

n%(Pi),  r\{P(,),  n\{Pn).  and  Ms(LR)  of  clustered  explosions.  This  covers  the  major  tests 

including  NTS,  DM,  SR,  NN2,  and  FS.  fow  hundred  sixty  two  events  were  used  in  the  GLM  inversion 
For  the  Lf  phase,  since  the  path  attenuation  has  been  accounted  for  in  computing  the  raw  station  mag¬ 
nitude,  Equation  [16]  was  used  to  determine  the  event  r%(Lg)  and  the  station  terms.  AH  the  remaining 
magnitudes  are  Wg  ■  That  is,  the  resulting  event  m^nitudes  are  the  maximum-Ukeiihood  netwoik  aver¬ 
age  of  the  station  magnitudes  that  have  been  corrected  lor  the  path  and  receiver  effects.  The  event 
magnitudes  based  on  these  6  phases  have  been  used  in  the  magnitude-magnlude  regressions,  and  the 
regression  results  are  summarized  in  a  later  section. 

As  mentioned  in  Section  1 ,  there  are  actually  many  other  events  lor  which  the  station  rr^  (or  other 
magnitudes)  were  trteasured  urKfer  Task  1,  and  yet  they  were  not  used  in  this  particular  GLM  inversloa 
The  model  descrt)ed  in  Equation  [17]  works  best  when  there  exists  a  suHe  of  events  clustered  in  the 
same  geological  /  geophysical  regioa  That  is,  the  redundancy  of  events  in  each  source  region  Is  the  key 
to  oon^rain  the  unknowns,  particularly  the  path  terms.  PNEs  and  explosions  in  Amchitka,  South 
Ncvaya  Zemlya,  Colorado,  and  Mississippi  were  not  fed  into  the  same  GLM  inversioa  The  magnitudes 
of  these  107  events  wrere  computed  with  anothur  algorthm,  and  the  results  are  listed  in  Appendb  C. 
These  event  magnHudes  are  called  the  EMILS  magnitudes,  /.e. ,  the  sir^le-event  network  averages,  witti 
the  station  terms  Inferred  from  the  GLM  inversion  removed.  They  are  furnished  by  Baumstark  and 
Wagner  (1994),  and  Ksted  in  this  report  merely  lor  reference. 

To  Hustrate  the  robustness  of  the  proposed  determination  scheme.  Figures  1  through  6  show 
the  raw  and  corrected  rr^  of  6  typical  explosions  from  Eastern  Kazakhstan  and  Novaya  Zemlya.  Fi(^ 
1  shows  the  scatter  plot  of  3  difterers  types  of  station  rr^  s  for  Soviet  Joint  Verification  Experiment  [JVE] 
explosion  SH680914.  The  26  good  recordings  artd  1  clipped  sigrtai  are  shown  wHh  fated  circles  arxf 
upward  arrow,  respectively.  The  raw  station  m^  s  (top)  have  a  standard  deviation  of  0.30  m.u.  Applying 
the  “primary*  station  ooaections  (l.e. ,  the  **Rcv“  cotunm  in  Table  6  or  the  *S*  term  in  Equation  [17]}  and 
the  '‘secondary*  corrections  {1.9. ,  the  “Path*  column  in  Table  6  or  “P'  term  in  Equation  [17])  reduces  the 
scatter  down  to  0.14  mot.  The  dashed  ines  of  1  o  range  around  the  network-averaged  dearly  Htos- 
trate  the  remarkablo  reduction  of  fluctuation  across  the  recording  stations.  The  mean  event  n\  Itself  Is 
not  significamy  changed,  however. 

Among  these  six  explosions,  the  Novaya  Zemlya  event  NZ^1014  (Figure  5)  has  the  smalest 
scatter  In  the  resulting  mo  values.  The  dramatic  reduction  of  variation  from  m,  to  itiq  shows  a  factor  of 
nearly  3.  Nde  that  the  path  correction  proposed  in  this  study  not  only  reduces  the  scatter  for  large 
everks,  I  also  Improves  the  data  consistency  tor  the  smal  events,  as  Mudrated  by  Figures  4  and  6. 
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Table  S.  QLM  Event  Magnitudes 


Event 


US660224 


US660414 


US670523 


US670526 


US680322 


US68061S 


US681208 


US690507 


US691008 


US730606 


US750514 


US750603 


US7S0619 


US7S1120 


US760214 


US760309 


US760317 


US780411 


US780411 


US780831 


US781216 


US790611 


US790926 


US800426 


US800612 


US800725 


US801217 


US810606 


US820212 


US820212 


US820425 


US820624 


US830326 


US830901 


US840725 


US841209 


4.91  ±0. 


S.OGtO. 


5.33±0.04 


6.1410.04 


5.6310.03 


.04 


5.7010.03 


5.9210.04 


5.3410.05 


5.9210. 


.07 


5.7910.03 


5.5710.03 


6.2110.03 


5.9910.03 


5.8110.03 


6.0510.03 


5.92t0. 


.04 


5.9410.04 


6.0510. 


5.4610. 


5.5210. 


5.6710. 


71 


5.5610.04 


5.5710.05 


5.5410.04 


5.5810.04 


5.5510.04 


5.1310.07 


5.5110. 


5.4710.07 


5.3610.06 


5.5710.04 


5.2610. 


5.3010.04 


5.3010. 


5.4510. 


4.4810.06 


4.7710.04 


.04 


5.7310.04 


5.4210. 


4.9210.04 


5.2610.03 


5.6110.04 


4.9410. 


5.6510. 


4.2810.07 


5.4910. 


5.3010. 


5.7710.03 


5.5310.03 


5.7710. 


5.7510.04 


5.7210. 


5.6910.04 


5.7410.04 


5.0510. 


5.2310. 


5.3910.06 


5.3410.06 


5.3210.04 


5.3010.05 


5.2510.04 


5.2610.04 


5.3110.04 


4.7210.07 


5.2210. 


5.0410. 


5.2210.07 


5.1010. 


5.3210.04 


4.3310. 


4.5010. 


4.7110.04 


5.4710.04 


5.1210. 


4.7410.04 


4.9610.03 


5.2510.04 


4.4810. 


5.3910. 


3.9210. 


5.2410 


5.0410.04 


5.6710.03 


5.5210.03 


4.7810.09 


5.4710.03 


5.4410.04 


5.4810.04 


5.4010.04 


5.5310.04 


4.8510. 


4.9610.03 


5.1810. 


4.9810. 


5.1110.04 


4.9910. 


4.9410.04 


4.9510.04 


5.0310.04 


4.2810.07 


4.8410. 


4.7510.06 


5.0410.07 


4.8710.07 


5.0510.04 


5.3910.25 


5.2010.11 


5.6210.14 


5.9710.25 


5.4410.25 


6.0710.14 


5.9610.14 


5.9210.12 


6.0010.12 


5.7210.12 


5.8310.12 


5.9110.12 


5.9010.12 


5.3810.10 


5.5710.14 


5.5710.14 


5.6110.12 


5.6010.12 


5.5310.18 


5.5510.11 


5.5610.11 


5.6310.10 


4.9910. 


5.42i0.10 


5.3110. 


5.4810. 


5.2310. 


5.3710.09 


5.2310.10 


5.5110.10 


5.6110.08 


5.7210.12 


5.7410.21 


5.7810. 


6.0010. 


5.1610.15 


6.1110.09 


6.2310.07 


6.0110.07 


6.2110.11 


6.0610.07 


6.3310.15 


6.1310.15 


6.0610.11 


5.6510.10 


5.6310.12 


5.6310.12 


5.9210.15 


5.6910.11 


5.9810.12 


5.7310.11 


5.2210.10 


5.6810. 


5.6210. 


5.7910.10 


5.4610.1  C 


5.7410.07 


5.1010.04 


.04 


5.1110. 


4.8410.04 


4.7910. 


5.3010. 


5.4010. 


5.4110.12 


5.6210.10 


5.6510.12 


5.4610.15 


4.8710.17 


5.2810. 


4.9510.10 


4.6710. 


a 


4.7510.07 


4.5610. 


5.2410.17 


5.2010. 


5.1510. 


4.7110.07 


5.3310.07 


5.1410.10 


5.3610.07 


5.32t0. 


4.9910. 


4.2010. 


4.5710. 


4.4710. 


4.4610. 


4.6810.03 


4.6010.04 


4.6410.04 


4.6210. 


4.6410.04 


21 

a 

a 


4.3910. 


4.6710.05 


4.5310.04 


4.5210.04 


4.5210.04 


06  j 

4.9910. 

06 

4.5810.06 

5.0010.10 

5.0710.12 

4.1110. 

4.2710.04 


4.1710.04 


4.3010. 


Event 

Site 

US84121S 

PMA 

US850502 

PMA 

US850612 

PMA 

US850725 

PMA 

US851228 

PMA 

US860422 

PMA 

US860625 

PMA 

US860717 

PMA 

US860930 

PMA 

US861016 

PMA 

US861213 

PMA 

US870418 

PMA 

US870430 

PMA 

US870924 

PMA 

US880215 

PMA 

US880817 

PMA 

US890622 

PMA 

US891208 

PMA 

US631016 

RNA 

US670626 

RNA 

US670831 

RNA 

US680229 

RNA 

US680924 

RNA 

US680115 

RNA 

US690212 

RNA 

US691205 

RNA 

US700505 

RNA 

US700526 

RNA 

US710629 

RNA 

US720502 

RNA 

US720720 

RNA 

US730605 

RNA 

US731012 

RNA 

US740619 

RNA 

US750405 

RNA 

US751024 

RNA 

US760512 

RNA 

US780913 

RNA 

US820923 

RNA 

US840215 

RNA 

5.38±0. 


5.39±0.04 


52410. 


5.2510. 


5.2910. 


5.4410.04 


5.5710. 


5.4210.04 


5.5210.04 


5.5510. 


5.4910.04 


5.4310.04 


5.5710.04 


5.3110.05 


5.4910.04 


5.1010.12 


5.3910. 


5.3510. 


5.0110. 


4.9910.10 


5.1410.07 


5.3310. 


5.3610.04 


5.0510.17 


5.1910.10 


5.1610.07 


4.8110. 


.10 


.12 


5.0610. 


5.0010.07 


4.8010.12 


4.9210.10 


.12 


4.9510.12 


4.8810.12 


.17 


5.3010.07 


5.2310.07 


Table  5.  GLM  Event  Magnitudes  (Cont) 


Pt 


5.1010.04 


5.3710.04 


5.1910.04 


5.0310. 


4.8710. 


4.9710.05 


5.1310.04 


5.3010.03 


5.1610.04 


5.1410.04 


5.1010.05 


5.1410.04 


5.1410.04 


5.2710.04 


5.0110.06 


5.1810.04 


5.0610.12 


5.3710.08 


4.7910.04 


5.0310. 


4.9910.04 


4.8210.05 


4.8010.06 


4.7510.05 


4.7810.04 


4.9910.04 


4.9010.04 


4.7810.04 


4.6110. 


4.8710.04 


4.8410.04 


5.0010.04 


5.0310.06 


4.6510.08 


4.7410.10 


4.9110.04 


4.9010.12 


5.1110.10 


5.0810.06 


5.2910.11 


5.5210.12 


5.4410.11 


5.1410.18 


5.2710.11 


5.3410.12 


5.3510.10 


5.4510.14 


5.3910.12 


5.3910.11 


5.4110.12 


5.3810.09 


5.3410.12 


5.4110.12 


5.2810.12 


5.4010.12 


5.2310.25 


5.6010.21 


6.0710.15 


5.6710.21 


5.4210.21 


5.4510.21 


5.6610.15 


5.6210.15 


5.8410.10 


5.6710.21 


5.8510.21 


5.7210.21 


5.7410.15 


5.6910.15 


5.7010.21 


Ms 


4.2310.04 


4.3010.05 


4.5110.04 


4.2810.04 


4.2810.04 


4.3310.04 


4.3010.03 


4.4510.04 


4.5010.04 


4.5310.04 


4.6010.05 


4.3710.04 


4.5110.04 


4.4510.03 


4.2710.04 


4.3910.04 


4.1210.10 


4.2210.09 


4.7710.10 


5.1510.04 


4.6710.10 


4.7910.08 


4.4410.08 


4.5810.10 


4.5510.12 


4.8010. 


4.6810.07 


4.5010.10 


4.7110.12 


4.6910.12 


4.5310.12 


4.9410.07 


4.8710.07 


4.6910.08 


4.8410.04 


4.4810.10 


4.4610.08 


4.2310.10 


4.5010.10 


4.2610.12 


4.7710.10 


4.5110.07 


4.1910.10 


4.5510.10 


4.5510.07 


4.9410.10 


4.4510.10 


4.3810. 


4.6910.09 


4.6010.11 


4.9410.12 


4.3110.12 


4.4810.10 


4.6310.14 


4.2310.12 


4.5010.11 


4.4910.12 


4.5810.12 


4.3010.12 


4.1510.12 


4.5710.12 


4.5710.12 


4.4310.12 


4.5710.14 


4.7010.10 


4.8210.11 


5.5110.09 


5.1610.11 


4.9810. 


5.1810.10 


5.2910.10 


5.5910.09 


4.9510.11 


5.0510.10 


4.9210.15 


4.5510.12 


4.8510.12 


5.1110.11 


5.1210.12 


5.2610.10 


4.8210.11 


4.9710.12 


4.9610.15 


5.0810.15 


4.9810.15 


4.8310.15 


5.2310.15 


5.1310.21 


3.6110.12 


3.7610.06 


3.6310.06 
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Event 


US620512 


US620627 


US620713 


US621005 


US630522 


US630913 


US631122 


US640116 


US640220 


US640424 


US640716 


US641009 


US641205 


US650218 


US650303 


US650326 


US650S14 


US650521 


US650723 


US650910 


US651203 


US651216 


US660118 


US660513 


US660519 


US660527 


US660602 


US660603 


US660615 


US660625 


US670119 


US670120 


US670223 


US670510 


US670520 


US670727 


US670907 


US670927 


US671018 


Table  5.  GLM  Event  Magnitudes 


Pb 


(Cont) 


4.9G±0. 


4.73±0. 


6.74±0. 


5.13±0.i 


4.85±0. 


5.81±0.i 


5.02±0. 


5.25±0.i 


5.1 1±0. 


4.99±0. 


4.5Q±0. 


4.76±0.i 


4.81±0.i 


4.47±0. 


5.31±0.i 


5.15±0.i 


4.73±0. 


4.69±0.i 


5.43±0.i 


S.05±0.i 


5.63±0.i 


5.22±0.i 


5.13±0.i 


5.56±0J 


5.76±0.i 


4.96±0.i 


5.55±0.i 


5.66±0.< 


5.07±0.i 


4.6110.1 


5.3110. 


5.1210.1 


5.7010.1 


4.9810.1 


5.9110.1 


4.7810. 


4.9910. 


5.8010. 


5.6810.1 


.8910.12 


.5710.10 


.4210.17 


i.1610.17 


.8110.07 


.7610.10 


1.5710.03 


.80±0.07 


.6410. 


.5510. 


.1010.12 


.5110.07 


.3110.06 


.2710.12 


.9610.04 


-.9410. 


.9710.17 


.0110.08 


.0410.04 


.7810.06 


i.3110.04 


i.8510.04 


1.8110. 


i.0810.10 


1.4610.03 


1.7110.06 


1.3010.03 


1.2410.03 


.7610. 


-.0710. 


1.9810. 


1.7010.07 


1.30±0.03 


.3710.07 


1.6010.03 


1.1910.10 


1.7810.06 


1.4810. 


1.4310. 


4.7210.12 


4.3010.10 


3.9010.17 


5.5910.17 


4.7510.07 


4.4910.10 


5.3810.03 


4.6810.07 


4.3210. 


4.1810.06 


4.3310.12 


4.5210.10 


4.2010.07 


4.4210.12 


4.7210.04 


4.8610.06 


4.1110.10 


4.8010.04 


4.4910.06 


5.0610.04 


4.6710.05 


4.5410.06 


4.9510.08 


5.1710.03 


4.1610.07 


4.9710.03 


4.9910.03 


4.4210.06 


3.8310.10 


4.6610.06 


4.5010.07 


4.1210. 


5.3710.03 


4.4710.10 


4.4010. 


5.2110.04 


5.1710.03 


4.9210.25 


4.8010.18 


4.3910.18 


3.9710.25 


4.8810.14 


4.4610.11 


5.7410.18 


4.5410.14 


4.9210.14 


4.8310.10 


4.9110.12 


4.2510.10 


4.5710.10 


4.6710.11 


4.1410.10 


5.0310.11 


5.1810.11 


4.4410.10 


4.4410.09 


5.2910.11 


4.9010.14 


5.5310.14 


5.0010.09 


4.9510.10 


5.42t0.14 


5.4010.12 


4.9010.14 


5.6110.14 


4.9510.14 


4.3410.11 


5.1810.11 


4.9810.14 


5.6310.18 


4.6010.11 


5.9610.25 


4.3610.12 


4.8710. 


5.7710.14 


5.6410.11 


5.4210.21 


4.8610.21 


4.7810.21 


5.1710.21 


6.1010.10 


4.7910.21 


5.3210.21 


5.6410.21 


4.9610.15 


5.1310.21 


5.1410.21 


4.9510.15 


5.6710.10 


4.7110.15 


5.1410.15 


5.7010.21 


4.2610.09 


5.6410. 


5.8310. 


6.0610. 


5.4410.21 


5.7910.10 


5.9510.07 


5.3210.21 


5.0410.21 


6.1010.12 


5.5710.15 


6.1710.12 


4.7910.15 


4.3910.06 


3.7910. 


4.8210.17 


4.3510.17 


4.5410.17 


3.5210.07 


4.6310.05 


4.4610.12 


3.6510.06 


.11 


6.0010.10 


3.8710.07 


4.7110.10 


Event 


US671108 


US680119 


US680221 


US680325 


US680410 


US680418 


US680517 


US680906 


US680917 


US690130 


US690321 


US690S27 


US690716 


US690716 


US691029 


US691121 


US691217 


US691217 


US700204 


US700205 


US700225 


US700323 


US700521 


US701014 


US701105 


US701216 


US701217 


US710624 


US710708 


US710818 


US720519 


US720921 


US721221 


US730308 


US730426 


US730628 


US740227 


US740523 


US740710 


Table  5.  GLM  Event  Magnitudes 


(Cont) 


4.84±0. 


4.68±0.17 


5.66±0.07 


4.45±0.12 


4.78±0.12 


5.08±0. 


4.89±0.05 


5.55±0.03 


4.81±0.07 


4.91±0.10 


4.9410.08 


5.33±0.04 


4.9110. 


4.4010.12 


5.4510.04 


5.6710.03 


4.7910. 


5.4510.04 


5.6110.12 


5.6710.04 


4.7210. 


5.0510.04 


5.4810.07 


4.9910. 


5.5710.04 


4.8310.10 


5.1810.04 


5.8710.03 


5.1310.06 


5.5110.03 


5.2610.04 


4.6510.10 


5.6810.03 


5.0510. 


5.3810.04 


5.5410.03 


4.9310.05 


5.6510.04 


4.8410. 


5.7410. 


4.4810.17 


5.5210.07 


4.4010.17 


5.3110.07 


3.9610.12 


4.7410. 


4.3210.05 


5.2510.03 


4.4910.07 


4.2510.10 


4.2910.12 


4.9410.04 


4.6310. 


5.1010.04 


5.2710.03 


4.4610.10 


5.1510.04 


4.8110.12 


5.3310.04 


4.3310. 


4.6710.04 


5.0510. 


4.4710.05 


4.3210.05 


5.0210.04 


4.1810.08 


4.2610.12 


4.6210.04 


4.2310.06 


4.8910.04 


4.9210.03 


4.4210.10 


4.9410.04 


4.6110.12 


4.0210.10 


4.5110.05 


4.6110.08 


4.2110.07 


5.0210.04 


4.4210.10 


4.7510.04 


5.5810.03 


4.4710.06 


5.1610.03 


4.8210.04 


5.2610.03 


4.6210.06 


5.0410.04 


5.1310. 


4.4610. 


5.3810.04 


4.4010.05 


5.1910.03 


4.3610.07 


4.9210.03 


4.4910.04 


4.9910.03 


4.4510.07 


4.7010.05 


5.0010.03 


4.2110.06 


4.7410. 


4.5310.18 


5.6010.14 


4.2110.12 


4.3710.10 


4.8910.18 


4.7510. 


5.5410.18 


4.5610. 


4.6910.10 


4.6510.11 


5.0410.14 


4.8110.11 


4.2910.14 


5.4110.11 


4.7710.12 


5.2210.18 


4.4610.14 


5.4710.14 


4.4010. 


4.8310. 


5.3310.14 


4.6510.10 


5.4510.18 


4.6610.11 


5.0410.14 


5.7810.25 


4.7410. 


5.4510.14 


4.9710.11 


4.4410.25 


5.5510.18 


4.7110.10 


5.0810.12 


5.3010.12 


4.7610. 


4.2810.07 

4.1110.07 

5.2910.03 

5.0310.03 

4.3810. 


6.1110.10 


4.7310.15 


5.3310.21 


5.0310.15 


6.0410.11 


5.4110.15 


4.9110.15 


5.5610.10 


5.3310.21 


5.7110.21 


5.8510.11 


5.3010.21 


4.9810.15 


5.4910.21 


5.7810. 


4.9110.15 


5.6910.15 


5.0810.21 


5.9510.12 


5.2110.21 


5.8310.07 


5.8310. 


5.0610.15 


6.0710.12 


5.6210.07 


5.8610. 


5.1810.11 


5.1210.21 


5.6610.12 


3.5710.07 


4.8710.12 


3.8610.17 


4.2110.17 


4.2710. 


4.6110.12 


4.9110.17 


3.7810. 


4.4810. 


.17 


4.2510.08 


4.4810.10 


4.5310.05 
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Table 


Event 

Site 

US740630 

YFT 

US740926 

YFT 

US750228 

YFT 

US750307 

YFT 

US750424 

YFT 

US750430 

YFT 

US750603 

YFT 

US751220 

YFT 

US760204 

YFT 

US760204 

YFT 

US760317 

YFT 

US760727 

YFT 

US760826 

YFT 

US761228 

YFT 

US770405 

YFT 

US770427 

YFT 

US770525 

YFT 

US770804 

YFT 

US770819 

YFT 

US770927 

YFT 

US771109 

YFT 

US771214 

YFT 

US780223 

YFT 

US780323 

YFT 

US780712 

YFT 

US780927 

YFT 

US780927 

YFT 

US781118 

YFT 

US790208 

YFT 

US790215 

YFT 

US790628 

YFT 

US790803 

YFT 

US790808 

YFT 

US790829 

YFT 

US790906 

YFT 

US800403 

YFT 

US800416 

YFT 

US810115 

YFT 

US811001 

YFT 

US811111 

YFT 

5.71±0.( 


5.54±0.l 


5.76±0. 


5.61±0. 


4.57±0. 


5.12±0. 


5.70±0. 


5.79±0.( 


5.e9±0.( 


5.73±0.( 


5.91±0.l 


5.21±0.( 


5.26±0. 


5.59±0.l 


5.71±0.( 


S.42±0.l 


5.43±0.( 


5.21±0.( 


S.66±0.l 


4.97±0.( 


5.81±0.l 


5.81±0.i 


5.6710.I 


5.66±0.i 


S.68±0J 


5.23±0.i 


5.92±0.i 


5.3610.I 


5.66±0J 


5.06±0.i 


5.18±0. 


4.74±0.( 


4.76±0.i 


5.06±0.l 


5.96±0J 


4.9910J 


5.45±0. 


5.70±0.l 


5.04±0.l 


4.96±0.i 


5.  GLM  Event  Magnitudes 


Pt 


48±0.05  5.33±0. 


(Cent) 


i.10±0.04 


>.43±0.03 


).27±0. 


I.96±0.10 


1.8310.05 


>.4210.03 


>.4810.03 


>.3610.03 


>.4910.03 


>.7110.04 


>0210.08 


1.8910.05 


>.2510.04 


>.4610.04 


>.1710.04 


>.0510.05 


t.8510.04 


>.4610.04 


t.54l0. 


i.5210.03 


>.5610.04 


>.3810. 


>.3110. 


>.3110. 


t.85l0.04 


>6210.04 


>.0610. 


>.3810.04 


t.6910.07 


1.8910.05 


1.4010.06 


1.4110.07 


1.6410.08 


>.6310. 


1.5910.07 


>.2610.08 


>.2910.04 


1.5210.07 


1.4710.07 


4.8410.04 


5.2310. 


4.9310. 


3.6910.10 


4.7010. 


5.1910.03 


5.1510.04 


5.2910.03 


5.4410.04 


4.7410. 


4.4510.06 


5.0510.04 


5.0810.04 


4.9210.04 


4.8110.05 


5.0710.04 


4.3410. 


5.2810.04 


5.2910.04 


5.0510.04 


5.0510.03 


5.0010.03 


4.5510.04 


4.8210. 


4.3810.07 


4.4610. 


4.2110. 


4.4010. 


4.5710. 


5.5610.25 


5.4310.14 


5.6310.12 


5.5510.14 


4.5110.10 


4.8510.12 


5.4410.14 


5.8310.11 


5.3010.12 


5.3210.18 


5.5010.25 


4.9110.18 


5.1510.12 


5.2710.14 


5.6010.14 


5.2410.10 


5.2910.11 


5.0610.08 


5.3710.25 


4.9110.11 


5.6410.18 


5.6110.14 


5.6410.11 


5.2710.12 


5.0110.11 


5.6210.25 


5.1910.11 


5.4510.12 


4.8710.10 


5.0510.10 


4.7010.11 


4.5910. 


4.8310.12 


5.8610.12 


5.6410.10 


5.8410. 


4.8310.15 


5.5810.10 


5.7810.15 


5.9210.11 


5.8610.09 


Ms 


4.5610. 


5.6710.15 


5.8410.10 


5.7810.07 


5.5510.10 


5.5910.10 


6.2310.21 


5.4610.15 


5.5810.21 


5.8810.10 


5.4710.21 


5.5410.15 


6.1510.10 


5.3310.15 


5.6610.11 


5.2310.12 


5.2910.11 


5.5110.15 


4.5810.07 


4.6810.10 


4.7410. 


4.0810.07 


4.3810. 


4.3210.10 


4.3410.07 


4.0610. 


4.5910.12 


4.8710.06 


4.8210.07 


4.3610. 


4.6610.07 


4.3710.04 


4.7210. 


4.3610. 


4.5310. 


03 

5.3710. 

[13 

5.6510.11 

6.0910. 

08 

4.7710.04 

4.5010. 


4.9710. 


5.0010.04 


4.3810.07 


4.3110.07 


4.6310.14 


5.2110.11 


5.3010.11 


4.9310. 


4.7110. 


5.3010.21 


5.8910.11 


5.8610.07 


5.4110. 


5.3110.10 


4.4710. 


4.7310.04 


4.1110.17 


3.5710.17 


Table  5.  GLM  Event  Magnitudes  (Cent) 


US811112 

YFT 

5.4310.04 

5.1110.04 

4.7810.04 

US811203 

YFT 

4.9410.17 

^1 

! 

US820128 

YFT 

5.9110.04 

5.5810.04 

5.3310.04 

US820507 

YFT 

5.6410.04 

5.3210.04 

5.0710.04 

US820729 

YFT 

4.5410.10 

3.8710.12 

4.0310.12 

US82080S 

YFT 

5.7710.04 

5.5110.04 

5.2610.04 

US821210 

YFT 

4.7210.07 

4.2710.07 

3.8010.08 

US830414 

YFT 

5.6510.04 

5.2610.04 

5.0810.04 

US831216 

YFT 

5.4910.07 

4.8410.10 

4.0410.17 

US840131 

YFT 

4.4410.12 

4.0810.12 

4.3710.17 

US840301 

YFT 

5.8310.04 

5.4610.04 

5.2510.04 

US840501 

YFT 

5.4810.04 

5.1810.04 

4.8810.04 

US840531 

YFT 

5.7510.05 

5.4710.05 

5.1610.05 

US840620 

YFT 

4.7410.07 

4.2810.07 

4.2310.07 

US840913 

YFT 

5.1010.07 

4.7410.07 

4.5110.07 

US850315 

YFT 

4.8410.06 

4.5110.06 

4.3110.07 

US850323 

YFT 

5.3410.04 

5.0410.04 

4.7510.04 

US850402 

YFT 

5.7410.04 

5.4210.04 

5.2310.04 

US851205 

YFT 

5.5510.04 

5.2610.04 

4.9510.04 

US860322 

YFT 

5.2910.05 

4.7710.05 

4.5110.06 

US860605 

YFT 

5.3310.05 

4.9810.05 

4.7610.05 

US861114 

YFT 

5.6910.04 

5.4010.04 

5.0010.04 

US870716 

YFT 

4.9410.07 

4.5610.07 

4.3110.07 

US870813 

YFT 

5.8610.04 

5.5810.04 

5.2710.04 

US871023 

YFT 

5.4610.06 

4.8110.06 

4.5010.06 

US881013 

YFT 

5.8810.06 

5.3610.06 

4.8910.06 

US890210 

YFT 

5.2610.12 

4.9510.12 

4.7810.12 

US890309 

YFT 

5.1910.17 

4.8210.17 

4.7110.17 

5.20±0.09 


4.58±0. 


5.54±0.09 


5.45±0.09 


4.44±0.07 


5.46±0.09 


4.45±0.08 


4.67±0.11 


4.24±0.08 


5.2210.11 


5.4710.09 


4.3810.11 


4.6610.11 


4.4910.10 


4.9910.12 


5.3810.12 


5.2110.11 


5.3810.11 


4.6910.18 


5.5210.14 


5.6710. 


5.1710.08 


5.9210.07 


5.8810.07 


4.9010.09 


5.7610.08 


4.7610.11 


5.8410.07 


5.3610.15 


4.9310.11 


5.1810.12 


5.8010.11 


6.1110.15 


5.5910.15 


5.9110.21 


5.5010.21 


6.1610.21 


5.4010.15 


SH680619 

BNE 

5.3910.03 

5.1310.03 

4.8910.03 

SH720210 

BNE 

5.4410.03 

5.1910.03 

4.9710.03 

SH721210 

BNE 

6.1010.03 

5.8910.03 

5.6910.03 

SH731214 

BNE 

5.8810.04 

5.7010.04 

5.4410.04 

SH741016 

BNE 

5.5310.03 

5.3010.03 

5.0610.03 

SH741227 

BNE 

5.6610.03 

5.4310.03 

5.2110.03 

SH751225 

BNE 

5.8010.03 

5.5010.03 

5.2410.03 

SH760609 

BNE 

5.2510.03 

5.0310.03 

4.7510.03 

SH760828 

BNE 

5.7610.03 

5.5410.03 

5.2310.03 

SH761123 

BNE 

5.8510.03 

5.6710.03 

5.40+0.03 

SH770629 

BNE 

5.3110.03 

4.9810.03  ! 

4.7710.03 

SH770905 

BNE 

5.8710.03 

5.6210.03 

5.3210.03 

5.4710.25 


Ms 


4.3510.04 


4.7710.04 


4.3710.04 


3.3210.07 


4.8410.04 


3.1010.07 


4.3510.04 


4.0410.04 


4.5910.04 


4.6010.04 


4.6310.04 


3.4510.05 


3.9110. 


3.3910. 


4.1910.04 


4.8110.04 


4.3910.04 


3.9710.07 


4.4110.04 


4.5910.04 


3.8310.05 


4.6210.04 


3.9010. 


4.7710. 


4.3510.12 


5.9910.10 


5.7810.11 


5.4210.21 


5.5710.12 


4.9710.25 


4.8610.25 


5.5710.25 


4.9410.25 


4.2510.25 


6.2010.18 


5.8310.12 


4.7210.04 


4.5910. 


3.6110. 


3.9410. 


3.2910.12 


3.9610. 


3.9810. 


4.0810. 


Event 

Site 

SH771029 

BNE 

SH780829 

BNE 

SH781104 

BNE 

SH790201 

BNE 

SH790707 

BNE 

SH791028 

BNE 

SH800612 

BNE 

SH801012 

BNE 

SH801227 

BNE 

SH810329 

BNE 

SH810527 

BNE 

SH821226 

BNE 

SH831120 

BNE 

SH840307 

BNE 

SH840526 

BNE 

SH841202 

BNE 

SH881112 

BNE 

SH890902 

BNE 

SH721102 

BSW 

SH760421 

BSW 

SH761207 

BSW 

SH770529 

BSW 

SH780611 

BSW 

SH780705 

BSW 

SH780915 

BSW 

SH781129 

BSW 

SH790623 

BSW 

SH790804 

BSW 

SH791202 

BSW 

SH791223 

BSW 

SH800629 

BSW 

SH800914 

BSW 

SH810422 

BSW 

SH811018 

BSW 

SH811129 

BSW 

SH811227 

BSW 

SH820704 

BSW 

SH820831 

BSW 

SH821205 

BSW 

SH831006 

BSW 

Table  5.  GLM  Event  Magnitudes 


Pb 


(Cent) 


5.80±0.C 


5.94±0.( 


5.64±0. 


5.42±0.( 


5.9e±0.C 


6.09±0.C 


5.66±0. 


5.96±0. 


S.91±0.( 


5.59±0.( 


S.40±0.( 


5.74±0.( 


5.52±0.( 


5.75±0. 


6.12±0. 


5.94±0.( 


S.36±0.( 


5.01±0. 


6.30±0. 


5.13±0.1 


5.85±0. 


5.83±0.( 


5.93±0. 


5.84±0. 


5.96±0.( 


6.04±0.( 


6.13±0. 


6.21±0. 


6.04±0. 


6.21±0. 


5.83±0. 


6.26±0.( 


6.00±0. 


6.05±0. 


5.63±0. 


6.22±0. 


6.25±0. 


5.27±0.( 


6.17±0. 


6.03±0. 


5.48dt0.04 


5.84±0.04 


5.47±0.03 


S.13±0.03 


5.7810.04 


5.8610.03 


5.4710.03 


5.7410.03 


5.7410.03 


5.3710.03 


5.1610.04 


5.5010.03 


5.2710.03 


5.5010.03 


5.9610.03 


5.7410.03 


5.1610.04 


4.8110.06 


6.1010.03 


4.9710.12 


5.6310.03 


5.5810.03 


5.7110.03 


5.6610.03 


5.7910.04 


5.8210.04 


5.9810.03 


6.0210.03 


5.8310.03 


6.0010.03 


5.6410.03 


6.0910.03 


5.7910.03 


5.8510.03 


5.4510.03 


6.0410.03 


6.0310.03 


5.0310.04 


5.9610.03 


5.8010.03 


5.4010.05 


5.6210.04 


5.2310.03 


4.8710.03 


5.5010.03 


5.5910.03 


5.1610.03 


5.4910.03 


5.4710.03 


5.0610.03 


4.8910.04 


5.2110.03 


4.9810.03 


5.2210.03 


5.6610.03 


5.4310.03 


4.9610.05 


4.6210.06 


5.8010.03 


4.7210.12 


5.3510.03 


5.2710.03 


5.3810.03 


5.3010.03 


5.4910.04 


5.5510.04 


5.7510.03 


5.5610.03 


5.7110.03 


5.3710.04 


5.8110.03 


5.5110.03 


5.5810.03 


5.1410.03 


5.7510.03 


5.7310.03 


4.7110.04 


5.6910.03 


5.5110.03 


Ms 


4.1510. 


5.0110.25 


5.9710.18 


5.7510.25 


5.6310.14 


5.5510.14 


5.2510.25 


5.6910.11 


5.7710.12 


5.5310.12 


5.5810.12 


5.7110.12 


5.6010.21 


5.3510.09 


4.6410.25 


6.3110.25 


5.9910.18 


5.8110.14 


5.8810.21 


5.8310.12 


5.8210.15 


.04 


3.3810. 


4.4610. 


4.4610. 


3.7010.04 


4.4710.04 


4.0910.04 


3.8510. 


3.5610. 


3.8210.04 


3.7010. 


4.0410.04 


4.5010. 


4.1210. 


3.2910. 


6.0210.10 


6.0610.25 


5.9010.18 


5.8910.25 


5.9310.21 


5.4710.11 


5.6510.12 


5.7010.15 


5.7910.12 


5.6710.11 


5.8710.15 


6.0010.12 


5.8810.21 


6.0110.25 


7.3110.25 


6.1910.25 


5.2310.25 


6.6310.25 


5.7410.25 


5.5510.25 


5.7810.15 


4.4810. 


2.9310.17 


4.1110. 


3.8210. 


4.4610.04 


3.8910.04 


4.2610.04 


4.2810.04 


4.4210.04 


4.5010. 


4.5210.04 


4.2910. 


3.8310.04 


4.3610. 


4.4010.04 


4.3910.04 


4.2510.04 


4.4410. 


3.8910.04 


4.4310. 


4.4710. 


Table  5.  GLM  Event  Magnitudes  (Cent) 


Event 


SH831026 

BSW 

SH840219 

BSW 

SH840425 

BSW 

SH840714 

BSW 

SH841027 

BSW 

SH841216 

BSW 

SH841228 

BSW 

SH850210 

BSW 

SH850615 

BSW 

SH850630 

BSW 

SH850720 

BSW 

SH870312 

BSW 

SH870403 

BSW 

SH870417 

BSW 

SH870620 

BSW 

SH870802 

BSW 

SH871115 

BSW 

SH871213 

BSW 

SH871227 

BSW 

SH880213 

BSW 

SH880504 

BSW 

SH880914 

BSW 

SH881217 

BSW 

SH890122 

BSW 

SH890212 

BSW 

SH890708 

BSW 

SH650115 

BT2 

SH691130 

BTZ 

SH710630 

BTZ 

SH730723 

BTZ 

SH740531 

BTZ 

SH750427 

BTZ 

SH751029 

BTZ 

SH760704 

BTZ 

SH771130 

BTZ 

SH790818 

BTZ 

SH800425 

BTZ 

SH801214 

BTZ 

SH810913 

BTZ 

SH820425 

BTZ 

6.16±0.03 


5.85±0.03 


5.93±0.03 


6.15±0.03 


6.29±0.03 


6.17±0.03 


6.01±0.03 


5.93±0.03 


6.10±0.03 


5.99±0.03 


5.98±0.03 


5.42±0.03 


S.95±0.03 


5.63±0. 


5.72±0.03 


5.93±0.03 


6.10±0.03 


5.9e±0.03 


S.86±0.03 


5.73±0.03 


5.92±0.03 


5.83±0.03 


5.79±0.03 


5.18±0.03 


S.65±0.03 


5.34±0.03 


5.49±0.03 


S.65±0.03 


5.81±0.03 


S.68±0.03 


5.54±0.03 


5.93±0.25 


5.80±0.18 


6.15±0.25 


5.81±0.18 


6.77±0.25 


5.7610.25 


6.0510.03 


6.1110.03 


5.8810.04 


6.0810.03 


6.1310.04 


6.1810.04 


6.1010.04 


5.8310.04 


6.0710.04 


5.8510.04 


5.6010.03 


6.0910.03 


5.2110.04 


6.3310.04 


5.8810. 


5.6010.03 


5.6410.03 


5.9410.03 


5.8510.03 


6.1310. 


5.4610. 


5.9110. 


6.1210. 


6.0810. 


5.8510.03 


5.9210.03 


5.6810.04 


5.9010.03 


5.9410.04 


5.9510.04 


5.8810.03 


5.9810.03 


5.9610.03 


5.6610.04 


5.8810.04 


5.6310.04 


5.3810.03 


5.8110.04 


5.8210. 


5.0010.04 


6.1110.04 


5.6410.03 


5.3310. 


5.3810.03 


5.7210.03 


5.6210. 


5.9310. 


5.2910. 


5.7510. 


5.9110. 


5.8710. 


5.6010.03 


5.4910.03 


5.4810.03 


4.9010.03 


5.8410.05 


5.5810.03 


5.6610.03 


5.4210.04 


5.8810.15 


5.6610.21 


5.9210.21 


5.9710.21 


5.8710.21 


5.9310.21 


5.8110.15 


6.0010.15 


5.7310.21 


5.7210.04 


5.7510.04 


5.7410.03 


5.6610.03 


5.3910.04 


5.6010.04 


5.4210.04 


5.1010.03 


5.5710.04 


5.5510.03 


4.7610.05 


5.8210.04 


5.3810.03 


4.9910.03 


5.1010.03 


5.4510.03 


5.9010.25 


5.8010.25 


6.1910.25 


6.1510.14 


5.9410.14 


6.1310.14 


6.0010.14 


5.9910.18 


6.1510.18 


6.2210.25 


5.8510.18 


5.7610.14 


5.9710.21 


5.9010.21 


5.6910.03 


5.0310.03 


5.4910.03 


5.6110.03 


5.5510.03 


5.7210.18 


5.9010.18 


5.9010.12 


5.9010.12 


5.3510.12 


5.9510.08 


5.7910.11 


5.4510.15 


5.6010.15 


5.5810. 


5.8810.12 


5.8610.12 


5.3210.15 


5.8410.11 


•a 

a 

•a 


Ms 


4.5010. 


4.3810.03 


4.5410. 


4.5710.03 


4.4810.03 


4.5910.03 


4.3710.03 


4.5110.03 


4.2310.03 


4.3310.03 


4.2210.03 


4.0910.04 


4.7210.09 


4.2610. 


4.2510. 


4.0210.04 


4.6410.03 


4.5010.04 


4.3110.04 


4.4910. 


4.4010.05 


4.5510.03 


4.4010.04 


4.4910.04 


4.5610.04 


4.0810.03 


4.5910.10 


4.3210. 


3.3210.12 


4.3310.05 


3.8610.07 


3.5610.17 


4.2910.04 


3.9510. 


4.1410.04 


3.3410.06 


4.2710.03 


4.5210.04 


4.4310.03 


24 


25 


Event 


DM690723 


DM691001 


DM700129 


DM700628 


DM700724 


DM700906 


DM701217 


DM710322 


DM710425 


DM711129 


DM711230 


DM720310 


DM720328 


DM720607 


DM720816 


DM721210 


DM730216 


DM730710 


DM731026 


DM740130 


DM740516 


DM740710 


DM740913 


DM750220 


DM750311 


DM750606 


DM750807 


DM751213 


DM760115 


DM760723 


DM761230 


DM770329 


DM770425 


DM770730 


DM770617 


DM771029 


DM780326 


DM780422 


DM760728 


DM781015 


Pnm 

5.4010. 

ai 

5.1610.03  1 

5.4510. 

03| 

5.7210. 

m\ 

5.2510.03 

5.5010. 

ai 

5.3710.03  1 

5.6610. 

5.9010. 

m\ 

5.3010. 

5.7010. 

5.4010. 

m 

5.1310. 

m 

5.3510. 

m 

5.0710. 

m 

5.3510. 

ai 

5.2110. 

mi 

5.1210.03 

5.3510. 

5.1710. 

m 

5.0710. 

m 

5.0610. 

m 

5.6810. 

a 

5.3410.03 

5.4210. 

m 

5.1410. 

m 

4.9510 

m 

5.0610 

m 

4.9310 

a 

4.8810.04 

5.1810 

4.8510.03 

5.0010 

4.8310 

m\ 

5.3610 

ai 

5.6010 

ai 

5.2010.03 

5.6310.03 

5.0010.03 

Table  5.  GLM  Event  Magnitudes  (Coni) 


Pt 


5.2010.03  I  5.0010. 


4.9410.03  4.6210.03  5.0010.25 


5.3010.03  5.0710. 


5.5110.03  5.2410.03  I  5.5010.25 


5.0010.03  4.7310.03  5.4510.25 


5.3010.03  5.0610.03  5.2910.25 


5.1910.03  4.9710. 


5.4610.03  5.2110.03  I  4.8710.25 


5.7210.03  5.5310.03  5.4210.18 


5.1610.03  4.9010. 


5.5010.03  5.2310. 


5.2310.03  4.9710.03  I  5.2010.25 


4.9610.03  4.7310.03  4.9410.14 


5.0110.03  4.7410.03  5.0710.18 


4.9010.03  4.6310.03  5.0810.18 


5.4710.03  5.2310.03  5.5710.25 


5.2110.03  5.0410.03  5.0910.25 


5.0510.03  4.7410.04  5.5710.25 


4.9810.04  4.7210.04  56110.25 


5.1710.03  4.9610.03  4.361025 


4.9610. 


4.8710.03  4.6310.03  5.0810.18 


4.8810.03  4.5710.03  5.1010.18 


5.4610.03  5.1310.03  5.3810.25 


5.2110.03  4.8610.03  5.5310.25 


5.1810.03  5.0010.03  5.2710.14 


4.8710.03  4.6410.04 


4.6810.04  4.4510.04  5.0310.18 


4.8610.03  4.6210.04  4.9510.14 


4.7610.03  4.4910.04  4.7210.25 


4.6310.04  4.4710.04  4.3810.18 


4.8710.03  4.6910.03  4.401025 


4.6710.03  4.4210.03  4.6610.18 


4.8110.03  4.4910.03  4.6610.12 


4.6410.03  4.3910.04  5.4610.25 


5.1410.03  5.0210.03  5.661025 


5.3710.04  5.1710.04  5.6110.14 


5.0310.03  4.8110.03  4.9910.18 


5.4910.03  5.2010.03  5.7410.18 


4.8610.03  4.7010.03  4.701025 


5.4010.12 


5.8210. 


3.6710.08 


3.6510. 


4.2110.05 


4.3410.04 


3.8510.08 


3.8910.10 


4.9810.04 


3.5110.17 


3.8210.09 


3.5810.12 


3.3510.08 


2.7110.10 


2.8510.10 


4.0510.05 


4.1610.05 


3.6010.04 


3.8410.04 


3.6110.06 


Table  S.  GLM  Event  MaonNudes  (Coni) 


Event 


DM781031 


DM781129 


DM790506 


DM790531 


DM791018 


DM800410 


DM800731 


DM810630 


DM810717 


DM820219 


DM820921 


DM831129 


DM831226 


DM840415 


DM870226 


DM870506 


DM670606 


DM870717 


DMe81123 


DM890217 


DM891004 


NZ641025 


NZ661027 


NZ671021 


NZ681107 


NZ681014 


NZ701014 


NZ710927 


NZ720828 


NZ730912 


NZ740829 


NZ750823 


NZ751021 


NZ760929 


NZ761020 


NZ770901 


NZ771009 


NZ780810 


Event 

Site 

NZ780927 

NN2 

NZ790924 

NNZ 

NZ791018 

NN2 

NZ801011 

NNZ 

N2811001 

NNZ 

NZ821011 

NNZ 

N2830618 

NNZ 

NZ830925 

NNZ 

NZ841025 

NNZ 

NZ870e02 

NNZ 

NZ880aU7 

NNZ 

NZ881204 

NNZ 

N2901024 

NNZ 

FS620501 

Sahara 

FS630318 

Sahara 

FS631020 

Sahara 

FS640214 

Sahara 

FS641128 

Sahara 

FS650227 

Sahara 

FS651001 

Sahva 

FS651201 

Sahara 

FS660216 

Sahara 

Table  5.  0114  Even!  Magnitudes  (Coni) 


S.S4±0.04 


$.6410. 


5.6510.03 


$.6110  03 


5.731003 


5.4010.03 


.03 


5.5710.03 


5.6710  02 


5  6010.03 


54310 


$6210.04 


5  6010 


5321007 


05 


56S10 


4  4210  17 


.12 


5.6410. 


4.5110  17 


53110  05 


50210  05 


5.3110  04 


54510  03 


5.4410.03 


54010.03 


5.5110  03 


.03 


5.4310.03 


5.4310  03 


54310  03 


52810.06 


.07 


.06 


44310.12 


$.0710.04 


5.2710. 


5.1810. 


5.1710  03 


5.2410  03 


5.0210.03 


5.2610.03 


6.3410.14 


6  7810  18 


6.1810.14 


6.5010.16 


6.0010.1 1 


5.9310.12 


5.7910.21 


52210  02 


5.1910  03 


6.0310.11 


62310  12 


.21 


.15 


6.1010.21 


5.9710.12 


52710.03 

5.1110  03 

6  0710.12 

.12 

5.2110  03 


5.0110. 


4.9010.07 


4.5710.05 


61910.12 


5.9210.12 


5.4110 

5.0710. 

.03 

5.35 

HD 

4.4510  12 


03 


Ms 


4.4010.04 


4.3410.04 


4.0110. 


4.1810. 


4.4210. 


4.0710.04 


4.3410.03 


.06 


4.4010.04 


4.3310. 


.04 


4.4510.04 


4.3310.06 


4.6710. 


VARIOUS  MAGNITUDES  OF  EVENT  SH8B091 4.pmax 


0  SO  100  ISO  200  2S0  300  3S0 


•  Bignal,  r  :  noiM.  *  :  cHp 


ngur*  1.  ScHiT  pM  of  3  tfitaroni  typoo  of  oMon  0%  o  for  Bolapon  ovoni  880914  (JVE).  Tho  26  good  rooorcfngo 
and  1  oip  or*  ohowm  «HM)  Mod  drcfn  otkI  upurard  orrow,  roopocdyofy.  Tho  raw  olalion  i%  •  (lop)  havo  a  atandard 
daalaten  of  0.30  itiu.  Appfyl)^  >Wnary”  aUdon  oorrocdona  radbooa  ttw  acaltar  to  0.21  lau.  Applying  tho  pro- 
poaad  “aaoondary**  oowaciona  to  oouni  tor  too  potfi  oflocto  rodboao  too  aattor  fudhar  down  to  0.14  m.u. 


MB 


VARIOUS  MAGNITUDES  OF  EVENT  SH891 01 0.pmax 


•  :  Y  ;  notM,  ^  ;  clip 


Rgur*  2.  Scatttr  plot  of  3  dHforonl  typo*  of  •Ution  a  for  Balapan  ovarH  891019. 
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VARIOUS  MAGNITUDES  OF  EVENT  DM71 0425.pmax 


•  ;  signaU,  Y  :  noiM,  *  ;  clip 


FIguraS.  Scatter  plot  of  3  dWaranttypMof  station  for  0«gt(*n«v*nt  71 0425. 
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VARIOUS  MAGNITUDES  OF  EVENT  DM761 230.pmax 


•  :  signal,  r  :  noiss,  *  :  clip 


Hgur*  4.  Scatt«r  p)ol  of  3  dHtoronl  typos  of  station  s  for  Dsgolsn  ovsnt  761230. 
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VARIOUS  MAGNITUDES  OF  EVENT  NZ691 01 4.pmax 


I 

I 

I 

I 

I 


Azimuth  (degrees) 


Tmst  SiU  :  NNZ 


•  :  sigrMl,  r  :  noism,  ^  :  clip 


*  Rgur*5.  Scattw’plotal3cfiff«’0n(typMoftUttionn%«forNovayaZ0nnlya*v«nl691O14. 

I 

I 
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IV.2  Rec«iv«r  and  Path  Effecta  on  Salamic  Wavaa 


Along  wHh  the  event  values,  the  station  terms  and  the  path  terms  are  generated  by  "rrHjgtm'at 
one  single  inversion.  Table  6  lists  the  station  corrections  arxj  path  correctiorrs  for  explosions  in  rWne 
nuclear  test  sites.  Note  that  the  station  terms  are  applicable  to  other  source  regions  of  the  \wor1d  as  wel. 
Applying  these  path  and  station  corrections  to  any  irxfividLial  explosions  wrouM  yield  a  reduction  in  the 
fluctuational  variation  of  station  magnitudes  with  a  factor  ranging  from  1.2  to  3.  Most  Novaya  Zemlya 
events  have  a  typical  reduction  factor  of  2,  as  noted  by  Jih  and  Wagner  (1992a). 

Figure  7  shows  our  receiver  terms  which  are  inferred  jointly  along  with  the  source-size  estimates 
and  path  terms  from  the  worldwide  explosions.  The  receiver  corrections  derived  with  our  approach 
rrratch  the  average  tectonic  structure  urvlerrreath  each  station  very  well,  maMy  due  to  the  broad  cover¬ 
age  of  azimuths  at  each  station.  Generally  speaking,  the  station  terms  are  positive  in  shield  regions 
such  as  Australia,  Canada,  Irrdia,  and  Scandinavia,  and  they  are  negative  in  the  east  Africa  rift  valleys, 
mid-ocean  ridges  (e.g. ,  Iceland  and  Azores  Islands),  island  arcs  (e.g. ,  Indonesia,  Japan,  and  Taiwan), 
arxl  Himalaya  Mountain  Rarrges  (Chaman  Fault,  rwrthem  Irrdia,  Nepal,  and  Burman  Arc).  Solomon  and 
Toksoz  (1970)  and  many  other  studies  (e.g. ,  Everrrden  and  Clark,  1970:  Booth  et  al. ,  1974)  found  that 
for  stations  in  the  U.S.,  the  attenuation  is  higher  between  the  Rockies  and  Cascades,  and  in  the 
northeastern  U  S.  This  pattern  is  also  observable  in  Figure  7  (see  also  North,  1977).  As  North  (1977) 
put  it,  it  is  gratifyirrg  that  a  simple  parameter  such  as  mu  can  be  utized  to  reveal  the  tectonics.  It 
should  be  noted,  however,  that  our  empirical  station  terms  also  irx:lude  the  effect  due  to  the  crustal 
amplification  if  such  local  site  effect  is  shared  by  all  ray  paths  from  different  test  sites  to  a  particular  sta¬ 
tion.  This  could  be  the  reason  of  a  few  outliers  such  as  HNR  (Honiara,  Solomon  Islands),  PMG  (Port 
Moresby,  East  Papua  New  Guinea),  arKl  RAB  (RabajI,  New  Britain),  which  do  not  show  negative  station 
terms  as  would  be  expected  from  the  strong  seismicity  in  that  region  (cf.  Figure  7).  Another  possbie 
reason  is  that  these  stations  have  relatively  poorer  azimuthal  sampling  In  our  data  set,  and  hence  the 
station  bias  at  these  three  stations  is  rwt  well  constrained.  The  minor  discrepancy  between  the  deter¬ 
ministic  corrections  by  Marshall  et  al.  (1979)  arxl  our  empirical  corrections  could  be  due  to  the  same 
reason.  Note  that  HNR,  PMG.  and  RAB  also  show  positive  station  terms  in  another  study  using  a 
different  data  set  (Jih  et  al..  1993). 

Figures  8  through  16  show  the  map  of  the  ”pure  path  effect"  (fop)  and  the  combined  station 
amplification  (bottom)  (defined  as  the  sum  of  the  receiver  term  and  the  path  effect)  for  all  9  test  sites 
used  in  the  GLM  inversion.  The  Semipalatirrsk  test  range  is  divided  into  four  test  sites  including 
Degelen  Mountain  [Deg]  and  three  subregiorrs  as  defined  in  Ringdal  et  al.  (1992):  southwestern 
Balapan  (BSWj,  northeastern  Balapan  [BNE],  and  the  transition  zone  |BTZ]  between  BSW  and  BNE. 
The  path  term  at  each  station  can  be  regarded  as  the  azimuthal  variation  (towards  the  various  source 
regions)  relative  to  the  averaged  station  ampUfication.  An  important  obsenration  is  that  all  these  five  test 
sites  exhibit  very  different  azimuthal  and  radial  amplitude  variations.  Events  from  Degelen  and  BTZ  are 
systematicaily  enhanced  in  the  western  U.S.  arxl  reduced  in  the  eastern  U.S.,  whereas  events  from 
BSW  and  BNE  are  enhanced  in  essentiafiy  the  whole  of  U.S.  Such  highly  dkection-dependenl, 
distance-dependent,  and  site-dependent  patterns  of  the  amplitude  fluctuation  could  be  a  diagnostic  for 
the  path  effects  in  the  proximity  of  the  test  sites.  Back  projectons  (e.g. ,  Lynnes  and  Lay,  1990)  of  the 
rOb  residuals  onto  the  upper  mantle  and  the  lower  crust  reveal  that  simiar  mi,  residuals  come  into 
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alignment  in  several  regions  partitioned  by  known  geological  features  (Jih  and  Vt/agner,  1991).  Degelen 
events  in  the  western  U.S.,  and  SW  Balapan  everts  at  western  European  stations  must  pass  through 
the  area  between  Chinrau  fault  and  Chingiz-Kat>a  shear  zone.  Alt  these  paths  show  positive  n%  residu¬ 
als.  The  area  north  of  ChirKau  fault  might  have  some  complex  features  that  result  in  negative  mean  n% 
residuals.  Paths  from  NE  Balapan  to  North  America  and  many  continental  European  statiorw  must  cross 
this  area  or  even  travel  along  the  Chinrau  fault  before  entering  the  deeper  mantle,  arxl  hence  the  com¬ 
plexity  in  the  waveforms  is  irtevitable.  it  seems  that  the  mean  m^-Lg  separation  of  0.07-0.17  m.u. 
(e.g. ,  Rir^al  and  Hokiand,  1987;  Ringdal  et  at..  1992;  Richards  ef  at..  1990;  Section  V.2  of  this 
report)  between  the  NE  and  SW  subregions  of  Balapan  could  be  due  in  part  to  the  path  effects,  in  addi¬ 
tion  to  the  difference  of  source  medium  postulated  previously  by  Marshall  at  at.  (1984).  A  detailed  dis¬ 
cussion  on  the  seismic  variabiity  within  Balapan  test  site  is  given  in  a  later  section.  Path  effects  can  also 
explain  why  the  SW  Balapan  waveforms  tend  to  be  more  complex  at  YKA  than  those  recorded  at  WRA, 
EKA,  and  GBA  arrays. 

The  initial  P  waves  from  the  three  adjacera  test  sites  have  virtually  the  same  incidera  angle  at 
each  teleseismic  station,  and  anything  in  common  across  all  events  (such  as  the  crustal  amplification  as 
weU  as  the  upper  mantle  attenuation  underneath  the  receiver)  would  have  been  lumped  into  the  constant 
station  term.  Thus  the  station  residuals  averaged  over  all  events  from  the  same  test  site  would  correlate 
very  little  with  the  receiver.  Instead,  they  should  reveal  more  site-dependent  information  about  the 
focusing/defocusing  pattern  underneath  E.  Kazakhstan. 

The  iarge.si  and  most  pronrtirtent  fault  in  the  region  is  the  southeast-trendirtg  Chingiz  right-lateral 
strike-slip  fault  that  passes  about  10  km  southwest  of  Degelen  Mountain  and  right  across  the  Murzhic 
test  area  (Rodean,  1979;  Bonham  ef  at. ,  1980;  Leith,  1987b).  Soviets  reported  that  this  fault  has  a  very 
steep  dip,  which  is  consistent  with  its  linear  expression  over  large  distance  as  seen  on  Landsat  imagery 
(Bortham  at  at.,  1980).  A  distinct  fault-line  scarp  is  developed  along  much  of  the  oldest  metamorphic 
rocks.  Chingiz  Fault  extends  for  a  total  length  of  aybout  700  km.  Soviet  reports  postulate  that  this  faul 
extends  down  to  the  boundary  of  the  granite  layer  of  the  crust  arxl  possbiy  into  the  upper  mantle. 

The  inferred  path  terms  for  Novaya  Zemtya  explosions  are  very  consistent  with  those  in  Jih  and 
Wagner  (1992a).  Jih  and  Wagner  (1992a)  also  compared  the  path  terms  with  the  travel-time  residuals. 
Their  results  irxlicate  that  paths  from  the  northern  test  site  in  Novaya  Zemlya  to  stations  in  North  Amer¬ 
ica  have  systematically  faster  arrivals  and  smaler  amplitudes,  suggesting  a  profound  defocusing  effect 
on  the  first  arrivals;  wNIe  stations  in  Ireland,  Scotland,  Spain,  Bangladesh,  northern  India,  Pakistan, 
Korea,  arxi  Kenya  report  slow  arrivals  and  large  amplitudes,  suggesting  a  focusing  effect.  Amplitudes  for 
paths  to  Greenland,  Iceland,  Alaska,  Turkey,  Germany,  Luzon,  Zimbabwe,  Italy,  Pueto  Rico,  Ethiopia, 
and  Hawaii,  however,  seem  to  be  controlled  by  the  anelastic  attenuation  with  slow  rays  also  associated 
with  small  amplitudes,  and  fast  rays  associated  with  large  amplitudes. 
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Table  6.  Receiver  and  Path  Terms  for  Eurasian  Nuclear  Test  Sites 


Station  Term  [S]  Path  Terms  IF] 


Code  I  Rev’  I  BNE  1  BSW  |  3TZ  |  Deg  |  NNZ  |  PMA  |  RNA  |  Sahara  |  YFT 


-.122  I  -.399 


.193 


.053 


-.101 


-.164 


-.140  -.049 


AAE 


AAM 


AFI 


AKU 


ALE 


ALQ 


ANMO 


ANP 


ANTO 


AQU 


ARE 


ARU 


AS023 


AS1027 


ATL 


ATU 


AZ043 


BAG 


BC045 


BCAO 


BOF 


BE051 


BEC 


BER 


BHP 


BJI 


BKS 


BLA 


BLC 


BM068 


BOCO 


BOD 


BOG 


BUL 


CAR 


CHG 


CHS 


CHTO 


1)  th*  MUion  biM  wNdt  naadi  ID  b*  ooffcud  (in  addiban  id 

2)  BSW  .  SW  (ubtilD.  BUlw;  BNE  .  NE  DubNlt.  BNiyw: 
PMA  .  PahuM  Maaa,  NTS;  YFT  .  Yucca  Flat,  NTS:  RNA  -  f 


>  tw  paih  affacQ. 

i;  BTZ  -  trandbon  zona.  Bat 

RaWar  Maaa,  NTS;  Sahara 


an;  Dag  •  Dagalan  MDuntain;  NNZ  -  northam  Wand,  Novaya  Zanbya; 
Ahaggar,  Franch  Sahara. 
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Table  6.  Receiver  arvJ  Path  Terms  for  Eurasian  Nuclear  Test  Sites  (cont) 


Station  Term  (S)  Path  Terms  [F] 


Code  I  Rev’  I  BNE  |  BSW  |  BT2  I  D«0  |  NNZ  I  PMA  I  RNA  I  Sahara  I  YFT 


CM074 


CMB 


CMC 


COL 


COP 


COR 


CTAO 


DAG 


DAV 


DUG 


EIL 


EPT 


ESK 


FCC 


FFC 


FL099 


FLO 


FRAO 


FRB 


FVM 


FX118 


GBABI 


GDH 


GEO 


GIE 


GOL 


GRFO 


GSC 


GUA 


GUMO 


GWe 


HA122 


HiA 


HKC 


HLW 


HNR 


IB132 


IL1154 


-.029 

-.192 

.004 

-.361 

1)  th«  tMion  bias  which  nasds  Id  ba  oonadad  (in  addition  id  tfia  path  aflacQ. 

2)  BSW  -  SW  subsits,  Balapan;  BNE  -  NE  aubailt.  Balapan;  BTZ  -  traniibon  zona.  Balapw;  Dag  -  Dagalan  Mountain:  NNZ  -  nonham  Mwid,  Noviya  Zamlya; 
PMA  >  Pahma  Masa.  NTS;  YFT  .  Yucca  Flat,  NTS;  RNA  .  RaMar  MMa,  NTS:  SNiva  -  Ahi«gw.  French  SNiwa. 


Table  6.  Receiver  and  Path  Terme  for  Eurasian  Nuclear  Test  Sites  (coni) 


Station  Term  [S]  Path  Terms  [F] 


Code  I  Rev'  I  BNE  BSW  BTZ  Deg  NNZ  PMA  RNA  Sahara  YFT 


.027 


IM182 


INK 


1ST 


JAS 


XT 


JER 


X191 


KAAO 


KBS 


KEV 


KIP 


KOD 


KON 


KRK 


KS213 


KTG 


KU238 


LA251 


LEM 


LHC 


LON 


LOR 


LPB 


LPS 


LUB 


LY252 


L2H 


MA265 


MAIO 


MAJO 


MAL 


MAT 


MBC 


MDS 


MI269 

.095 

.247 

.221 

.205 

-.175 

MSO 

-.247 

-.154 

MUN 

.257 

-.004 

-.107 

-.026 

.001 

N04C2 

.511 

.744 

.645 

.315 

1)  th«  ituian  bias  wtvdi  naadt  B  ba  corraciad  (in  addWon  B  tf«a  pa 

2)  BSW  .  SW  MibaiB.  Balapan:  BNE  -  NE  aubaiB,  BN^an;  BTZ  • 

PMA  .  PahuB  Maaa,  NTS;  YFT  •  Yucca  Flat.  NTS;  RNA  .  Rainiar 


I  bananon  zona.  Ban 
Maaa,  NTS;  Stfirn 


an:  Dag  >  Dagolan  Mountain;  NNZ  -  northam  Wand,  Novaya  ZamVa; 
Ahaggar,  French  Sahara. 
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Station  Term 


Code 


Table  6 


ISI 


Receiver  and  Path  Terms  for  Eurasian  Nuclear  Test  Sites  (cont) 


Path  Terms  [F] 


BNE  BSW  BTZ  I  Deg  I  NNZ  I  PMA  |  RNA 


Sahara 


NAO 


NAI 


NAT 


NOI 


NNA 


NOR 


NRAO 


NUR 


NVS 


NWAO 


OBN 


OGO 


OTT 


OXF 


PDA 


PEL 


PHC 


PTO 


CUE 


RAB 


RAR 


RES 


RSCP 


RSNY 


RSON 


RSSD 


SCH 


SCP 


SOB 


SEO 


SH276 


SHA 


SHI 


SHK 


SHL 


1)  ih*  UMkm  biM  which  nwdi  id  b*  oorrKMd  (In  addibon  id  tw  path  aHacQ. 

2)  BSW  -  SW  Kibaila,  Balapan;  BNE  -  NE  tubtila,  Baiapw;  BTZ  -  trwwiton  nna.  Balapai;  Dag  -  Dagaltn  MDuniain:  NNZ  • 
PMA  >  Pahuia  Mata.  NTS;  VFT  .  Yuoea  Flat.  NTS;  RNA  -  RaMar  Maaa,  NTS;  Sahara  -  Ahaggar,  Franch  Sahara. 


rmiham  iaiand,  Naaaya  Zandya; 
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Table  6.  Receiver  and  P^h  Terme  for  Eurasian  Nuclear  Test  Sites  (cont) 

Station  Term  [S] 

Path  Terms 

IF] 

Code 

BNE 

BSW 

BTZ 

Oeg 

NNZ 

PMA 

RNA 

Sahara 

YFT 

SJG 

.177 

-.586 

-.368 

-.297 

SLR 

-.482 

■■■I 

■■■ 

SNG 

-.153 

-.126 

.033 

S0286 

-.068 

-.116 

-.101 

-.079 

.264 

STJ 

.388 

_ 

-.011 

STU 

.071 

K2I 

.015 

.175 

-.262 

-.295 

-.054 

-.282 

TAB 

.520 

TATO 

-.322 

.189 

-.130 

.062 

.345 

TLY 

-.289 

-.034 

TOL 

.131 

-.213 

-.181 

-.043 

-.144 

.502 

.072 

.087 

TRI 

-.106 

-.053 

.308 

.231 

-.473 

-.159 

-.649 

TRN 

.008 

.341 

-.016 

-.121 

-.280 

-.175 

TT297 

-.261 

.039 

.002 

.028 

-.095 

TUC 

-.051 

-.383 

-.023 

.210 

TUP 

-.369 

.142 

-.230 

UME 

.110 

.297 

.356 

-.071 

-.188 

-.452 

-.424 

.034 

UPA 

-.389 

.020 

-.204 

UZH 

-.135 

-.330 

-.125 

-.191 

-.211 

.599 

VAL 

.031 

Hd 

.079 

.179 

-.354 

-.325 

-.215 

WES 

-.251 

.113 

.037 

-.414 

WMQ 

.375 

-.031 

WRAR1 

.398 

-.016 

.102 

-.193 

-.131 

YKAB5 

-.032 

.345 

.401 

.155 

-.231 

-.211 

-.241 

YKC 

-.009 

-.068 

20BO 

-.101 

.039 

1)  ih«  tiMian  bias  ««hich  naad*  lo  b*  corraobd  (in  adtfiion  b  tfw  path  aflacQ. 

2)  BSW  .  SW  wbaiia.  Balapan;  BNE  •  NE  aubaila.  Balapan:  BTZ  .  ttmibon  lona.  Balapan;  Dag  -  Dagalan  Mountbrc  NNZ  -  nodbam  iaiwid.  Novaya  Zandya; 
Pli4A  .  Pahuia  Maaa.  NTS:  YFT  .  Yucca  Flat.  NTS;  RNA  .  Ranar  Mata.  NTS;  SN«m  -  Ah^g».  Frarch  Stfaaa. 
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MEAN  STATION  AMPUFICATION  ON  MB 


Potar Khmithal •quMMtnt prolacUon,  71.00,60.00 


FIgura  7.  Map  of  ir^  station  tarma  Marrad  from  a  GLM  invaraion  solvat  for  3,306  unknowns  wKh  36,211  Nnoar 
aquationa.  Tha  high  oorralation  batwaan  tha  taoUxtic  typo  ond  tha  station  tarma  suggstts  that  thaas  ampiricai 
oorractiona  do  raflact  tha  uppar  mantla  oonditiona  undamaath  tha  racaivara.  Tha  4  darkanad  atara  rapraaani  aoma 
of  tha  nudaar  taal  aNaa  uaad  in  thia  atudy. 
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Path  Effect 
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Flgura  8.  Tho  map  showing  tho  "pura  propagation  affocT  (lop)  and  ttia  combined  station  ampificalion  (bottom) 
defined  as  the  sum  at  the  receiver  term  (Figure  7)  and  the  path  effect  for  NE  Bab;  sn  explosions. 
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Figure  9.  and  the  combined  station  ampiification  (bottom)  defined  as  the  sum  of  the  receiver  term  (Figure  7)  and 
the  path  effect  for  SW  Balapan  test  site. 
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FIgura  11.  Tha  map  showing  tha  "pura  piopagatbn  affacT  (top)  and  tha  oombinad  station  amplification  (bottom) 
dafinad  as  tha  sum  of  tha  racaivar  tarm  (Figura  7)  and  tha  path  affect  for  Dagalen  test  site. 
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FIgura  12.  Tho  map  showing  tha  "pure  propagation  affacT  (top)  and  tha  oombinad  station  ampiirication  (bottom) 
dafinad  as  tha  sum  of  tha  racaivar  farm  (Figure  7)  and  tha  path  attact  for  Northarn  Novaya  Zamlya- 
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FIgura  13.  The  map  showing  the  “pure  propagation  effecT  (top)  and  the  combined  station  ampUTication  (bottom) 
defined  as  the  sum  of  the  receiver  term  (Figure  7)  and  the  path  effect  tor  Pahute  Mesa. 


Path-fRecelver  Path  Effact 


STATION  AMPUHCATION  OF  MB  FOR  RAINIER  SHOTS 
Polar  azimuthal  •quMIstaiM  projection,  *116.40,37^ 


Figure  14.  The  map  showing  the  "pure  propagation  affocT  (top)  and  the  combinod  station  amplification  (bottom) 
defined  as  the  sum  of  the  receiver  term  (Figure  7)  and  the  path  effect  for  Rainier  Mesa. 


Path4>Raceiv6r  Peth  Efftct 


Ui 


X 

0L500 

X 

OJMO 

X 

0.167 

X 

0.125 

X 

0.100 

X 

0.083 

a 

•0i»3 

o 

•0.100 

o 

•0.125 

O 

•0.167 

O 

•0.250 

o 

•0.500 

X 

0.500 

X 

0.250 

X 

0.167 

X 

0.125 

X 

0.100 

X 

0.083 

a 

•0.083 

o 

•0.100 

o 

•0.125 

O 

•0.167 

o 

•0.250 

o 


•0^ 


STATION  AMPURCAT10N  OF  MB  FOR  SAHARA  SHOTS 
Polar  azimuthal  aquidiatant  projaction,  5.05 , 24.00 


FIgura  16.  Tha  map  showing  tha  "puia  propagation  affacT  (lop)  and  tha  oombinad  station  amplification  (bottom) 
dafinad  as  tha  sum  of  tha  racaivar  tarm  (Figura  7)  and  tha  path  affact  for  Franch  Sahara. 
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V.  COMPARISON  OF  VARIOUS  MAGNITUDES 
V.1  Comparison  of  Various  GLM  Magnitudes 

Table  7  compares  and  n\(Pt)  relative  to  mtiPa)  at  several  nudear  test  sftes. 


Table  7.  rnbiP^m^  and  mi,(P*)  vs.  n^(Pa)  | 

Test  Site 

m^fP m«*)~m^(P i)  # 

mb{Pb)-mb{Pb)  # 

PMA 

0.569±0.016  58 

0.28010.010  58 

RNA 

0.590±0.047  15 

0.22210.042  15 

YFT 

0.601±0.015  141 

0.23110.013  141 

NTS 

0.59210.011  214 

0.24310.009  214 

Deg 

0.44610.008  98 

0.24810.007  98 

BNE 

0.47810.011  30 

0.25510.009  30 

BT2 

0.49510.014  19 

0.28210.007  19 

BSW 

0.48110.005  48 

0.28110.004  48 

KTS 

0.46410.005  195 

0.26110.004  195 

NN2 

0.40710.012  30 

0.22610.007  30 

Sahara 

0.42210.050  7 

0.12510.053  7 

Before  Bocharov  at  al.  (1989)  published  the  yields  and  other  source  information  of  98  historicai 
Soviet  events,  several  attempts  had  been  made  to  investigate  the  characteristics  of  cratering  explosions 
in  that  region.  For  instance,  McLaughlin  at  al.  (1985)  studied  the  ratio  of  the  P»  phase  and  Pmn 
phase  of  presumed  Balapan  contained  and  cratering  explosions  by  comparing  the  WWSSN  station 
mb's.  The  rrwtivation  was  that  the  logarithm  of  amplitude  ratio  of  PmuiP,  of  event  SH650115  was 
signiTicadly  smaller  than  other  presumed  contained  explosions  in  the  vidnity.  Assuming  the  phase  P.  is 
unaffected  by  the  influence  of  the  non-linear  tree-surface  interference,  then  an  adjustmeit  to  the 
tnt  (PiMx)  should  be  able  to  convert  that  to  a  contained  explosion  of  the  same  yield.  McLaughlin  ef  aL 
(1985)  concluded  that  a  correction  between  0.17  and  0.27  is  needed  for  this  conversion,  assuming  a 
yield  of  125KT. 

Der  at  al.  (1985)  deconvolved  four  contained  and  the  cratering  Balapan  events  SH650115 
recorded  at  EKA.  and  then  they  convolved  the  Green’s  fundions  with  an  appropriate  attenuation  opera¬ 
tor  as  well  as  the  source-time  function  of  various  yields  of  interest.  By  comparing  the  phases  P,  and 
Pfflu  d  the  synthetics,  they  obtained  a  cratering-to-contained  correction  of  0.15,  0.15,  and  0.18  at  60. 
125,  and  300KT,  respectively- 

Day  at  al.  (1986)  did  a  theoretical  study  with  nonlinear  source  calculations  to  account  for  coupling 
variations  with  depth.  Their  results  are  summarized  as  follows: 

mb  (contained)  -  (cratering)  s  -  (0.1  to  0.15)  direct  coupling 

*  (0.1  10  0.25)  r  rrface  interaction  effects 
-  (0  to  0.15)  total  bias. 
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Based  on  46  Balapan  explosions  recorded  at  EKA,  Ringdai  ef  al.  (1992)  derived  a  value  of  0.75 
as  their  niean  toQ(P„^/P, )  across  the  EKA  array  using  the  same  techniques  as  used  in  McLaughlin 
et  al.  (1985).  The  cratering  event  SH650115  had  m/,(Pmu)  *  ^(^a)  -  0.62  at  EKA,  and  hence  they 
apply  a  correction  of  5.87  -t-  (0.75  -  0.62)  -  6.00  for  a  hypothetical  contained  explosion  with  equivaient 
yield. 

For  event  SH650115,  our  /n«(Pmu}  *  ■  0-36  (cf.  Table  5),  which  is  about  0.10  to  0.14 

m.u.  lower  than  the  average  /n^(Pmu)  -  mt,{Pa)  value  of  0.464  (KTS)  or  0.495  (BTZ,  where  this  event 
is  located)  shown  in  Table  7.  Thus  statistics  in  Ta^e  7  suggest  that  our  correction  is  in  better  agreement 
with  those  derived  by  Ringdai  et  al.  (1992)  and  Day  et  al.  (1986).  It  is  also  in  agreement  with  the  result 
in  Jih  et  al.  (1993)  using  primarily  WWSSN  data. 

The  cratering-to-oontained  conversions  cited  above  typically  require  extra  information  about  the 
general  behavior  of  contained  explosions  in  the  same  source  region.  For  the  purpose  of  estimating  the 
yield  of  a  cratering  shot  in  an  isolated  region,  directly  using  P,  could  be  a  much  easier  approach  (Jih 
et  al.,  1993). 

Without  access  to  official  yields,  the  simplest  way  to  examine  the  transportability  of  scaling  laws 
between  different  test  sites  would  be  to  relate  the  seismic  phase  of  interest  to  the  generally  more  tran¬ 
sportable  Ms .  Table  8  gives  the  results  of  regressing  mi,{Pmu) ,  f^biPb)  •  and  mb(P,)  on  Ms  (LR)  at 
major  test  sites  with  the  slope  fixed  at  1.  It  is  interesting  to  note  that  rnb(Pmui-Ms  at  KTS  is  0.55  m.u. 
larger  than  that  at  NTS,  and  that  NNZ  has  an  intermediate  value  between  those  of  KTS  and  NTS. 

Tables  9  and  10  compare  against  Lg  and  P„,  respectively.  The  mg-Lg  bias  between  KTS 
and  NTS  is  0.38  m.u..  which  is  in  good  agreement  with  the  comnx>nly  quoted  ntg  bias  of  0.35  m.u.  at 
150  kt  level  between  KTS  and  NTS  (e.g..  Jih  et  a/.,  1993;  U.S.  Congress/Office  of  Technology  Assess¬ 
ment,  1988;  and  many  others).  In  his  pioneering  Lg  study,  Nuttli  (1988,  1986ab,  1987)  suggested  that 
Lg  excitation  from  explosions  may  be  independent  of  the  source  material,  arxl  hence  mgiLg)  can  be 
used  to  calibrate  which  is  susceptible  to  the  upper-mantle  attenuation.  Nuttli  derived  his  ntig^Lg)- 
yield  relationship  with  NTS  explosions  and  applied  the  same  formula  to  estimate  the  yields  of  Semipala- 
tinsk  explosions,  and  the  resulting  yield  estimates  are  very  good.  This  has  led  to  numerous  follow-up 
research  in  recent  years  using  Lg  as  an  excellent  relafive  yield  estimator.  It  should  be  noted,  however, 
that  the  absolute  nib(^)  (or  Lg )  scale  is  still  subject  to  the  source  region  bias.  As  Table  9  indi¬ 
cates.  mb-ntgiLg)  at  NNZ  is  even  smaller  than  that  at  NTS.  The  same  mb(Lg)  formula  was  used  for  all 
test  sites,  arxf  the  path  corrections  we  used  are  either  based  on  Nuttli’s  calculation  or  based  on  a 
scheme  which  produces  path  attenuation  equivalent  to  what  Nuttli  would  give.  In  Nuttli's  earlier  Lg  work 
for  earthquakes,  the  time-domain  approach  was  use  in  determining  the  Lg  attenuation.  Assuming  the 
path  attenuation  parameters  obtained  with  the  coda-Q  method  is  appropriate  for  NNZ  explosions,  then 
Table  9  would  suggest  that  a  Lg  bias  of  approximately  0.3  rau.  is  necessary  to  bring  Lg  scaling  in 
alignmert  with  that  of  Ms  {cf.  Table  8).  That  is  to  say,  NNZ  mi,{Lg)  are  probably  somewhat  over¬ 
estimated. 
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Test  Site 


PMA 


RNA 


YFT 


NTS 


Deg 


BNE 


BTZ 


BSW 


KTS 


NNZ 


Table  8.  n\  (QLM)  vs.  Mg  (GLM)  at  Various  Sites 


0.395±0.031  53 


0.854±0.093  2 


0.52710. 


0.48110.019  139 


1.10810.041 


1.27210.040  26 


/J1^(P(^)  — LR,  i 


0.67610.02B  53 


1.11210.121  3 


023 

84 

0.79910.021  84 

1.14810. 

.95910.027  53 


1.46610.124  3 


84 


1.08310.019  140 


1.54510. 


1.74710.043  26 


1.71410.059  18 


1.67010.029  47 


1.63710.022  154 


1.26410.035  29 


1.19010.031  47 


1.17410.022  154 


0.85610.036  29 


060 

18 

1.50310. 

058 

0.75910.018  140 


1.35210.041  63 


1.53110.042  26 


18 


1.47110.031  47 


1.43610.022  154 


1 .08310.036  29 
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Table  9.  mg  (QLM)  vs. 


mg(P t)—Lg ,  # 


-0.79110.024  48 


-0.58510.054  15 


-0.87310.027  107 


-0.82510.020  170 


-0.25510.041  12 


-0.34610.039  15 


-0.28810.024  27 


-0.31010.017  74 


-0.62010.043  24 


Lg  (GLM)  at  Various  Sites 


mg(Pingg^  ~Lg ,  § 


-0.21510.020  48 


.032 


-0.25810.027  113 


-0.21710.019  183 


0.16710.035  12 


0.12710.033 


0.14510.038  15 


0.19710.024  27 


0.16310.015  74 


-0.23710.041  24 


-0.51210.023  48 


0.37010.035  18 


•0.63810.027  108 


0.57610.019  174 


-0.01510.037  12 


-0.08810.037  20 


15 


0.00010.023  27 


.016  74 


-0.40110.046  24 
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Table  10.  /n»(GLM)  vs 


mgiPs)-P„.i 


-0.51510.018  54 


-0.05510. 


P„  (GLM)  at  Various  Sites 


n\{Pb)-Pn.* 


0.23410.014  54 


-0.38510.018  207 


-0.33410.047  55 


-0.176ct0.123 


•0.562t0.110  6 


.45210.064  32 


-0.35110.037  113 


-0.14010.017  211 


.043  55 


0.08410.119  20 


-0.29110.113  6 


-0.176±0.063  32 


.035  113 


.012  54 


15 

0.16310.026  18 

0.52210.032 

138 

-0.14310. 

022 

139 

0.23110. 

022 

058  20 

-0.72110.058  20 

-0.52510. 

223 

145 


0.21810.017  221 


0.10310.043  55 


.123 


-0.09010.111  6 


0.02210.064  32 


113 
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V  J  n\*Lg  VartabiUty  wHhIn  Balapan  Tast  Sit# 

Marshall  et  al.  (1984)  found  that  explosiortt  in  the  northeast  and  southwest  portions  of  Balapan 
test  site  produce  distinctly  different  waveforms  when  recorded  at  the  UK  seismological  array  statiorw, 
suggestirx)  that  Balapan  test  site  can  be  subdivided  Mo  two  areas  characterized  by  different  geophysi¬ 
cal  properties.  Ringdal  and  Hokland  (1987)  find  that  this  pattern  is  persistently  present  whether  n\ 
based  on  worldwide  network  or  (Pooda)  of  NORSAR  is  used.  They  inferred  the  average  -Lg 
between  SW  and  NE  subregions  as  0.17  m.u.  In  a  foHow-up  study,  Ringdal  and  Fyen  (1988)  suggest 
that  there  appears  to  be  a  transition  zotm  between  the  NE  and  SW  subregbns.  Ringdal  at  al.  (1992) 
recomputed  the  SW-NE  bias  as  0.15  m.u.  with  101  Balapan  events  recorded  at  ISC  stations  and  NOR¬ 
SAR.  AKiiough  Ringdal  at  al.  (1992)  agree  that  the  possibilily  of  a  mb(Lp)  bias  contributing  to  this 
differerx»  between  SW  and  NE  cannot  be  entirely  njled  out,  they  propose  an  empirical  approach  to 
correct  for  this  bias  by  assuming  this  bias  is  solely  due  to  a  relative  bias  between  these  two  areas. 

We  followed  the  zoning  of  Rirrgdal  at  al.  (1992)  in  partitioning  Balapan  test  site  Mo  three  regions: 
southwest  (SW),  transition  zone  (TZ),  and  northeast  (NE).  Figure  17  shows  the  spatial  pattern  of  mg  - 
Lg  residuals  of  Semipalatinsk  explosions  based  on  Qeotech’s  /?%  values  and  RMS  Lg  values  reported 
at  NORSAR.  There  is  a  significant  difference  in  the  source  medium  across  the  CNnrau  fault  separating 
the  northeastern  and  southwestern  portion  of  Balapan  test  site,  as  reported  by  Ringdal  at  al.  (1992)  and 
Marshall  at  al.  (1984).  The  mean  mg  -Lg  bias  between  SW  and  NE  Balapan  is  about  0.07  m.u. 

Note  that  the  mg(Pmui-f-g  bias  of  0.08  m.u.  between  SW  and  NE  (cf.  Table  12)  is  signiflcantly 
smaller  than  that  of  previous  studies.  Regressing  the  RMS  Lg  furnished  by  Israelson  (1992)  and  our 
/Hg  on  the  yields  published  by  Bocharov  at  al.  (1989)  (and  Vergino,  1989)  shows  that  NE  explosions 
have  positive  Lg  residuals  and  negative  mg  residuals,  whereas  SW  explosions  show  the  opposite  trend 
(Figure  10  of  Jih  at  al.,  1993).  A  three-dimensionai  geological  model  of  the  Balapan  test  site  by  Leith 
and  Unger  (1989)  shows  a  distinct  difference  between  the  NE  and  SW  portions  of  the  test  site,  with  the 
granites  closer  to  the  surface  and  the  alluvium  thinner  in  the  southwest.  The  thicker  alluvium  layer  in  NE 
region  could  increase  the  waveform  complexity  and  reduce  the  magnitudes  measured  with  .  The 
first  motion  should  be  least  affected  by  this  factor,  however.  Jih  at  al.  (1993)  suggest  that  the  mg  -Lg 
bias  between  SW  and  NE  Balapan  can  be  tentatively  decomposed  into  several  parts: 

[I]  Difference  in  pP  between  SW  and  NE, 

[II]  Difference  in  mg  coupling,  l.a. ,  mg  (SW)  >  mg  (NE), 

[III]  Difference  in  Lg  coupling.  La. ,  Lg  (NE)  >  Lg  (SW), 

[iV]  Effects  due  to  the  station-station  correlation  structure, 

[V]  Effects  due  to  the  uneven  geographical  clustering  of  stations,  as  well  as  any  path  effect  which  is  not 
fully  accounted  for  through  the  network  averaging. 

The  bias  of  0.08  m.u.  for  mg{Pm^-RMS  Lg  between  SW  and  NE  Balapan  (Table  12)  is  essentialy  the 
sum  of  [I]  through  [III].  However,  if  the  time-domain  Lg  magnitude  is  used  instead,  the  bias  between  SW 
and  NE  Balapan  is  insignificant  (Table  11). 

For  ISC  data,  we  estimate  that  [V]  is  about  0.02  m.u.  if  derived  by  the  conventiortal  LSMF 
are  used.  When  /Hg  is  used,  this  term  is  eliminaled,  and  hence  a  smaller  mg  -Lg  bias  is  obtained.  [II] 
and  [III]  can  be  easily  illustrated  with  regressions  on  Bocharov's  published  yields,  as  explained  earlier. 
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There  are  only  a  handftil  of  Balapan  events  wih  published  yields  in  Bocharov  at  al.  (1989).  Howrever, 
the  5  large  historical  events  (lor  which  the  yields  were  exchanged  during  JVE)  can  also  provide  some 
supplementary  clue  in  support  of  our  postulated  hypotheses  [I]  through  [III].  The  yield  estimate  based  on 
Pmn  tor  two  (out  of  three)  historical  events  in  SW  subregion  (790804B  and  791223B)  is  larger  than  that 
based  on  P.  •  On  the  other  hand,  the  two  events  in  NE  subregion  (791028B  and  840526B)  have  a 
smaller  yield  estimate  based  on  P„^  as  compared  to  P« .  The  larger  bias  of  0.15  m.u.  that  Ringdal 
ef  al.  (1992)  obtained  with  (ISC)  could  have  been  slightly  “enhanced”  due  to  (IV]  and  [V].  The 
determination  procedure  presented  in  this  study  does  not  correct  lor  [IV]  either.  However,  the  cortribu- 
tion  of  inter-station  correlation  alone  is  believed  to  be  insignlicant  if  WWSSN  is  used. 

In  Figure  19  we  show  the  differertce  of  path  effects  between  BSW  arvJ  BNE  at  each  WWSSN  sta¬ 
tion,  which  is  a  measure  of  the  relative  bias  between  BSW  and  BNE  atong  each  path.  Positive  symbols 
represers  the  stations  where  BSW  everss  are  enhanced  relative  to  BNE  events.  K  the  raw  station  mag¬ 
nitudes  are  used  in  the  network  averagirtg  without  fully  accounting  for  such  path-effect  differeraial, 
significant  bias  (relative  to  the  magnitude)  will  be  present.  ISC  network  is  dominated  by  western 
European  stations,  and  hence  the  effect  due  to  [V]  would  be  more  severe  than  that  on  WWSSN. 

Nutty  (1987,  1988)  suggests  that  there  is  a  mi>  bias  of  about  0.2  m.u.  between  Degelen  and 
Balapan,  with  Degelen  explosions  having  even  larger  excitation  (relative  to  ).  We  do  not  see  such 
Degelen-Balapan  bias  with  Nuttli's  (Table  11)  or  RMS  Lg  measured  at  NORSAR  (Table  12). 


Table  11.  /flc  (GLM)  vs.  m„{Lg)  (NvnW)  at  Various  Sites 

Test  Site 

mb(Pi,) -Nutty.  # 

nii>(PmM)-Nuttli,  # 

mi,{Lg)-HiM 

PMA 

-0.894±0.054  16 

-0.63O±0.0S3  16 

-0.326±0.044  16 

-0.06210.042  13 

RNA 

-0.627±0.066  6 

-0.458±0.035  7 

-0.057l0.035  8 

-0.07110.062  9 

YFT 

-0.841±0.046  56 

-0.636±0.0^  57 

-0.27310.048  59 

0.16710.121  44 

NTS 

-0.836±0.036  78 

-0.619±0.037  80 

-0.26310.036  83 

0.09010.062  66 

BSW 

-0.406±0.020  26 

0.08010.020  26 

-0.10410.032  18 

BNE 

-0.415±0.026  21 

-0.150±0.026  21 

0.05910.026  21 

BT7 

-0.415±0.029  12 

-0.133±0.025  12 

0.07710.030  12 

-0.03610.059  9 

Deg 

-0.350±0.054  20 

-0.100±0.058  20 

0.10210.057  20 

_ 

KTS 

-0.395±0.017  79 

-0.125±0.017  79 

0.07910.017  79 

-0.08110.025  45 

NNZ 

-0.589±0.035  25 

-0.363±0.037  25 

-0.18810.032  25 

0.02110.024  20 

Table  12.  me  (GLM)  vs.  RMS  Lg  (NORSAR)  at  Various  Sites 

Test  Site 

mt,(P,)-RMS  Lg ,  # 

mt(Pb)-RMS  Lg.* 

mg(Pmtji-RMS  Lg ,  # 

BSW 

-0.08410.008  39 

0.11410.009  39 

-0.10310.019  24 

BNE 

-0.19210.012  23 

0.02810.014  23 

-0.13110.032  18 

BTZ 

-0.46610.018  13 

-0.17810.015  13 

0.03310.015  13 

-0.14310.033  9 

Deg 

-0.10610.028  5 

0.07410.029  5 

-0.12210.118  2 

KTS 

-0.40410.006  80 

-0.13210.008  80 

0.07410.008  80 

-0.12010.015  53 

NNZ 

-0.57710.019  15 

-0.34610.022  15 

-0.14310.023  15 

0.12410.051  10 
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Figure  17.  The  spatial  patism  dl  n%-Lf  residuals  of  Semipaialinsk  sspbsions  using  the  mg  (Pm)  of  this  study 
(Table  5)  and  RMS  Lg  values  teported  by  Ringdal  ef  siL  (1002).  The  locations  are  based  on  Bocharov  ef  af. 
(1000)  and  Thuiber,  Quin,  and  Rich«ds  (1003),  and  the  tectonics  are  based  on  Bonham  ef  al  (1080)  and  Le«h 
(1087).  The  pattern  suggests  some  dWerenoe  in  ttte  source  medium  across  the  Chinrau  faul  separating  the 
northeastern  and  aouthwestam  portion  of  the  test  ske.  The  mean  -Lg  bias  between  SW  and  NE  Balapan  is 
about  0.08  m.u. 
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Tn^(Tnb[ML4,Pmax\) 


vs.  rrib\ML4,P7ncLx\ 


D^geUn:  mb[ML4.Pmax]=Lg[N0RSAJi]+0. 07*(0. 029) 
Shagan:  mb[ML4.Pmax]=Lg[NOJiSAR]+0. 074(0.  OOB) 
SR.NE:  mb[ML4.Pmca]  =Lg[N0RSAR]+0. 028(0. 0 1 4) 
SR.  TZ:  mb[ML4.pTnaxhLg[N0RSAR]+0. 033(0. 0 1 5) 
SR.SW:  mh[ML4.Prnax]=Lg[N0RSAR]-¥0. 1 14(0.009) 


Fiflur*  IS.  Comparison  of  irio  (Pmn)  against  RMS  Lg  (NORSAR)  raportad  by  Ringdal  at  a/.  (1992).  Tha  slopa  is 
fixad  at  1 . 
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AVERAGED  SW-NE  BIAS  AT  WORLDWIDE  SEISMIC  STATIONS 
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FIgura  10.  Avaragad  SW-NE  bias  at  aach  WWSSN  station.  Positiva  symbols  raprasant  tha  stations  whara  ampn- 
tuda  of  BSW  avants  is  anhancad  ralativa  to  that  of  BNE  avants  of  tha  sama  souica  strangth.  This  pattam  raflacts 
tha  diffaranoa  of  path  affacts  on  thasa  two  adjaoant  fast  sitas.  For  nalwoik  with  an  unavan  gaographical  distrlMJ- 
tion  of  stations  (such  as  ISC),  tha  simpla  natwork  avaraging  of  station  magnitudas  can  only  aliminala  tha  path  affact 
to  certain  axtanL 


V  J  Maanltud«-fnagnKu(l0  Regression  wtth  Uncertain  Data 


The  standard  approach  of  yield  estimation  is  to  use  known-yield  events  to  construct  a  magnitude- 
yield  relationship  which  is  then  utilized  to  estimate  the  yield  of  other  events.  Typically  either  the  yield  or 
the  is  assumed  to  be  exact  in  the  regression.  The  magnitude-magnitude  regressions  of  Sections  V.1 
and  V.2  are  such  that  the  independent  variable  is  assumed  to  be  exact  arxf  that  the  slope  Is  fixed  at  1 . 
In  reality,  however,  both  the  yields  and  the  magnitudes  are  subject  to  error.  The  regression  result  could 
be  misleadirx)  if  we  simply  assume  that  the  yields  of  19  Semipalalir^  explosions  published  in  Soviet 
literature  are  exact.  It  has  been  speci  oviets  might  have  rounded  8  of  the  armounced  19 

yields  to  the  nearest  5  kt  or  10  kt.  Thus  an  announced  yield  of  100  Id  could  mean  something  actually 
nieasured  between  95  kt  and  104  kt  if  rounding  is  the  only  source  of  error.  It  could  also  indicate  that 
perhaps  100  kt  was  the  designed  energy  release,  arxl  the  actual  yield  was  somewhere  nearby.  UnISte 
the  magnitude-yield  regression,  the  rounding  errors  can  be  controlled  in  the  magnitude-magnitude 
regression  models.  The  starxlard  errors,  however,  stilt  can  not  be  ignored. 


A  more  general  regression  routine  is  given  in  this  section  to  take  the  standard  errors  in  the  magni¬ 
tudes  into  account.  For  each  (X,  Y)  pair  (where  X  and  Y  are  two  different  magnitudes),  we  use  a  ran¬ 
dom  number  generator  to  produce  a  perturbed  (X.  Y)  pair  according  to  their  uncertairAy  distribution.  A 
standard  least-squared  regression  is  then  performed  for  each  data  set  of  perturbed  pseudo¬ 
observations.  The  procedure  is  repeated  for  several  hurrdred  iterations,  and  all  the  resulting  calibration 
curves  are  then  used  to  infer  the  ensemble  behavior.  This  "doubly-weighted  least-squares  scheme” 
[OWLSQ]  is  an  extension  to  the  “ordinary  weighted  least-squares"  [OWLS]  in  which  only  errors  in  the  Ys 
would  be  used  to  adjust  the  inferred  parameters. 

The  “upper  95%  confidence  limit"  of  the  predicted  mi,  at  a  given  X  (say,  Xq)  can  be  computed  as 
follows: 


?(max)  +  t(D.O.F..0.975)Io*(Y)  +  o2(regression)(-^  + 


(Xo-X)^ 

Z(Xj  - 


[19] 


where  N  -  number  of  data  points  used  in  the  regression,  D.O.F.  -  N-2,  a(  Y  )  -  the  mean  S.E.  in  the 
network  m^  [2]  used  in  the  regression,  o(regression)  «  the  o  of  residuals,  tTx,(max)  -  estimate  of  the 
largest  possible  mean  [2]  at  the  given  X,  X  is  the  mean  X  used  in  the  regression,  and  t(D.O.F.. 
0.975)  is  the  97.5  percentile  of  Student's  t  distribution  at  "D.O.F.”  degrees  of  freedom.  Most  statistics 
textbooks  have  a  tabie  of  such  values  after  Fisher  and  Yates  (1963).  Several  comrmnly  quoted 
t(D.O.F.,  0.975)  values  are  listed  in  Table  13. 


Table  13. 97.5  Percentile  of  t  Distribution 


D.O.F. 

5 

10 

20 

30 

40 

60 

oo 

t(D.O.F.) 

2.571 

2.228 

2.086 

2.042 

2.021 

1.960 

The  "lower  95%-confidence  limit”  can  be  computed  in  a  similar  way: 
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t(min)  -  t(D.O.F..0.975)[<j2(Y)  +  <j2(fegression)(^  + 


(Xo-X)^ 

Z(X,-X)2'' 


[20] 


In  the  case  where  each  error  range  in  the  X  and  the  Y  is  reduced  to  zero,  then  all  the  random 
resai.iplings  will  simply  produce  identical  replica  ol  the  original  data  set.  Consequently,  the  several  hurt- 
dred  regressions  will  all  give  the  result  identical  to  a  single  call  of  the  standard  least-squares.  This  Illus¬ 
trates  how  the  “doubly  weighted  least-squares”  [OWLSQ]  would  degenerate  to  the  starxlard  least- 
squares  when  the  urKertairties  in  the  X  and  Y  shrink  to  zero.  By  the  same  reasoning,  this  “doubly 
weighted  least  squares”  regression  method  is  also  an  extension  to  the  conventional  “weighted  least- 
squares”  in  which  only  the  errors  in  the  Y's  would  be  used  to  adjust  the  inferred  parameters. 

Table  14  summarize  the  magnitude-magnitude  regression  models  for  all  test  sites  with  the  algo¬ 
rithm  described  above.  is  used  as  the  independent  variable  in  all  regresskxts.  The  standard 

errors  associated  with  the  event  magnitudes  as  listed  in  Table  5  are  taken  into  account.  Outliers  are 
rejected  and  then  the  regression  was  repeated  unti*  no  outlier  was  detected.  Figures  20  through  24 
show  the  regressions  of  different  phases  on  at  major  test  sites.  Figure  21  show  the 

mbfPmu)  regression  with  the  outliers  retained.  The  uncertainties  in  both  X  and  Y  variables  were 
accounted  for  by  the  same  bootstrap  resampling  procedure. 

For  comparison.  Figures  25  and  27  show  the  standard  least-square  regression,  ignoring  the  starv 
dard  errors  in  and  M$ .  The  regression  results  using  other  phases  are  summarized  in  the 

Appencfx  D. 
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Table  14. 


Uncertain 


Magnitude  inagnitude  Regression  Results 


X  &  Y,  Bootstrap  Used,  CXitlier  Rejected 


N  Slope  Intercept  o(Y) 


t.043±0.021  -0.55810.1 16  0.073 


1.05010.031  -0.57410.171  0.065 


0.75610.167  0.87710.855  0.102 


1.04710.028  -0.60510. 153  0.095 


0.98210.010  -0.09610.057  0.036 


0.96510.024  -0.00810.130  0.043 


1.05110.040  -0.51310.228 


0.99910.049  -0.20910.288  0.043 


0.99010.053  -0.1 4010.322  0.025 


1.05010.017  -0.46910.101  0.034 


0.99110.125  -0.255l0.683  0.113 


BSW 


NNZ 


Sahara 


Pmax 


Pmax 


Pmax 


Pmax 


Pmax 


Pmax 


Pmax 


Pmax 


Pmax 


Pmax 


NTS 

Pmax 

Pa 

192 

1.03610.019 

-0.78610.104 

0.108 

0.96 

PMA 

Pmax 

Pa 

53 

1.00910.031 

-0.61110.173 

0.085 

0.96 

RNA 

Pmax 

Pa 

15 

0.90210.201 

-0.08711.034 

0.185 

0.73 

YFT 

Pmax 

Pa 

126 

1.03110.023 

-0.76610.129 

0.113 

0.96 

KTS 

Pm« 

Pa 

184 

0.95410.009 

-0.20810.050 

0.050 

0.99 

Deg 

Pmax 

Pa 

91 

0.98010.017 

-0.336l0.092 

0.056 

0.98 

BNE 

Pmax 

Pa 

30 

1.01410.040 

0.058 

0.98 

BTZ 

Pmax 

Pa 

19 

0.97710.049 

-0.35910.293 

0.062 

0.97 

BSW 

Pmax 

Pa 

48 

1.00310.054 

-0.49710.328 

0.036 

0.99 

NNZ 

Pmax 

Pa 

29 

1.02610.017 

-0.553l0.099 

1.00 

Sahara 

Pmax 

Pa 

7 

0.73310.108 

0.98010.591 

0.066 

0.98 

NTS 

Pmax 

Pn 

204 

1.05410.029 

-0.46410.153 

0.137 

0.95 

PMA 

Pmax 

Pn 

54 

0.96210.068 

0.15310.371 

0.090 

0.95 

RNA 

Pmax 

Pn 

22 

0.63810.183 

1.31510.930 

0.148 

0.72 

YFT 

Pmax 

Pn 

141 

1.00210.033 

-0.21610.172 

0.125 

0.96 

KTS 

Pmax 

Pn 

100 

0.92510.048 

0.37410.268 

0.243 

0.86 

Deg 

Pmax 

Pn 

52 

0.83410.100 

0.82210.528 

0.263 

0.72 

BNE 

Pmax 

Pn 

20 

0.68510.177 

1.50011 .012 

0.542 

0.35 

BTZ 

Pmax 

Pn 

6 

0.41810.362 

3.49312.141 

0.170 

0.67 

BSW 

Pmax 

Pn 

30 

0.74210.209 

1.47111.261 

0.218 

0.58 

NNZ 

Pn 

20 

0.92110.086 

0.96610.465 

0.248 

0.84 

NTS 

Pmax 

Lg 

168 

0.84910.029 

1.06610.156 

0.151 

0.91 

PMA 

Pmax 

Lg 

47 

0.80510.070 

1.30210.391 

0.118 

0.89 

RNA 

Pmax 

Lg 

22 

0.77610.183 

1.14610.928 

0.153 

0.78 

YFT 

Pmax 

Lg 

107 

0.83410.036 

1.17810.191 

0.151 

0.92 

KTS 

Lg 

74 

0.72510.066 

1.45910.390 

0.113 

0.86 

Table  14.  Magnitude :magnitude  Regression  Results  (Cont) 


Uftcertain  X  &  Y.  Bootstrap  Used.  CXitlier  Rejected 


Site 

X 

Y 

N 

Slope 

Intercept 

o(Y) 

P 

Deg 

Pmax 

Lg 

12 

0.826±0.148 

0.111 

0.91 

BNE 

Pmax 

Lg 

20 

0.647±0.138 

1.92010.806 

0.126 

0.78 

BT2 

Pmax 

Lg 

15 

0.64410.126 

0.107 

0.88 

BSW 

Pmax 

Lg 

27 

0.55710.221 

2.48611.335 

0.107 

0.63 

NNZ 

Pmax 

Lg 

24 

0.71010.051 

1.95010.304 

0.109 

0.97 

NTS 

Pmax 

LR 

134 

1.17310.031 

-2.04210.174 

0.189 

PMA 

Pmax 

LR 

51 

1.11310.049 

0.169 

0.88 

YFT 

Pmax 

LR 

80 

-1.81810.214 

0.170 

0.92 

KTS 

Pmax 

LR 

138 

0.97210.017 

-1.50110.102 

0.191 

0.89 

Deg 

Pmax 

LR 

55 

1.31210.045 

-3.31 010.242 

0.196 

0.89 

BNE 

Pmax 

LR 

26 

1.21210.069 

-2.96210.399 

0.209 

0.86 

BTZ 

Pmax 

LR 

18 

-3.20210.595 

0.238 

0.83 

SSW 

Pmax 

LR 

42 

0.47910.033 

1.50610.197 

0.126 

0.66 

NNZ 

Pmax 

LR 

29 

0.99610.023 

-1.23810.140 

0.189  1  0.95 

At  M5  -  4.0,  the  corresponding  mt,{P^  at  NTS,  KTS,  and  NNZ  are  5.151,  5.659,  and  5.259, 
respectively.  Once  again,  under  the  assumption  that  Ms  is  relatively  less  susceptible  to  source  medium 
and  propagation  effect  than  other  phases,  then  Us  -mt,  bias  would  reflect  primarily  the  mi,  bias,  and  our 
results  would  support  the  claim  that  excitation  at  NNZ  is  intermediate  between  that  of  NTS  and 
KTS,  with  nit,  the  largest,  given  the  same  Us  level.  This  rt  suil  is  In  agreement  with  that  based  on 
fixed-slope  regression  as  did  in  the  previous  section. 

The  :mf,{Lg)  regression  models  appear  to  lead  some  inconsistency  with  the 

n^biPmwd  'Ms  models.  At  mt,[Lg)  of  5.5,  the  corresponding  nibiPmaxi  at  NTS,  KTS,  and  NNZ  are 
5.223,  5.574,  and  5.000,  respectively.  This  is  in  conflict  with  the  excitation  relative  to  Ms .  We  have 
noticed  this  inconsistency  in  the  previous  section.  It  would  seem  that  our  mi,(Lg)  values  are  too  high  for 
NNZ  explosions.  To  remove  this  ifKX>nsistency,  at  least  two  simple  approaches  can  be  adopted.  The 
first  one  is  to  revise  Nuttli's  path  attenuation  parameters  for  NNZ  explosions.  Perhaps  his  O  values 
based  on  the  coda-Q  method  are  too  tow  which  caused  the  over-compensation  of  the  attenuation  effect. 
An  10%-15%  increase  in  the  Q  values  would  bring  the  mi,{Lg)  values  down  by  0.2-0.3  m.u.,  which 
would  be  sufficient  to  remove  the  inconsisterxry.  Alternatively,  assuming:)  the  Q  values  derived  by  the 
coda-Q  method  do  match  perfectly  with  those  derived  with  the  time-domain  amplitude  decay  relation¬ 
ship,  then  probably  the  apparently  strongest  m^(Lp)  excitation  at  NNZ  is  due  to  erroneous  scaling  con¬ 
stant  in  the  mi,(Lg)  formula.  In  this  study,  we  have  used  a  single  formula  to  compute  the  station 
m^{Lg) ,  Equation  (1],  which  is  the  same  formula  that  Nuttli  used  in  his  Lg  study  of  NTS,  KTS,  and  NNZ 
explosions: 

.  4.0272 *  loom  4  ilooWkm))  4  ^l°a|sln(,„  . 
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This  formula  is  scaled  so  that  a  seismic  source  with  l-sec  Lg  amplilude  of  110  pm  at  10  tun  epi- 
central  distance  would  correspond  to  a  fng(Lg)  of  4.0272  2.0414  *  0.3333  •  1.4019  *  0.0000  •  5.000, 
which  was  suggested  to  be  appropriate  for  eastern  North  America  and  Semipalatinsk.  However,  NuttK 
(1980)  also  suggested  that  a  seismic  source  in  Iran  with  the  ISC  bulletin  5.0  should  excite  Lg  ampl- 
tude  of  270  microrrs  at  a  10-km  extrapolated  distance.  Possible  explanations  for  INs  apparera  bias  of 
0.26-0.39  m.u.  inciude  [1]  uncertainty  in  the  ISC  buletin  values  (Nuttli,  1981),  [2]  differences  in  Lg 
excitation  reiative  to  ntb .  or  [3]  differences  in  the  upper  mantle  property  which  could  cause  a  bias  In  rr^ . 
Many  Lg  studies  have  used  the  ISC  for  “normaizing”  their  Lg  magnitude  scale  (a.g. ,  page  2146  of 
Nuttli,  1986a;  page  128  of  Israeison,  1992).  Without  careful  re-processing,  the  ISC  bulletin  values 
typically  would  be  associated  with  iarge  uncertainty.  The  used  in  this  study  are  those  that  have  been 
correrded  for  the  station  amplificatiorts  as  weU  as  the  path  effects,  and  therefore  we  can  rule  out  the 
possibility  of  targe  uncertairty  in  mu .  It  would  seem  plausible  that  NNZ  might  have  a  Lg  scale  similar  to 
that  in  Iran: 

.  3.6372  *  fogA(A)  .  ilog(A<km))  »  7  . 

With  a  reduction  of  0.39  m.u.  in  the  /T^(Lp)  scaling  constant,  it  is  sufficient  to  make  the  ntg  :Lg 
and  mg  :Ms  modeis  consistent  across  the  three  test  sites  we  discussed.  This  adjustment  would,  how¬ 
ever,  raise  a  doubt  whether  the  Lg  scaiing  appropriate  lor  eastern  U.S.  and  Semipalatinsk  is  really  tran¬ 
sportable  to  Novaya  Zemiya  as  is  Ms . 
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Figura  23.  RagrcMing  n^{Lg)  on  n%,(P^)  with  outlior*  oxciudod. 
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Figur*  24.  Ragratting  Ms  on  with  outlior*  oxdudod. 
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VI.  CONCLUSION  ANC  RECOMMENDATIONS 


The  most  important  goal  of  this  contract  was  to  provide  AFTAC/TT  a  huge  explosion  database  of 
quality  measurements  with  which  AFTAC/TTR  can  conduct  various  follow  on  research  and  analysis. 
This  goal  has  been  accomplished  under  Task  1  (cl.  Baumstark  and  Wagner,  1994).  Furthermore,  the 
delivered  product  also  includes  all  the  digital  waveforms,  fully  documented  and  converted  to  CSS  3.0 
format,  which  were  used  in  making  the  measurements.  The  measurements  made  by  Teledyne  analysts 
can  be  verified  or  revised  by  bringing  up  the  waveforms  with  an  interactive  display  routine  “geotoor 
developed  by  Teledyne  and  SAIC  (Henson  and  Coyne,  1993).  The  detailed  description  of  the  measured 
parameters  as  well  as  the  procedure  used  in  making  those  measurements  can  be  found  in  an  accom¬ 
pany  final  report  (Baumstark  and  Wagner,  1994). 

To  fully  digest  the  huge  explosion  dataset  would  require  a  significarl  amount  of  time  and  labor 
which  are  certainly  beyond  the  level  of  efforts  allocated  to  this  Task  3.  The  results  presented  in  this 
report  should  be  regarded  as  preliminary  in  nature  with  many  topics  remaining  to  be  explored.  That  is, 
this  part  of  the  final  report  is  a  pilot  study  which  only  serves  to  illustrate  how  these  carefully  measured 
parameters  can  be  used  in  some  magnitude^ield  analysis  with  the  statistical  software  package 
developed  at  Teledyne.  In  Sections  II  and  III.  we  presented  the  methodlogies  used  in  determining  the 
station  magnitudes  arxl  the  event  magnitudes,  respectively.  The  GLM  algorithm  not  only  produces  the 
network-averaged  magnitudes,  which  are  optimal  in  the  maximum-likelihood  sense,  it  also  generates  two 
types  of  by-products:  the  station  corrections  and  the  path  corrections.  For  explosions  clustered  in  a 
region  of  the  “same'*  geological  /  geophysical  characteristics,  the  residuals  (after  the  mean  event  magni¬ 
tude  and  the  station  effect  have  been  removed)  typically  exhibit  some  coherent  pattern  reflecting  the 
propagation  effect  (cf.  Jih  and  Wagner,  1991,  for  detailed  discussion).  The  GLM  algorithm  has  been 
upgraded  during  the  past  several  years  to  include  the  path  tenns  in  the  joint  inversioa  This  GLM  routine 
has  been  fully  tested  and  it  has  been  applied  to  all  time-domain  magnitudes  used  in  this  study,  although 
much  attention  was  paid  to  for  illustrative  purpose  as  explained  before.  The  explosion  data  set 
established  under  this  contract  has  a  very  broad  spatial  coverage.  As  a  result,  the  station  correc¬ 
tions  derived  by  GLM  routine  appear  to  match  with  the  (receivers’)  tectonics  fairly  well.  These  station 
corrections  have  been  used  in  computing  the  event  magnitudes  of  those  isolated  explosions  that  were 
not  included  in  the  GLM  inversion.  Thus  the  event  magnitudes  of  PNEs  and  explosions  in  Amchitka. 
South  Novaya  Zemiya,  Colorado,  and  Mississippi  are  “sub-optimal”  in  the  sense  that  only  receiver 
effects  were  removed,  but  not  the  path  effects.  Another  non-standard  algorithm  we  used  in  this  study  is 
a  magnitude-magnSude  regression  routine  which  accounts  for  the  uncertainties  In  the  event  magnitudes 
through  Monte-Carlo  resampling.  A  feature  of  automatic  outlier  rejectton  has  been  added  to  this  routine. 
This  feature  should  be  used  carefully,  since  “outliers"  of  the  magnitude-magnSude  regression  (or 
magnitude-yield  regressnn)  may  sometimes  be  diagnostic  of  the  peculiar  testing  practice  or  the 
emplacement  media. 

The  mu  results  presented  in  this  report  are  very  consistent  with  those  of  our  previous  study  (cf. 
Jih  et  a/..  1993)  using  essentially  the  same  technologies  and  different  explosion  dataset.  In  comparing 
mt,.  Ms,  and  mi,(Lg) .  we  observe  some  Inconsistency  in  the  scaling  relationships.  The  systematic  bias 
in  the  scale  can  be  renK>ved  either  by  raising  Nuttli's  path  attenuation  coefficients  (which  were 
derived  with  the  coda  O  method),  or  by  revising  the  absolute  Lg  excitation  level  for  Novaya  Zemiya 
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explosions.  Whichever  approach  adopted,  there  would  be  some  question  about  whether  the  Lp  scale  is 
really  as  transportable  as  Nuttli  suggested.  The  good  results  of  applying  Nuttli's  NTS  Lp  yield  formula  to 
Semipalatinsk  explosions  may  turn  out  to  be  a  special  case  rather  than  a  general  rule,  if  the  Lp  scale 
suitable  for  Iranian  Plateau  were  used  in  computing  the  /nj,(Lp)  of  Novaya  Zemlya  explosions,  then  the 
magnitude-yield  relationships  based  on  mb.  Ms  >  and  Lg  would  be  consistent. 

Under  this  contract,  the  Moscow-Prague  formula  has  been  used  in  computing  Ms  when  A  is  larger 
than  25  degrees.  Recently  Herak  and  Herak  (1993)  suggest  that  the  amplitude  of  20-second  Rayleigh 
wave  should  decrease  as  A~^^,  instead  of  A~'^  as  in  the  Moscow-Prague  formula.  Herak's  new  for¬ 
mula  is  very  similar  to  that  of  von  Seggem  (1977).  It  would  also  be  more  consistent  with  the  formula 
Nuttli  and  Kim  (1975)  suggested  for  10*  <  A  <  25*. 
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APPENDIX  A:  PREREQUISITE  MATHEMATICS  FOR  MAXIMUMUKELIHOOO  ESTIMATOR 

Proposition  1.  Let  X  be  a  Gaussian  random  variable  wth  the  probability  density  function  [p.d.f.]  g  and 
the  cumulative  distribution  function  [c.d.f]  G,  respectively.  Its  mean  and  variance  are  denoted  by  p  and 
o^,  respectively.  Then 

a 

Jxg(x)dx  =  pG(a)-o2g(a). 


Proof. 


=  pG(a)  -  a^(a). 


In  particular,  when  a  =  «o,  this  integral  gives  the  mean  of  X,  namely,  p. 

Proposition  2.  Let  X  be  a  Gaussian  random  variable  with  mean  p  and  variance  o^,  then  E  [  X  I  X  <  a  ] 
»  p  -  o^g(a)/G(a) . 


Proof.  Let  Y  be  the  random  variable  X  I  x  <  •>  then 

P(Y<b)  =  P(X<blX<a)  =  ^  ^ 

which  is  1  if  b  >  a,  and  G(b)/G(a)  if  b  ^  a.  Therefore,  the  p.d.f.  of  Y  is: 

h(x)  =  0  if  X  >  a  ,  h(x)  =  g(x)/G(a)  if  x  <  a, 
and  the  expectation  of  Y  is 

••  a 

E(Y)  =  J  xh(x)dx  =  J  xg(x)/G(a)dx 

_  ^^y  proposition  1) 

G(a) 

-  p  -  <j2g{a)/G(a) . 

Similarly  it  can  be  shown  that  EIXlX>a]»p-f  o^(a)/G(-a) .  Note  that  the  conditional  expec¬ 
tation  E  [  X  I  X  >  a  ]  is  the  “best”  guess  of  X  under  the  constraint  that  one  knows  only  that  X  >  a. 

In  computing  E[XlX>a],  itis  generally  more  convenient  to  transform  the  random  variable  X 
into  the  standard  random  variable,  Z  ~  N(0,  1),  for  which  the  p.d.f.  and  c.d.f.  are  typically  available  as 
system-furnished  functions  or  pail  of  some  utility  Ibraries  in  the  public  domain.  Let  O  and  ^  be  the  c.d.f. 

and  p.d.f.  of  Z,  respectively,  then  G(a)  » and  g(a)  -  a  ♦(^^].  Therefore. 
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Eixix<al  =  u-o#(-^l^^l. 

O  0 
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APPENDIX  B.  LINEAR  REGRESSION  WITH  CENSORED  OBSERVATIONS 

Consider  the  situation  where  the  independent  variable  can  be  precisely  measured,  and  that  we 
want  to  regress  the  dependent  variable  as  a  linear  function  of  the  independent  variable.  This  is  an 
extension  to  the  smgie-eveni  network  magnitude  determination  discussed  in  Sedion  1.3. 

Suppose  we  have  a  linear  model  of  the  comrrKm  form: 

Y  =  a  +  pX  +  V  . 

where  X  is  the  irxJependent  variable  which  has  a  precision  relatively  much  better  than  that  of  Y,  the 
dependent  variable,  a  and  p  are  the  intercept  and  slope,  respectively,  to  be  determined,  and  v  is  an 
error  term,  v  is  assumed  to  be  a  Gaussian  rarxlom  variable  with  mean  zero  and  standard  deviation  o. 
Furthermore,  assume  that  there  are  four  types  of  data  available: 

[0]  the  observed  measurement,  Y,  is  known  as  yo, 

[1]  Y  is  orriy  known  to  be  less  than  certain  level. 

[2]  Y  is  only  known  to  be  larger  than  certain  level,  and 

[3]  Y  is  only  known  to  lie  between  two  bounds. 

Type  3  data  are  not  uncommon.  The  maiority  of  Soviet  yields  recently  published  by  Bocharov 
et  a/.  (1989)  and  Vergirto  (1989)  actually  fail  in  this  category. 

Elegant  maximum-likelihood  theory  can  be  derived  for  this  model.  Suppose  there  are  no.  n^,  nj, 
and  na  measurements  for  each  type,  respectively.  The  conditional  likelihood  function  of  the  censored 
obsenrations  (  yo,  ti,  ta,  ta)  given  the  intercept  a.  slope  p.  and  o  is 

no  ft, 

L  (  Yo.  ti.  t2.  ta  I  a,  p.  o  )  =nP(  Yj  =  y^  I  a,  p,  o )  •  nP{  Yj  <  t,,  I  a,  p.  a )  • 

i-i  j.i 


nP(  Yj  >  t2j  I  a.  p,  o )  •  nP(  1^  <  Yj  <  ttj  I  a,  p,  a )  . 

j-i  j-i 

and  the  log-likelihood  function  is 

In  L  (  Yo.  ti,  I2.  la  I  a.  p,  a  )  =  -  -^  ln(2jio2)  -  X(yoi  -  a  -  P  Xoj)^  + 

«  2<r  j « 1 

Ri  Rj  03 

£ln  <b(zij)  -H  £ln  <IK-Z2i)  +  Z'n  [‘*KZbi)-«K2iij)l  . 

j-i  i-i  j.i 

where  Zj  ■  [yj  -  a  -  ptj/a;  Yo,  ti,  t2,  and  ta  are  the  vectors  of  the  four  data  types. 

Solving  >  0  implies  immediately  that  the  ft,  the  optimal  estimate  of  o,  must  satisfy  the  follow- 
do 

ing  necessary  condition: 
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^  yo,  -  a  -  p  Xo(  )* 

- tl - 

^  ,  54(2bj)Zbi-^(2.j)z.i 

dInL  dInL 

Solving  m  0  implies  that  the  sum  of  the  “rermed  residuals”  should  be  zero.  Solving  ■  0 
oo.  op 

implies  that  the  vector  of  refined  residuals  should  be  orthogonal  to  the  vectors  of  means,  it  follovM  that 

the  optimal  estimate  of  a  and  p  can  be  obtained  by  the  “standard  least  squares”  inversion  with  the  oerv 

sored  data  all  replaced  by  their  conditional  expectations,  Le. ,  the  “refined  observations.”  Thus  a  can  be 


solved  iteratively  with  [6]  along  with  a  and  p  using  the  EM  algorithm.  In  the  noivcensored  case,  this 
“doml02‘’ oo6e  gives  results  Ueniical  to  those  derived  by  the  standard  least  squares. 


I 
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APPENDIX  C:  EMILS  EVENT  MAGNITUDES 


Table  C.  EMILS  Event  Magnttudes' 


Event 


US651029 


US691002 


US711106 


NZ730927 


NZ731027 


NZ741102 


NZ751018 


Site 


Alaska 


Alaska 


Alaska 


5.9Q±0. 


6.55±0.06 


6.8e±0. 


5.54±0.06 


6.2710.06 


6.6110.07 


US631026 

Other 

5.3210.21 

5.1310.45 

US641022 

Other 

4.6210.00 

4.3410.00 

US671210 

Other 

5.0910.07 

US690910 

Other 

4.7810.00 

4.1310.00 

US730517 

Other 

5.1310.14 

4.9210.15 

5.2010.07 


6.0410. 


6.4010. 


5.2910.53 


4.3010. 


4.0210.41 


4.4710. 


4.2410.53 


4.0610.26 


4.8310. 


5.8010.13 


5.5710. 


6.1710. 


4.3010. 


PN660930 

PNE 

PN671006 

PNE 

PN680521 

PNE 

PN690902 

PNE 

PN690908 

PNE 

PN690g26 

PNE 

PN691206 

PNE 

PN700625 

PNE 

PN701212 

PNE 

PN701223 

PNE 

PN710323 

PNE 

PN710702 

PNE 

PN710710 

PNE 

PN710919 

PNE 

PN711004 

PNE 

PN711022 

PNE 

PN720411 

PNE 

PN720709 

PNE 

PN720820 

PNE 

PN720904 

PNE 

PN720921 

PNE 

PN721003 

PNE 

PN721124a 

PNE 

PN721124b 

PNE 

PN730815 

PNE 

PN730828 

PNE 

5.8310.06 


7.1410.10 


6.8110.10 


.08 


5.0110.09 


4.7410.11 


5.3610.07 


5.6110.06 


6.9710.12 


6.6410.12 


6.5310.08 


4.2610.16 


5.1610.09 


4.4910.10 


4.9410.19 


5.6210.09 


4.8710.17 


6.0910.1 1 


6.2610.10 


5.7210.10 


4.5110.34 


5.1110.12 


5.3110. 


4.5510.10 


5.3710.11 


4.7910.1 1 


4.4810.15 


5.6710. 


5.1010.00 


4.5810.18 


5.7810.09 


5.94+0.10 


5.5410.08 


4.2710.30 


4.7510.12 


4.2010.13 


4.2210.11 


4.9910.09 


5.4110.07 


6.7710.14 


6.5510.14 


6.2310.10 


4.7210.12 


3.7910.25 


5.0410.10 


4.1410.11 


4.0410.03 


4.9710.14 


4.6110.21 


5.4910.08 


5.5610.10 


5.3110. 


4.2210.42 


4.3210.14 


4.0410.11 


4.3110.10 


4.9210.12 


4.4310.12 


4.3510.30 


5.7710.17 


5.8210. 


6.8710. 


6.9610. 


6.7410. 


4.3510.12 


5.7010.19 


6.0010. 


5.6010. 


5.4410. 


3.3910.34 


a 


4.8910.16 


5.6010.26 


5.9810. 


4.5410.49 


4.9510.15 


4.1710.10 


5.4411.15 


5.2210.01 


4.4610. 


4.5310.07 


08 

5.3010.07 

5.0510. 

08 

5.5010.00 

r  5.8510.10 

3.5610.12 

4.9910.00 


4.4410.10 


5.2410. 


4.7410.13 

5.1110. 

4.7810.10 

5.3010. 

08 

5.0610.06 

5.0910. 

09 

4.8510.07 

4.8010.00 


4.4910. 


4.8810.13 


4.5210.17 


4.6210.13 


4.7410.11 


4.5310.23 


4.8410. 


4.9410. 


5.9210. 


4.77+0.00 


4.4310.00 


5.2610. 


4.6710.15 


4.4410. 


1)  from  Baumstarfc  and  Wagner  (1994). 


Event 


PN730919 


PN730930 


PN731026 


PN740708 


PN740814 


PN740829 


PN741002 


PN750812 


PN750929 


PN761105 


PN770726 


PN770810 


PN770820 


PN770910 


PN780809 


PN780824 


PN780921 


PN781007 


PN781017 


PN790812 


PN790906 


PN791004 


PN791007 


PN801008 


PN801101 


PN801210 


PN610525 


PN810902 


PN81 0926a 


PN81 0926b 


PN811022 


PN820730 


PN820904 


PN820925 


PN821010 


PN821016a 


PN821016b 


PN821016C 


PN8210l6d 


PN830710a 


PN830710b 


Site 


PNE 


PNE 


PNE 


PNE 


PNE 


PNE 


PNE 


PNE 


PNE 


PNE 


PNE 


PNE 


PNE 


PNE 


PNE 


PNE 


PNE 


PNE 


PNE 


PNE 


PNE 


PNE 


PNE 


PNE 


PNE 


PNE 


PNE 


PNE 


PNE 


PNE 


PNE 


PNE 


PNE 


PNE 


PNE 


PNE 


PNE 


PNE 


PNE 


PNE 


PNE 


Table  C.  EMILS  Event  Magnitudes' 


Pt 


(Cent) 


4.9710.11 


5.3110.08 


4.8110.18 


4.6610.16 


5.5210.06 


4.9910. 


4.5910.21 


4.9810.11 


5.0710.11 


5.1210.06 


5.0810.06 


4.9410.09 


5.0510.07 


.06 


5.6810.07 


5.2810.11 


5.4110.15 


5.0610.11 


5.7610. 


5.0010.09 


5.0710.11 


5.5710.11 


5.0110. 


5.1810.12 


5.1110. 


4.6710.16 


5.4210.11 


5.2210.17 


5.2610.10 


5.1310.13 


4.8710.11 


5.1610.10 


5.4010.15 


5.1210.08 


5.0410.10 


5.0810. 


5.1610.10 


5.4210.07 


5.3210.09 


5.3410.07 


4.6910.10 


4.9910.09 


4.5110.22 


3.9110.44 


5.1610.06 


4.6910.08 


3.8110.18 


4.6210.12 


4.6610.13 


4.8310.10 


4.8510.06 


4.6910.06 


4.7910. 


4.6610.10 


5.3210.07 


4.8910.10 


4.9810. 


4.8110.13 


5.3810.08 


4.6810.08 


4.8010.12 


5.3310.10 


4.7510.10 


4.8310.14 


4.9110.10 


4.3310.21 


5.1610.10 


4.7710.11 


4.8510.10 


4.6410.10 


4.7210.10 


4.7710.10 


4.9810.13 


4.8510. 


4.6810.11 


4.7010.13 


4.7610.10 


4.9410.10 


5.0610.07 


5.0710.07 


4.4710.10 


4.7010.13 


4.3810.21 


4.1110.37 


4.8610.09 


4.5810.14 


3.9110.03 


4.1810.29 


4.4510.11 


4.4210.11 


4.5810.07 


4.4310.12 


4.4110.13 


4.9610. 


4.5310.15 


4.7110.08 


4.4810.16 


4.9510.11 


4.0310.30 


4.3910.12 


4.9710.15 


4.4210.11 


4.5910.14 


4.6410.12 


4.4610.22 


4.8910.12 


4.5010.12 


4.6710.11 


4.3710.13 


4.4610.10 


4.4810.11 


4.5810.13 


4.5710.10 


4.5710. 


4.3710.12 


4.6510. 


4.7210. 


4.7210.07 


5.6310. 


5.7110.00 


5.1910.00 


4.7510.00 


5.3610.00 


5.1810.31 


4.9510.48 


5.1510.04 


4.6710.00 


4.3310.00 


5.0410.41 


4.9210.67 


2.5910.00 


3.5410.31 


5.1310.00 

5.4210.23 

4.8010.00 

4.0610.00 

5.0210.00 


5.6110.00 


5.9710. 


5.2610. 


5.1810. 


5.0110.00 


5.2610.10 


5.1410.04 


3.3310.35 


3.6610. 


3.1810.13 


3.6110. 


2.9610. 


3.4010.17 


3.4910.13 


3.4010.12 


3.2310. 


3.2910.13 


3.5810.11 


3.2010. 


3.1410.15 


3.5110.11 


4.0010.07 


3.6210. 


3.4410. 


3.6610. 


2.9110.27 


3.2010.12 


3.7510. 


3.0310.10 


2.9310.11 


3.0010. 


3.3610.07 


3.1210.07 


3.1410. 


Table  C.  EMILS  Event  Magnitudes'  (Cent) 


PN830710C 

PNE 

5.2110.07 

4.8210.07 

4.5010.12 

PN830924a 

PNE 

5.2210.12 

5.0010.13 

4.6510.10 

PN830924b 

PNE 

4.9710.20 

4.7110.18 

4.3010.11 

PN830924C 

PNE 

4.6610.14 

1  4  4110.18 

PN830924d 

PNE 

5.1610.20 

4.9510.14 

4.5610.18 

PN830924e 

PNE 

5.1810.22 

4.8810.20 

4.6210.18 

PN830924f 

PNE 

5.2710.15 

4.9710.12 

4.4610.14 

PN840721a 

PNE 

5.2410.09 

5.0110.09 

4.6810.10 

PN840721b 

PNE 

5.2210.07 

5.0110.07 

4.6410.13 

PN840721C 

PNE 

5.2710.07 

4.9410.08 

4.6310.10 

PN840811 

PNE 

5.0810.14 

4.9110.08 

4.3610.12 

PN840825 

PNE 

5.4010.08 

5.1310.09 

4.8010.08 

PN840827 

PNE 

4.6910.09 

4.7210.12 

4.3610.21 

PN840828a 

PNE 

4.3410.23 

4.3210.26 

3.7510.36 

PN840828b 

PNE 

4.3110.18 

4.3010.22 

4.0610.14 

PN840917 

PNE 

5.3710.15 

5.0910.11 

4.6910.12 

PN841027a 

PNE 

4.9010.16 

4.6310.13 

3.9510.37 

PN841027b 

PNE 

4.9010.18 

3.8110.50 

PN850718 

PNE 

4.9810.12 

4.6710.13  1 

4.3310.13 

PN870419a 

PNE 

4.5410.30 

3.7910.36 

4.0410.10 

PN870419b 

PNE 

4.4510.22 

4.2910.07 

3.7610.20 

PN870706 

PNE 

5.1010.25 

4.8910.22 

4.8610.26 

PN870724 

PNE 

5.1710.11 

5.0010.09 

4.6610.16 

PN870812 

PNE 

4.9610.15 

4.6910.13 

4.5610.16 

PN871003 

PNE 

5.1310.13 

4.8510.16 

4.5910.14 

PN880822 

PNE 

5.3010.14 

4.9310.14 

4.5610.17 

PN680906 

PNE 

4.9610.23 

4.5410.25 

4.3&t0.21 

5.39±0. 


4.98±0.14 


6.29±0. 


5.33±0. 


5.39±0. 


4.30±0. 


5.02±0. 


4.76±0. 


5.17±0.00 


5.1710.00 


5.0010.00 


5.3010.49 


4.7710.16 


4.8010.39 


5.2510. 


4.6710.07 


Ms 


3.2610.09 


3.9910.12 


3.8210.11 


2.3110.15 


2.4210.16 


2.6310.20 


2.4010.18 


2.3210.16 


3.0710. 


3.3010.06 


1)  from  BaumstBffc  and  Wagner  (1M4). 
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APPENDIX  D:  MAGNITUDE-MAGNfTUDE  REGRESSIONS  WITH  ALTERNATIVE  SCHEMES 


NTS 

Pmax 

Pb 

216 

PMA 

Pmax 

Pb 

58 

RNA 

Pmax 

Pb 

18 

YFT 

Pmax 

Pb 

142 

KTS 

Pmax 

Pb 

195 

Oeg 

Pmax 

Pb 

98 

BNE 

Pmax 

Pb 

30 

BT2 

Pmax 

Pb 

19 

BSW 

Pmax 

Pb 

48 

NNZ 

Pmax 

Pb 

30 

Sahara 

Pmax 

Pb 

7 

NTS 

Pmax 

PMA 

Pmax 

RNA 

Pmax 

YFT 

PfTtax 

KTS 

Pmax 

Deg 

Pmax 

BNE 

Pmax 

BT2 

Pmax 

BSW 

Pmax 

NNZ 

Pmax 

Sahara 

Pmax 

NTS 

Pmax 

PMA 

Pmax 

RNA 

Pmax 

YFT 

Pmax 

KTS 

Pmax 

Oeg 

Pmax 

BNE 

Pmax 

BTZ 

Pmax 

BSW 

Pmax 

NNZ 

Pmax 

NTS 

Pmax 

PMA 

Pmax 

RNA 

Pmax 

YFT 

Pmax 

KTS 

Pmax 

Table  D1.  Magnitude  :magnitude  Regression  Results 


Uncertain  X  &  Y,  Bootstrap  Used.  Outlier  Retained 


Slope  Intercept 


1.061±0. 


1.015±0.033 


0.756±0.167 


1.052±0.026 


0.98710.011 


0.96110.024 


1.05110.040 


0.99910.049 


0.99010.053 


1.05210.018 


0.99110.125 


0.92110.021 


-0.67510.121 


-0.37510.184 


0.87710.855 


-0.64710.141 


-0.12910.061 


0.00610.129 


-0.51310.228 


-0.20910.288 


-0.14010.322 


-0.48310.101 


-0.255l0.683 


-0.16610.116 


-0.53810.197 


-0.08711.034 


0.00110.144 


0.044 


0.048 


0.047 


0.043 


0.025 


0.042 


3 


66 


26 


0.185 


0.177 


0.066 


0.078 


0.058 


62 


0.036 


0.066 


66 


0.90210.201 


0.88710.027 


0.96510.010 


0.99210. 


1.01410.040 


0.97710.049 


1.00310.054 


1.02910.017 


0.73310.106 


0.96110.033 


0.96210.06B 


0.63810.183 


0.93710.047 


0.79710.102 


0.68510.177 


0.41810.362 


0.70210.202 


0.92110.086 


-0.40510.116 


-0.55910.226 


-0.35910.293 


-0.497l0.328 


-0.57410. 


-0.009l0.173 


0.15310.371 


1.31510.930 


0.62410.209 


0.25710.265 


0.97810.538 


1.50011.012 


3.49312.141 


1.78211.217 


0.96610.465 


0.090 


0.148 


0.259 


0.384 


0.315 


0.542 


.170 


0.359 


0.248 


0.77610.031 

1.41010.168 

0.83110.072 

1.15010.402 

0.77610.183 

1.14610.928 

0.69310.038 

1.88010.206 

0.72510.066 

1.45910.390 

.135 

0.87 

.153 

0.78 

.259 

0.77 

.113 

086 

i.  Magnitude rmagnitudd  Regression  Results  (Cont) 


Uncertain  X  &  Y,  Bootstrap  Used,  Outlier  Retained 


N  I  Slope  Intercept 


0.826±0.148 


0.647±0.138 


0.64410.126 


0.55710.221 


0.71010.051 


1.14210.035 


1.06410.047 


o(Y) 


0.111 


0.126 


0.107 


0.107 


0.109 


0.220 


92 


0.203 


0.269 


Deg 

Pmax 

Lg 

BNE 

Pmax 

i-g 

BTZ 

Pmax 

Lg 

BSW 

Pmax 

L-g 

NNZ 

Pmax 

Lg 

NTS 

Pmax 

LR 

PMA 

Pmax 

LR 

YFT 

Pmax 

LR 

KTS 

Pmax 

LR 

Deg 

Pmax 

LR 

BNE 

Pmax 

LR 

BTZ 

Pmax 

LR 

BSW 

Pmax 

LR 

NNZ 

Pmax 

LR 

.82510.846 


1.92010.806 


1.95110.748 


2.48611.335 


1.95010.304 


-1.86210.195 


-1.31710.260 


-1.59810.249 


-1 .23510.098 


-2.47610.233 


-2.96210.399 


-3.20210.595 


-0.81210.345 


-^238l0.140 


0.3,<i9l0.017 


•;  7i»-L3.043 


1.21210.069 


1.25210.099 


0.85710.057 


0.99610.023 


0.238 


0.200 


0.189 


Table  D2.  Magnitude 


Standard  Least 


N 


imagnitude  Regression  Results 


Squares,  CXitlier  Rejected 


NTS 

Pmax 

Pb 

PMA 

Pmax 

Pb 

RNA 

Pmax 

Pb 

YFT 

Pmax 

Pb 

KTS 

Pmax 

Pb 

Deg 

Pmax 

Pb 

BNE 

Pmax 

Pb 

BTZ 

Pmax 

Pb 

BSW 

Pmax 

Pb 

NNZ 

Pmax 

Pb 

Pmax 


Pmax 


Pmax 


Pmax 


Pmax 


Prrax 


Pmax 


Pmax 


Pmax 


Pmax 


Pmax 


Pmax 


Pmax 


Pmax 


Pmax 


PiTiax 


Pmax 


Pmax 


Pmax 


Pmax 


Pmax 


Pmax 


Slope 


1.098±0.011 


1.100±0.025 


1.066ct0.131 


1.082±0.015 


0.987±0. 


0.987±0.013 


1.047±0.021 


1.01 110.038 


1.00810.014 


1.05210.011 


1.00110.113 


1.03310.017 


0.99010.038 


1.09710.281 


1.07710.023 


0.96410.006 


1.01910.018 


1.04710.036 


0.98810.055 


1.04310.019 


1.02810.019 


0.74010.066 


1.08510. 


0.98210.043 


1.02610.161 


1.01810.025 


0.93110.047 


0.81710.094 


0.68910.431 


0.40910.229 


0.50110.178 


0.89210.116 


0.85910.025 


0.79110. 


0.97110.176 


0.84010.028 


0.73410.051 


Irrtercept 


-0.857l0.057 


-0.84010. 136 


-0.69110. 


-0.78410. 


-0.1 2510.033 


-0.11810. 


-0.48210.122 


-0.27410.225 


-0.245l0.087 


-0.48110.064 


-0.30010.595 


-0.76310.090 


-0.49910.211 


-1.08511.440 


-1.02110.121 


a(Y) 


0.047 


0.052 


0.085 


0.057 


0.032 


0.034 


0.030 


0.043 


0.025 


0.034 


0.113 


0.071 


0.074 


0.182 


NTS 


PMA 


RNA 


YFT 


KTS 


Deg 


BNE 


BT2 


BSW 


NNZ 


Sahara 


NTS 


PMA 


-0.53810.097 


-0.75410.206 


-0.422l0.327 


0.94610.345 


-0.62510.120 


0.03910.235 


-0.67010.818 


-0.30310.130 


0.28710.265 


0.88510.498 


1.46712.458 


3.54211.340 


2.91111.077 


1.08810.643 


1.02310.133 


1.37010.328 


0.15410.891 


1.15910.150 


1.39810.300 


0.042 


0.047 


0.047 


0.062 


0.028 


0.057 


0.065 


0.123 


NTS 

Pmax 

PMA 

Pmax 

RNA 

Pmax 

YFT 

Pmax 

KTS 

Pmax 

0.121 


.117 


0.192 


0.184 


0.542 


0.169 


0.176 


0.202 


0.118 


0.148 


0.114 


0.113 


Table  D2.  Magnitud6:magnttude  Regression  Results  (Cont) 


Standard  Least  Squares,  Outlier  Rejected 


N  Slope  I  Intercept 


Deg 

Pmax 

Lg 

12 

0.83610.121 

0.76210.688 

0.111 

0.91 

BNE 

Pmax 

Lg 

20 

0.65710.125 

1.85610.724 

0.126 

0.78 

BT2 

Pmax 

Lg 

15 

1.86510.583 

0.107 

0.88 

BSW 

Pmax 

Lg 

27 

0.58410.146 

2.31210.881 

0.106 

0.63 

NNZ 

Pmax 

Lg 

24 

0.71110.039 

Mill  1  I’M 

0.109 

0.97 

Pmax 


Pmax 


Pmax 


Pmax 


Pmax 


Pmax 


Pmax 


Pmax 


1.008±0.042 


1.07910.075 


1.06110.051 


0.91410.040 


1.20410.084 


1.50410.119 


1.26710.213 


0.48710. 


1.07710.045 


-1.11910.232 


-1.44210.414 


-1.47810.278 


-1.15510.230 


-2.73410.455 


-4.68310.686 


-3.29211.255 


1.46010.533 


-1.70410.268 


0.131 


0.130 


0.147 


0.166 


0.163 


0.138 


0.238 


0.126 


0.112 


Table  03.  Magnitude 


Standard  Least 


N 


218 


58 


18 


142 


195 


98 


30 


19 


48 


30 


magnitude  Regression  Results 


Squares.  Outlier  Retained 


Slope  Ir^ercept 


a(Y) 


0.120 


BSW 


NNZ 


Sahara 


Pmax 


Pmax 


Pmax 


Pmax 


Pmax 


Pmax 


Pmax 


Pmax 


Pmax 


Pmax 


1.08710.021 


1.04010.037 


0.95710.139 


1.07710.027 


0.99610.006 


0.98010.017 


1.06410.031 


1.01110.038 


1.00810.014 


1.05510.014 


1.00110.113 


-0.81710.112 


-0.1 4610.708 


-0.179l0.043 


-0.09110. 


-0.58910.1 76 


-0.27410.225 


-0.245l0.087 


-0.502l0.079 


-0.300i0.595 


0.131 


0.044 


0.048 


0.046 


0.043 


0.025 


0.042 


0.113 


NTS 

Pmax 

Pa 

214 

0.94410.029 

-0.28910. 157 

0.165 

0.91 

PMA 

Pmax 

Pa 

58 

1.01910.056 

-0.675i0.311 

0.125 

0.92 

RNA 

Ptnax 

Pa 

15 

1.09710.281 

-1.08511 .440 

0.182 

0.73 

YFT 

Pmax 

Pa 

141 

0.90710.036 

-0.105i0.193 

0.177 

0.90 

KTS 

Pmax 

Pa 

195 

0.97310.012 

0.066 

0.99 

Oeg 

Pmax 

Pa 

98 

-0.50710. 151 

0.077 

0.96 

BNE 

Pmax 

Pa 

30 

1.02710.038 

-0.63110.218 

0.057 

0.98 

BTZ 

Pmax 

Pa 

19 

0.98810.055 

-0.42210.327 

0.062 

0.97 

BSW 

Pmax 

Pa 

48 

1.02010.021 

0.60410.125 

0.036 

009 

NNZ 

Pmax 

Pa 

30 

rnmmsm 

-0.59310.124 

0.066 

0.99 

Sahara 

Pmax 

Pa 

7 

0.74010.066 

0.94610.345 

0.065 

0.98 

NTS 

Pmax 

Pn 

221 

0.98610.044 

-0.14410.232 

0.256 

064 

PMA 

Pmax 

Pn 

54 

0.98210.043 

0.03910.235 

0.090 

0.95 

RNA 

Pmax 

Pn 

22 

0.80110.171 

0.48610.868 

0.145 

0.72 

YFT 

Pmax 

Pn 

145 

0.85710.051 

0.52610.269 

0.259 

061 

KTS 

Pmax 

Pn 

113 

0.94410.085 

0.20710.478 

0.73 

Deg 

Pmax 

Pn 

55 

0.81410.137 

0.87810.727 

0.315 

lEISI 

BNE 

Pmax 

Pn 

20 

0.68910.431 

1.46712.458 

0.542 

0.35 

BTZ 

Pmax 

Pn 

6 

3.54211 .340 

0.169 

0.67 

BSW 

Pmax 

Pn 

32 

0.72510.332 

1.63511.998 

0.358 

067 

NNZ 

Pmax 

Pn 

20 

■KEninilES 

0.97510.788 

0.248 

064 

NTS 

Pmax 

Lg 

183 

0.79910.045 

1.28310.239 

0.242 

060 

PMA 

Pmax 

Lg 

48 

0.85310.070 

1.02810.387 

0.135 

067 

RNA 

Pmax 

Lg 

22 

0.97110.176 

0.15410.891 

0.148 

0.78 

YFT 

Pmax 

Lg 

113 

0.70710.056 

1.80410.298 

0.259 

0.77 

KTS 

Pmax 

Lg 

74 

0.73410.051 

1.39810.300 

0.113 

066 

Table  03.  Magnitude  .magnitude  Regression  Results  (Cent) 


Standard  Least  Squares,  Outlier  Retained 


Intercept 


0.762±0.688 


1.856±0.724 


1.865±0. 


2.312±0.881 


1.935±0.232 


N 

Slope 

12 

0.83610.121 

20 

0.65710.125 

15 

0.^10.099 

27 

0.58410.146 

24 

0.71110.039 

140 

1.16710.058 

53 

1.09310.101 

84 

1.10510.064 

154 

0.93910.056 

63 

1.18510.137 

26 

1.23410.152 

18 

1.26710.213 

47 

0.87610.118 

29 

0.99910.066 

-1 .47810.561 


-1.71810.350 


-1.28910.318 


-2.535l0.732 


o(Y) 


0.111 


0.126 


.107 


0.106 


0.109 


0.220 


.192 


0.203 


0.269 


.311 


0.209 


-0.92610.704 


-1.25610.386 


0.200 


0.189 


APPENDIX  E:  SMALLEST  EVENTS 


The  following  Table  lists  those  events  with  the  smaller  f<%{Pmaid  in  each  source  regioa  Al  events 
with  n\{P^  (GLM)  i  4.80  are  also  induded. 


Table  E.  Smallest  GLM  Event(s)  of  Each  Test  Site 


Event 

Site 

Pnmx 

Pt 

Pm 

Pm 

^0 

US681208 

PMA 

4.69±0.07 

4.28±0.07 

3.9210.08 

4.5610.09 

5.1610.15 

1 

US731012 

RNA 

4.80±0.12 

_ ± _ 

_ ^1 _ 

4.1510.12 

4.8210.11 

_ ± _ 

US780913 

RNA 

4.69±0.17 

_ 

1 

4.5710.14 

4.8310.15 

1 

US620706 

YFT 

4.73±0.17 

4.42±0.17 

3.9010.17 

4.3910.18 

4.8610.21 

4.2610.09 

US640716 

YFT 

4.50±0.i2 

4.10±0.12 

4.3310.12 

4.2510.10 

4.9610.15 

1 

US641009 

YR 

4.76±0.07 

4.51±0.07 

4.5210.10 

4.5710.10 

5.1310.21 

1 

US650218 

YFT 

4.47±0.12 

4.27±0.12 

4.4210.12 

4.1410.10 

4.9510.15 

1 

US650514 

YFT 

4.73±0.17 

3.97±0.17 

1 

4.4410.10 

4.7110.15 

_ 1 _ 

US650521 

YFT 

4.69±0.08 

4.01±0.08 

4.1110.10 

4.4410.09 

5.1410.15 

1 

US660625 

YFT 

4.61±0.07 

4.07±0.08 

3.8310.10 

4.3410.11 

5.0410.21 

1 

US670727 

YFT 

4.78±0.10 

4.19±0.10 

4.4710.10 

4.3610.12 

4.7910.15 

_ ^1 _ 

US680119 

YFT 

4.48±0.17 

4.4010.17 

4.5310.18 

1 

1 

US680325 

YFT 

4.45±0.12 

± 

1 

4.7310.15 

1 

US680410 

YFT 

4.78±0.12 

3.96±0.12 

4.0110.12 

4.3710.10 

_ 1 _ 

1 

US690716a 

YFT 

4.40±0.12 

± 

1 

1 

1 

US691121 

YFT 

4.79+0.08 

4.46±0.10 

4.4210.10 

4.7710.12 

5.3010.21 

1 

US700205 

YFT 

4.72±0.08 

4.33±0.08 

4.0210.10 

Mffgy.ytai 

4.9810.15 

1 

US720519 

YFT 

4.65±0.10 

± 

1 

4.4410.25 

5.0610.15 

1 

US750424 

YFT 

4.57±0.10 

3.9e±0.10 

3.6910.10 

4.5110.10 

4.8310.15 

1 

US790803 

YFT 

4.74±0.06 

4.40±0.06 

4.2110.08 

4.7010.11 

5.2310.12 

1 

US790808 

YFT 

4.76±0.06 

4.41±0.07 

4.4010.06 

4.5910.09 

5.2910.11 

1 

US820729 

YFT 

4.54±0.10 

3.87±0.12 

4.0310.12 

4.4410.07 

4.9010.09 

3.3210.07 

US840131 

YFT 

4.44±0.12 

4.0810.12 

4.3710.17 

4.2410.08 

4.9310.11 

1 

US840620 

YFT 

4.74±0.07 

4.2810.07 

4.2310.07 

4.3810.11 

^1 

3.4510.05 

SH890902 

BNE 

5.01±0.06 

4.8110.06 

4.6210.06 

5.6710.18 

1 

"  3."(^.05 

SH760421 

BSW 

5.13±0.12 

4.9710.12 

4.7210.12 

1 

1 

2.9310.17 

SH710630 

BTZ 

- : - 1 

5.2110.04 

-i: - 1 

5.0010.04 

4.7610.05 

5.6810.25 

5.3510.12 

3.3210.12 

DM890217 

Deg 

4.77±0.17 

4.6310.12 

4.0110.10 

5.1310.14 

+ 

1 

NZ641025 

iKlB 

4.67±0.04 

4.4710.04 

4.2810.04 

5.1910.11 

5.3410.11 

1 

NZ771009 

NNZ 

4.2610.04 

5.1010.12 

1 

3S210.10 

FS640214 

Sahara 

4.42±0.17 

1 

± 

1 

1 

1 

FS641128 

Sahara 

4.4310.12 

4.4510.12 

1 

1 

1 

FS651001 

Sahara 

4.51±0.17 

1 

1 

1 

1 

1 
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